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The  impedance  of  an  antenna  probe  surrounded  by  ionospheric  plasma 
is  examined.  A  review  of  several  physical  plasma  parameters,  research 
techniques,  some  past  experimental  results,  and  various  theories  is  pre¬ 
sented.  A  need  is  shown  for  further  theoretical  sophistication.  The 
theory  develops  the  impedance  of  an  electrically  short,  cylindrical  an¬ 
tenna  probe  immersed  in  a  warm,  lossy,  compressible,  magneto- ionic, 
electron  fluid.  The  results  are  examined  and  compared  to  previous  work 
and  are  appLied  to  several  ionospheric  experiments.  Experimental  re¬ 
sults  as  interpreted  by  the  theory  are  examined  and  the  usefulness  of 
the  theory  is  thus  demonstrated.  In  particular,  the  theory  is  shown  to 
be  of  value  in  determining  electron  density,  temperature,  and  collision 
frequency  from  RF  probe  measurements. 
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Chapter  I 


INTRODUCTION 

It  is  the  purpose  of  this  report  to  examine  the  interaction  between 
an  antenna  probe  and  the  ionospheric  plasma  surrounding  it.  This  intro¬ 
ductory  chapter  discusses  plasma  and  ionospheric  parameters  and  several 
research  techniques.  Chapter  II  considers  past  theoretical  work  that 
relates  antenna  probe  impedance  to  the  plasma  parameters.  The  need  for 
further  effort  is  then  demonstrated.  A  theory  that  relates  the  imped¬ 
ance  of  an  electrically  short,  cylindrical,  antenna  probe  to  the  magnetic 
field  and  plasma  parameters  of  a  warm,  lossy,  compressible  electron  fluid 
is  developed  in  Chapter  III.  The  theory  is  explored  as  a  function  of  its 
parameters  and  is  compared  to  other  accepted  results  in  Chapter  IV. 

Chapter  V  describes  the  application  of  the  work  to  various  experiments. 

In  particular,  the  theory  is  useful  in  interpreting  measurements  of  electron 
density,  electron  temperature,  and  electron  collision  frequency.  The  re¬ 
sults  of  the  analysis  then  demonstrate  the  usefulness  of  the  theory  and 
its  broad  application  in  ionospheric  research. 

It  should  be  noted  that  rationalized  MKS  units  are  used  throughout. 

The  Plasma  State 

A  gaseous  plasma  is  an  electrically  neutral  system  of  particles  con¬ 
taining  primarily  free  electrons  and  positive  ions,  but  sometimes  negative 
ions  and  neutral  molecules  as  well.  It  has  been  estimated  that  perhaps 

more  than  99.9  percent  of  the  universe  exists  in  such  an  ionized  form 

-2 

[Baohynskij  1961],  While  approximately  10  electron  volts  of  energy 


2 


per  particle  are  required  to  change  from  the  solid  to  liquid  or  from  the 
liquid  to  gaseous  state,  one  to  ten  electron  volts  of  energy  per  particle 
are  required  to  obtain  the  plasma  state.  Coulomb  forces  are  dominant  in 
plasmas  as  they  are  much  more  effective  than  gravitational  or  short-range 
nuclear  effects.  Thus,  the  ionized  particles  are  endowed  with  considerable 
energy  and  exhibit  strong  electrical  interactions  that  determine  the  pri¬ 
mary  character  of  the  plasma. 

A  system  of  ionized  particles  is  considered  to  be  a  plasma  only  if 
collective  effects  dominate  any  single  particle  effects.  Therefore,  the 
dimensions  of  a  plasma  must  be  at  least  as  large  as  the  range  of  inter¬ 
action  of  single  particles.  This  range  is  called  the  Debye  length  A^ 

[ Brown  et  al.j  1963]  and  can  be  expressed  by  the  equation 


A 


D 


£o  KT 


1/2 


(1) 


Where 

K  =  Boltmann's  Constant 
T  =  kinetic  temperature  (electrons) 
N  =  ionization  density  (electrons) 
e  *  electronic  charge 
£  =  permittivity  of  free  space 


The  Debye  length  is  a  very  useful  parameter  and  is  often  used  as  a  length 
scaling  constant  in  normalized  equations.  The  Debye  length  can  also  be 
thought  of  as  the  maximum  distance  over  which  some  charge  separation  can 
occur.  If  any  macroscopic  region  (large  dimensions  compared  to  A^) ,  free 
of  external  influence,  should  have  a  net  charge,  then  the  charge  would 
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create  an  electric  field.  This  field  would  then  cause  the  free  ionized 
particles  to  redistribute  into  a  minimum  energy  condition  of  zero  field 
and  thus  zero  net  charge  density  would  again  result  It  is  not  surpris¬ 
ing  then  to  find  that  the  characteristic  shielding  distance  a  is  directly 
dependent  upon  the  energy  of  the  electrons  (electron  temperature)  and  is 
inversely  related  to  the  average  number  of  electrons  that  shield  any  dis¬ 
continuities. 

From  the  above  discussion,  it  is  apparent  that  disturbed  plasma  tends 
to  return  to  equilibrium.  Because  the  particles  also  have  mass,  the  sys¬ 
tem  exhibits  properties  characteristic  of  second  order  systems,  including 
a  natural  resonant  frequency,  known  as  the  plasma  frequency,  f^.  This 
quantity,  the  frequency  of  the  "quivers"  of  the  "plasma  jelly,"  can  be 
directly  related  to  the  electron  density  N  and  the  physical  constants 
electronic  charge  e,  electronic  mass  m,  and  permittivity  of  free  space 
£  by  the  equation 


[Brown  et  at  3  1963].  The  only  variable  in  this  equation  is  N,  the 
electron  density.  Since  the  feature  that  distinguishes  a  plasma  from 
other  types  ot  matter  is  its  ionization,  it  is  natural  to  expect  that 
the  electron  density  is  the  major  parameter  that  describes  the  behavior 
of  a  plasma. 

At  this  point,  perhaps  it  is  only  natural  to  ask,  "Why  are  the 
variables,  electron  density  and  temperature,  the  only  ones  specified  in 
Ap  and  u)N?  Shouldn't  the  positive  ion  density  and  temperature  also  be 
included?"  Of  course  in  the  strict  sense,  the  answer  is  that  and 
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are  defined  in  this  manner,  but  it  is  true  that  shielding  and  naturaL 
oscillation  in  plasmas  are  indeed  affected  by  other  types  of  particles. 
However,  the  electrons  are  much  less  massive  than  the  other  particles 
(by  about  four  orders  of  magnitude)  and  since  they  have  the  same  charge, 
they  respond  to  electric  and  magnetic  fields  to  a  much  larger  degree. 
Then  coo,  approximate  thermal  equilibrium  usually  exists  in  a  plasma; 
therefore,  electrons  have  the  same  energy  per  particle  as  the  other 
heavier  particles.  Thus  the  average  electron  velocity  is  at  least  two 
orders  of  magnitude  larger  than  the  ion  velocity.  This  is  easily  seen 
by  equating  the  mean  kinetic  energies  as  follows 


1/2  mv2  *  1/2  Mv2 


v  102  v 


where 


(3) 

(4) 


l5> 


v  =  mean  electron  velocity 

v  =  mean  molecular  velocity 

m 

M  =  molecular  mass 
m  =  electronic  mass 

Thus  the  electrons  can  be  visualized  as  mobile,  rapidly  moving  particles 
among  relatively  fixed  positive  ions  and  neutral  molecules.  For  most  of 
our  work,  it  will  be  possible  to  consider  a  plasma  as  compressible  elect¬ 
ron  gas,  diffused  among  fixed  ions  and  molecules  whose  only  effect  is  to 
maintain  average  charge  neutrality  and  to  possibly  thermalize  the  rapidly 
moving  electrons  through  collisions.  This  will  require  that  any  exciting 
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frequency  a  be  at  least  as  large  as  10  to  10  a  .  Otherwise  ion 

•  *  A 

and  molecular  motions  must  also  be  considered. 

A  useful  measure  of  the  collisions  of  electrons  with  other  paricles 
is  the  effective  electron  collision  frequency  j .  It  is  defined  as  the 
average  number  of  times  a  second  that  an  electron  suffers  a  collision 
that  effectively  changes  its  direction  of  motion  [ Carribel ,  1963] .  Pfister 
(196&J  considers  several  definitions  of  electron  collision  frequency  and 
relates  these  to  the  mean  collision  frequency  as  derived  from  the  average 
electron  velocity. 

If  the  exciting  frequency  a  is  considerably  larger  than  then 
essentially  no  collisions  occur  during  one  cycle  of  excitation.  Thus 
only  one  degree  of  freedom  is  possible  for  the  electrons,  that  is,  the 
direction  of  the  original  excitation.  Thus,  the  ratio  of  specific  heats 
for  our  compressible  electron  gas  must  be  taken  as  3  \3pitzer-,  1962*. 

It  was  shown  earlier  that  steady  state  charge  separation  does  not 
generally  exist  in  a  plasma  However,  one  important  exception  occurs 
at  any  sharp  boundary  where  the  ionization  density  changes  markedly 
over  a  few  Debye  lengths.  At  such  a  surface,  in  the  absence  of  electro¬ 
magnetic  fields,  the  electron  flux  impinging  upon  the  boundary  will 
greatly  exceed  the  positive  ion  flux  because  of  the  much  larger  average 
electron  velocity.  Thus  in  the  steady  stare,  electric  fields  are  set 
up  by  the  accumulation  of  electrons  on  the  surface,  such  that  the  ion 
and  electron  currents  are  equal-  There  results,  therefore,  a  depletion 
of  electrons  near  the  boundary  and  a  sheath  (free  of  electrons)  is  said 
to  'ixist  between  the  surface  and  the  plasma.  The  thickness  of  this 
sheath  is  on  the  order  of  a  Debye  length  and  can  be  modified,  even  dis¬ 
pelled,  by  the  external  application  of  electric  potentials  to  the  surface. 
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This  generally  draws  current  from  the  plasma  and  thus  upsets,  to  some 
degree,  the  natural  plasma  conditions.  Nevertheless,  it  will  sometimes 
be  assumed  that  sheaths  around  antenna  probes  can  be  collapsed  so  that 
sheath  effects  can  be  ignored. 

An  external  magnetic  field  greatly  modifies  the  behavior  of  a  plasma. 
The  general  effect  is  to  "stiffen"  the  plasma  in  directions  normal  to  the 
direction  of  the  field.  Magnetic  fields  are  thus  responsible  for  the  an¬ 
isotropic  character  of  such  plasmas.  It  is  possible  to  consider  the  field 
as  a  parameter  of  the  medium,  however,  a  related  quantity,  the  gyrofreq- 

uency  f  ,  will  often  be  more  useful.  Ihe  relation  is 
h 


CO 


H 


2”£h 


eB 

m 


(6) 


where 


B  *  magnetic  flux  density. 

In  summary,  our  physical  model  of  the  plasma  state  may  be  visualized 
as  a  warm,  compressible  gas  of  electrons,  diffused  among  fixed  ions  and 
molecules  so  that  macroscopic  charge  neutrality  is  maintained.  The  proper¬ 
ties  of  the  medium  can  be  expressed  by  specifying  the  plasma  parameters, 
that  is,  by  stating  plasma  frequency  f  ,  gyrof requency  fu,  electron  col- 
lisional  frequency  v,  electron  temperature  T  and  the  ratio  of  specific 
heats  of  the  electron  gas. 

The  Ionosphere 

The  ionosphere  is  a  particularly  significant  system  which  is  com¬ 
posed  of  dilute,  weakly  ionized  plasma  in  the  earth’s  magnetic  field. 

It  is  produced  above  the  earth,  primarily  by  the  ionizing  radiation  from 
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the  sun  incident  upon  the  earth's  upper  atmosphere  The  structure  of  the 
ionosphere  is  thus  strongly  dependent  upon  the  local  sun  angle  and  the 
level  of  solar  activity  [Johnson,  1961],  It  is  customary  to  identify 
several  distinctive  levels  of  ionization  as  the  D,  E  and  F  layers  as 
shown  in  Figure  1, 

The  parameters  of  the  ionospheric  plasma  are  illustrated  in  Figures 
1,  2,  3  and  4.  Figure  1  shows  electron  density  and  plasma  frequency  as 
a  function  of  altitude.  Figure  2  is  a  similar  plot  of  electron  temperature, 
while  the  expected  range  of  Debye  length  is  shown  in  Figure  3.  The  expected 
electron  collision  frequency  is  given  in  Figure  4.  The  gyrof requency  in  the 
ionosphere  is  completely  determined  by  the  strength  of  the  earth's  magnetic 
field.  This  does  not  change  radically  as  a  function  of  height  in  this  re¬ 
gion,  arid  the  value  1,4  MHz  is  a  representative  value. 

Since  the  time  (1902)  Ileavyside  and  Kennelly  [ Mitra ,  1952,  p.  177] 
first  postulated  the  ionosphere  to  explain  transatlantic  radio  communi¬ 
cation,  measurements  of  the  properties  of  the  ionosphere  have  been  of 
great  interest.  Not  only  are  the  ionospheric  parameters  of  interest  for 
communication  studies,  but  they  are  also  important  to  studies  of  other 
geophysical  phenomena  such  as  aurora,  airglow,  radio  noise,  magnetic  field 
fluctuations  and  weather.  More  recently,  perturbation  of  the  ionosphere 
by  space  vehicles  and  atomic  weapons  tests  have  emphasized  the  need  for 
studies  of  localized  regions  of  the  ionosphere.  Thus  it  can  be  seen  that 
a  knowledge  of  the  properties  of  the  ionosphere  is  very  important  for 
many  activities. 


ALTITUDE  ( km)  ALTITUDE  (km) 


PLASMA  FREQUENCY  (MHz) 


DAYTIME  ELECTRON  CONCENTRATION  (•Isclrorw/m8) 


PLASMA  FREQUENCY  (MHz) 


NIGHTTIME  ELECTRON  CONCENTRATION  UI«ctraA«/m9) 


Fig.  1.  Normal  electron  distribution  of  the  ionosphere* 
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Fig.  2.  Normal  electron  temperature  Fig.  3.  Normal  Debye  length  in 
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Fig.  4.  Normal  electron  collision 
frequency  in  the  ionosphere. 

Ionospheric  Measurement  Techniques 

Ionospheric  parameters  have  been  traditionally  measured  from  the 
earth's  surface  by  several  techniques,  notably  by  ionospheric  sounders 
[Mitra,  1952,  p.  219].  However,  fine  structure  of  localized  regions  and 
disturbed  areas  of  the  ionosphere  can  only  be  examined  by  direct,  in  sitv. 
measurements.  Sounding  rockets  are  often  the  most  useful  for  examining 
these  localized  or  disturbed  regions.  Satellite  vehicles  are,  on  the 
other  hand,  especially  useful  in  obtaining  worldwide  coverage,  often 
over  an  extended  period  of  time. 

Three  general  plasma  measurement  techniques  have  been  extensively 
exploited  for  rocket  and  satellite  application.  The  first  technique 
was  developed  from  the  concepts  of  the  ground-based  ionospheric  sounder. 
That  is,  the  electromagnetic  propagation  characteristics  of  the  ionos¬ 
pheric  plasma  are  measured  by  transmitting  probing  signals  between  the 
carrier  vehicle  and  the  ground.  Many  variations  of  this  basic  technique 
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have  been  utilized  even  though  extensive  instrumentation  and  complex 
analyses  are  required. 

The  second  technique  is  the  application  of  the  Langmuir  Probe  [Lang¬ 
muir,  1923]  to  ionospheric  measurements.  One  or  more  electrodes  ate  ex¬ 
tended  from  the  vehicle  into  the  ionospheric  plasma  and  excited  by  various 
potentials.  The  resulting  dc  currents  carry  the  measurement  information. 
Related  devices  such  a c  ion  traps  and  resonance  probes  have  also  been 
developed.  These  all  ise  dc  potentials  to  collect  and  examine  the  ionized 
plasma  particles  by  measuring  dc  currents  from  the  collector  electrodes. 
These  techniques  have  nearly  dominated  the  measurements  of  electron  and 
ion  temperatures  and  ion  densities.  These  methods  often  suffer  from  dis¬ 
turbing  vehicle  potential  perturbations;  nevertheless,  they  have  provided 
the  only  direct  measurements  of  electron  temperature. 

The  third  method  is  the  measurement  of  the  plasma  perturbated,  RF 
probe  impedance.  Although  previously  well  known  in  the  laboratory,  iono¬ 
sphere  measurements  using  this  technique  seem  to  have  started  when  ap¬ 
parent  antenna  impedance  changes  were  noted  in  a  dipole  being  used  in  a 
rocket  propagation  experiment.  This  effect  was  then  used  to  estimate 
the  local  electron  density.  Since  that  time,  many  variations  of  the 
basic  technique  have  been  used.  The  interaction  of  an  antenna  with  the 
ionospheric  plasma  is  the  prime  concern  here,  thus  we  will  examine  some 
of  the?'  IF  impedance  type  experiments  in  details.  Later,  flight  data  from 
some  of  these  experiments  will  be  analyzed  by  applying  the  theory  developed 
he rein. 

Three  RF  probes  will  be  considered.  Each  technique  utilizes  a  pair 
of  2.5  to  3-meter  probes,  sometimes  referred  to  as  a  dipole  antenna,  moun¬ 
ted  normal  to  the  longitudinal  vehicle  axis  and  extended  into  the  iono- 
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Fig.  5. 


spheric  plasma  (see  Figure  o).  In  each  case,  the  antenna  is  excited  with 
a  low-level  RF  signal  (approximately  0.5  volt  rms)  in  the  0.5  to  15-MHz 
range;  in  addition,  a  slowly  varying  bias  voltage  from  0  to  +5  volts  is 
placed  on  the  antennas  with  respect  to  the  vehicle  skin.  The  bias  modi¬ 
fies  the  sheath  and  allows  the  effect  of  the  sheath  to  be  determined. 

An  analysis  of  the  effects  of  the  plasma  on  the  antenna  RF  impedance  or 
on  the  dc  current  variations  as  a  function  of  frequency  provides  a  mea¬ 
sure  of  the  plasma  parameters. 

Standing  Wave  Impedance  Probe  (SWIP) 

Perhaps  one  of  the  most  successful  impedance  probes  is  the  Standing 
Wave  Impedance  Probe  [Ulwick  et  al. }  1964].  This  system  measures  the  RF 
impedance  of  a  dipole  antenna  which  is  being  driven  by  one  or  two  time- 
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multiplexed  RF  frequencies  (see  Figure  6).  The  impedance  probe  frequencies 
are  set  so  that  one  or  the  other  is  near  the  local  plasma  frequency  through¬ 
out  the  experiment.  The  RF  fields  from  the  antenna  excite  motion  in  the 
plasma  electrons  which,  in  turn,  modify  the  antenna  impedance.  This  change 
of  impedance  is  indicated  by  a  change  in  the  standing  wave  pattern  of  the 
voltage  along  an  artificial  transmission  line  connecting  the  driving  RF 
oscillator  and  the  antenna.  The  voltage  standing  wave  pattern  is  sampled, 
detected  and  telemetered  to  a  ground  station  where  it  is  recorded  for 
further  analysis.  This  experiment  can  provide  reliable  impedance  measure¬ 
ments  at  one  or  two  fixed  exciting  frequencies. 

A  very  simple  analysis  is  usually  employed  [Ulwick  et  al. ,  1964],  and 
this  generally  produces  good  results  [Baker,  et  al. ,  1966].  Occasionally, 
electron  temperature  and  other  effects  become  important,  thus  pointing  out 
the  need  for  further  theoretical  development  [Baker  et  al. ,  1965,  and 
Whale ,  1962]. 

Plasma  Frequency  Probe  (PFP) 

The  Plasma  Frequency  Probe  [Haycock  and  Baker,  1962]  excites  the  anten¬ 
na  with  a  variable  frequency  that  sweeps  from  about  .1  to  10  MHz  (see 
Figure  7).  The  resonance  frequencies  of  the  antenna  impedance  are  sensed 
by  determining  the  frequency  where  the  phase  angle  between  the  RF  antenna 
current  and  voltage  is  zero.  A  parallel  resonance  condition  (high  imped¬ 
ance)  is  expected  to  occur  near  the  natural  frequency  of  the  plasma,  that 
is,  the  plasma  frequency.  Experimentally,  it  has  been  found  that  a  series 
resonance  condition  (low  impedance)  also  occurs  [Baker  et  al. ,  1966].  It 
will  be  shown  later  that  this  series  resonance  is  related  to  the  electron 
temperature  of  the  plasma. 


,  'Wrr  , 


Standing  Wave  Impedance  Probe. 


Fig.  7.  Plasma  Frequency  and  Resonance  Rectification  Probe. 
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Resonance  Rectification  Probe  (REP) 

The  Resonance  Rectification  Probe  [ Takayama 3  Ikegami  and  Miyazaki, 
1960;  Despain,  1964]  is  a  combination  of  a  Langmuir  Probe  and  an  RF 
sweep  oscillator  similar  to  that  used  in  the  Plasma  Frequency  Probe  (see 
Figure  7) .  A  dc  bias  (0  to  +5  volts) ,  stepped  with  successive  RF  oscil¬ 
lator  sweeps,  is  applied  to  the  antenna.  The  resulting  dc  current  is 
monitored.  Because  the  RF  antenna  impedance  changes  with  frequency  and 
the  dc  characteristics  of  the  antenna  probe  exhibit  nonlinearities,  the 
dc  current  is  a  function  of  the  RF  frequency.  It  had  been  contended 
that  the  frequency  where  the  dc  current  reaches  a  maximum  is  the  local 
electron  plasma  frequency  [ Takayama ,  1960;  Miyazaki  et  al 1960;  Ichi¬ 
kawa  and  Ikegami ,  1962;  Hirao  and  Muraoka,  1964];  however.  Baker,  Des- 
pain  and  Ulwick  [1966]  have  shown  experimentally  that  this  conclusion 
is  not  true;  in  fact,  the  maximum  is  more  closely  associated  with  the 
series  resonance  frequency.  Thus  additional  theoretical  work  on  antenna 
impedance  is  also  needed  to  aid  in  the  interpretation  of  this  experiment. 

Preliminary  Conclusion 

The  ionosphere  is  of  general  interest  and  exhibits  interesting  and 
complex  behavior.  A  knowledge  of  the  ionospheric  parameters,  especially 
in  localized  regions,  is  prerequisite  to  an  understanding  of  this  be¬ 
havior.  Some  of  the  probe  techniques  for  ionospheric  parameter  measure¬ 
ment  require  an  analysis  of  the  interaction  of  an  RF-driven  probe  with 
the  surrounding  ionospheric  plasma.  The  theories  of  such  interactions 
will  now  be  discussed  in  Chapter  II. 


Chapter  II 


PROBE  THEORIES 

The  problem  of  interaction  of  a  small  RF  excited  probe  with  a  plasma 
has  been  studied  for  at  least  sixty  years  [Wilson  and  Gold t  1906].  Most 
investigators  have  assumed  that  the  impedance  change  is  due  solely  to  a 
change  in  the  effective  complex  permittivity  e  of  a  locally  homogeneous 
ionized  gas  as  expressed  in  the  Appleton-Hartree  Theory  [Ratcliff 1962]; 
that  is,  the  electrons  are  considered  as  mobile  charges,  free  to  move  un¬ 
der  the  influence  of  electromagnetic  fields.  The  electric  current  that 
results  from  such  movement  is  accounted  for  by  a  change  in  the  permittivity 
of  the  medium.  If  thermal  motions  and  magnetic  fields  are  ignored,  the 
result  of  such  a  calculation  is  shown  in  Appendix  A  to  produce  the  following 
expression  for  the  total  antenna  impedance  Z^. 


where 

Cq  =  free  space  capacitance  of  the  probe 
u)  =  radian  operating  frequency 
=  radian  plasma  frequency 
v  =  electron  collision  frequency 

Rr  =  free  space  radiation  resistance 
X  .  uj/u 

Z  =  v/u) 
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For  and  v  equal  to  zero,  a  particularly  simple  interpretation  results 

in  the  probe  being  considered  as  a  capacitor  C  of  value 


This  simple,  idealized  model  has  been  extensively  applied  to  interpret  im¬ 
pedance  probe  results  [Ulwiak  et  al.3  1964;  Jackson  and  Pickar3  1957;  Whale 
1963]. 

A  more  rigorous  approach  by  King  and  Harrison  [I960]  results  in  some 
highly  involved  expressions;  while  an  expression  developed  by  DeChamps 
[1961]  is  somewhat  simpler.  However,  neither  of  these  works  consider 
magnetic  field  or  electron  temperature  effects  and  hence  add  little  to  the 
general  formulation  that  we  seek  here.  If  magnetic  field  effects  are  to 
be  considered,  then  the  dielectric  constant  becomes  a  tensor  [ Ratcliff e 3 
1962]  and  the  geometry  of  the  probe  must  be  considered.  The  prime  interest 
of  this  paper  is  cylindrical  geometry;  therefore,  only  solutions  for  this 
case  are  considered. 

Perhaps  the  most  obvious  approach  for  a  tensor  dielectric  constant  is 
to  apply  the  tensor  [e]  as  given  by  the  usual  Appleton-Hurtree  formu¬ 
lation  [Ratcliffe3  1962]  or  the  formulation  of  Sen  and  Wyler  [I960]. 

The  Sen-Wyler  approach  allows  for  collisional  frequency  variations  with 
electron  energy  but  is  otherwise  identical  to  the  Appleton-Hartree  tensor. 
Shkarofksy  [1961]  has  generalized  the  Appleton-Hartree  equations  even  further 
by  considering  the  degree  of  plasma  ionization  as  well. 

The  question  of  how  a  given  [e]  should  be  applied  to  determine  an¬ 
tenna  impedance  has  been  considered  by . several  workers.  Herman  [1963] 
and  Crouse  [1964]  investigated  conical  geometries,  while  Kaiser  [1962] 
considered  biconical  antennas .  These  workers  all  assumed  a  simple  capa¬ 
citor  model.  In  Appendix  A,  the  impedance  of  a  cylindrical  monopole, 
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using  the  capacitor  model  and  the  Appleton-Hartree-  formulation  for  [e] 
is  also  developed.  The  result  is 


Z 
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jwC 


1  +  X 


1  -  1/2  (Y/U)' 


Sin2y 


U((Y/U)  -1)  J 


-1 


(9) 


where 

U  =  1  -  jZ 

Y  = 

Y  =  Angle  between  antenna  axis  and  the  magnetic  field 


The  capacitor  models  discussed  above  ignore  source  current  distri¬ 
butions  and,  in  fact,  all  magnetic  fields  except  the  steady  field  as  it 
appears  in  [e].  The  classical  approach  to  antenna  theory,  on  the  other 
hand,  generally  assumes  a  source  current  distribution  and  all  the  para¬ 
meters  are  then  determined  from  this  specification.  Bromley  [1962]  used 
this  approach  to  first  obtain  results  in  an  isotropic  region  similar  to 
those  of  King ,  Harrison  and  Denton  [1961].  Then  he  extended  the  results 
to  the  magneto-ionic  case  by  deriving  an  approximation  valid  for  weak 
magnetic  fields  or  small  electron  densities.  Ament  et  at.  [1964]  con¬ 
sidered  the  general  case  and  produced  an  expression  for  impedance  that 
requires  extensive  numerical  in*  ‘•gration  and,  in  general,  seems  to  be 
very  difficult  to  apply,  Balmain  [1964],  however,  has  developed  some 
very  useful  impedance  relations.  His  results  appear  to  be  the  most  phy¬ 
sically  satisfying  of  all  the  work  mentioned  above  and  yet  are  quite 
tractable  in  terms  of  computation.  His  expression  is 


Z  = 


jw2ire  LK'/f 
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where 
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None  of  the  above  derivations  includes  the  explicit  effects  of 
finite  electron  temperature.  That  is,  a  cold,  yet  lossy  medium  is  always 
assumed.  However,  electron  temperature  can  appreciably  perturbate  the 
impedance  of  an  antenna  probe  in  a  plasma,  other  than  its  effect  on  col- 
lisional  frequency.  Whale  [1963]  observed  such  effects  in  an  impedance 
probe  experiment,  primarily  as  an  energy  loss  mechanism.  It  will  also 
be  shown  later  that  the  experimentally  determined  antenna  series  resonance 
frequencies  of  Baker t  Despain  and  Ulwick  [1966]  are  indeed  related  to 
the  electron  temperature  as  they  suspected.  Thus  electron  temperature  must 
be  considered  in  a  derivation  of  antenna  impedance. 

Several  workers  have  attempted  to  include  electron  temperature  as  a 
parameter.  Cohen  [1961]  presented  a  thorough  theoretical  discussion  of 
sources  in  isotropic  warm  plasmas.  Whale  [ 1963 [  developed  some  theory 
to  explain  his  experimental  results  mentioned  above.  Balmain  [1964] 
applied  Cohen's  equations  to  his  problem  for  the  case  of  no  magnetic  field 
and  derived  an  approximate  impedance  relation.  In  a  similar  manner,  Fejer 
[1964]  theoretically  investigated  spherical  geometry  without  magnetic  field 
and  predicted  a  series  resonance  below  the  plasma  frequency  directly  depen¬ 
dent  upon  the  electron  temperature,  Despain  [1965],  in  preliminary  work 
to  this  report,  investigated  cylindrical  geometry  with  similar  results. 

The  work  of  Wait  [1964]  should  also  be  mentioned  as  it  concerns  slot 
antennas  in  warm  plasmas.  An  attempt  to  develop  a  permittivity  tensor  [e] 
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valid  for  finite  electron  temperature  and  to  apply  this  [ e ]  to  the  capa¬ 
citor  model  of  the  antenna  met  with  failure  (see  last  section  of  Appendix  A) . 
The  impedance  expression  that  resulted  did  not  agree  even  approximately 
with  the  other  theoretical  work  nor  did  it  exhibit  behavior  known  experi¬ 
mentally  to  be  characteristic  of  such  systems.  Further  discussion  is  included 
in  Appendix  A. 

In  reviewing  the  work  mentioned  above,  it  is  apparent  that  further 
theoretical  work  is  needed  on  the  electron  temperature  effects  on  cylin¬ 
drical  antennas,  especially  in  the  presence  of  magnetic  fields  and  elec¬ 
tron  collisions.  Therefore,  the  present  investigation  was  initiated  in 
order  to  relate  electron  density,  temperature  and  collisional  frequency 
to  antenna  impedance. 


Chapter  III 


DERIVATION  OF  ANTENNA  IMPEDANCE 


Introduction 

The  basic  assumption  of  this  derivation  is  that  the  electrons  behave  as 
a  continuous  fluid  and  the  effect  of  all  electron  interactions  may  be  rep¬ 
resented  by  various  forces  on  the  electron  fluid.  This  view  of  a  plasma  is 
similar  to  that  of  Spitzer  [ 1962] ,  Cohen  [1961]  and  Fejer  [1964].  This  approach 
allows  the  electron  temperature  effects  to  be  taken  into  account  as  a  property 
of  the  electron  fluid.  Besides  the  effect  of  electron  density  and  temperature, 
the  effects  of  an  external  constant  magnetic  field  and  electron  collisions  are 
also  included.  It  is  further  assumed  that  the  undisturbed  plasma  is  neutral 
and  homogeneous.  Later  some  restricted  space  variations  of  the  plasma  para¬ 
meters  are  allowed  near  the  antenna  probe  surface  by  using  a  simple  sheath 
model . 

The  antenna  probe  is  taken  to  be  an  isolated,  electrically  short  cylinder, 
whose  length  to  radius  ratio  is  large  enough  so  that  end  effects  can  be  ignored. 
Because  of  the  long  wavelength  of  the  excitation,  electromagnetic  effects  are 
assumed  to  be  small  in  comparison  to  the  plasma  effects,  and  thus  the  antenna 
impedance  may  be  calculated  from  the  effective  quasi-static  electric  fields  of 
the  antenna.  It  is  further  assumed  that  the  excitation  level  is  very  low  so 
that  linearized  equations  can  be  used  and  only  first  order  effects  considered. 
Exp  (jut)  time  dependence  and  MKS  units  are  assumed  throughout. 

The  strategy  of  the  derivation  is  to  first  develop  a  force  equation  that 
relates  all  the  forces  on  an  element  of  the  electron  fluid.  From  this  equa¬ 
tion,  the  velocity  of  the  electron  fluid  is  found.  The  velocity  expression  is 
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then  substituted  into  the  continuity  equation  and  a  second  order  partial  dif¬ 
ferential  equation  is  derived  in  terms  of  the  scalar  parameters  n  (electron 
density  variation)  and  $  (electric  potential).  Poisson’s  equation,  also  a 
second  order  partial  differential  equation  in  4>  and  n,  is  used  to  simplify 
the  first  equation.  A  solution  of  this  equation  for  n  leads  directly  to  a 
solution  of  Poisson's  equation.  Boundary  conditions  can  be  applied  to  de¬ 
termine  a  solution  for  the  potential  <£.  Once  4>  is  known,  the  fields  sur¬ 
rounding  the  antenna  can  be  easily  determined,  and  the  antenna  impedance 
can  be  found  as  the  complex  ratio  of  antenna  potential  to  antenna  current. 
Sheath  effects  are  then  added  as  an  additional  effect. 

Fundamental  Equations 

Consider  the  forces  that  are  exerted  on  a  volume  element  dV  of  the 
electron  fluid.  The  force  Fp  due  to  the  variations  of  the  pressure  p  is 

dFp  =  -(Vp)dV  (11) 

The  Lorentz  force  F^  on  the  element  is 

dF  ■  -Ne[E  +  v  x  B]  dV  (12) 

L 

where 

N  =  average  electron  density  (constant) 
e  =  electronic  charge 

v  =  average  velocity  of  the  electrons  in  the  element 
B  =  external  constant  magnetic  field 

E  *  electric  field 

p  =  pressure  variation  from  its  mean  value 
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and  the  force  due  to  electron  collisions  is  giver  in  terms  of  the  effective 
electron  collision  frequency  v  as  follows 

dF  »  -vNm  v  dV  (13) 

c 

where  m  is  the  electronic  mass. 

Applying  Newton's  Law,  we  equate  these  forces  to  the  mass  of  the  element 
times  its  acceleration  to  produce  the  force  equation 

Nm  ~~  dV  »  -(Vp)dV  -  Ne[E  +  v  x  B]  dV  -  vNmvdV  (14) 

or 


juNmv  *  -Vp  -  Ne  [E  +  v  x  B]  -  vNmv 


(15) 


Now  Vp  may  be  written  according  to  the  chain  rule  of  differentiation  as 

vp  -  vp  (16) 

where  p  is  the  mass  density.  The  adiabatic  law  relates  p  and  p 

P  -  C,pY'  (17) 

where  y*  *  ratio  of  specific  heats.  As  discussed  In  Chapter  I,  y'  has  the 
value  of  3.0  for  an  electron  gas  provided  w>>v: 

Now 

fir  =  3cp2  us) 

dp 


3p/p 


(19) 


but  from  the  equation  of  state 


p  =  KTn 


(20) 
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where 

K  «  Boltzmann’s  constant 
T  *  electron  temperature 

n  *  variation  of  electron  density  from  its  mean  value, 
and  by  definition 

P  *  mn  (21) 


thus 


I2-  =  3KT/m 
dp 

At  this  point,  it  is  interesting  to  note  that  lE.  is  just  the  square  f  the 
propagation  velocity  y  that  describes  the  velocity  of  longitudinal  (sound) 
waves  in  compressible  fluids.  Thus 


(22) 


p  =  y  V(mn) 


(23) 


where 


y  =  3KT/m 


Thus  (15)  is,  in  linearized  parameters 


jwNmv  =  -y  V(mn)  -Ne  [E  +  v  x  B]  -  vNmv 


(24) 


or 


N[ jwv  +  a)uv  x  C  +  vv]  *  Jy  Vn  +  (Ne/m)  V$ 
H. 


(25) 


where 


OJ 


H 


eB/m,  electron  gyrofrequency 
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C  *  direction  of  magnetic  field  (unit  vector) 

4>  *  scalar  potential 

Up  to  this  point,  the  force  equation  is  valid  in  any  system  of  coordinates. 

A  further  attack  on  the  problem  will  now  require  that  the  cylindrical  coor- 

A  *  A 

dinates  r,  0  and  z  be  applied.  The  axis  of  the  antenna  probe  has  been  chosen 
to  coincide  with  the  z  axis  (see  Figure  8) . 


Z ,  ANTENNA  AXIS 


Fig.  8.  Coordinate  system. 


The  next  step  is  to  determine  the  velocity  v  in  terms  of  n  and  $ .  Because 
v  is  entangled  in  a  curl  operation,  this  is  not  a  straight  forward  process. 
We  begin  by  solving  for  the  velocity  components  separately.  In  the  r 
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direction. 


jwN  [Uv^.  -  jY(ve^z  ”  vz^e 


2  _3n  Ne  _3$ 

y  3r  m  3r 


(26) 


In  the  0  direction 


jt»N  [UvQ  -  jY(vzCr  -  vr?z)] 


2  i  Ne  jL 

11  r  30  m  r  30 


(27) 


In  the  z  direction 


j„N  [Uvz  -  jY(vrCe  -  v#tr)]  -  -P2  f  +  ?£  (28) 

where 

Y  *  (dh/u) 

U  -  1  -  jZ 

Z  *  v/uj 

However  because  of  symmetry 


3n 

3z 


0 


and 


3z 


0 


and  thus  from  (28) 


v 

z 


3K(Ve  -  Vr> 


(29) 


or 


v  =  jv  W  -  jv.W 
z  J  r  0  J  0  r 


(30) 
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and  (26)  and  (27)  become 


ju>UN  lvr  -  jWzv0  +  jW0(jvrW0  -  jvQWr)) 


-V 


2  3n  Ne  _3i 

Sr  m  Sr 


(31) 


where 


or 


ju>UN  [v0  +  jWzvr  +  jWr(jvrWQ  -  jvQWr)]  = 


2  1  3n  Ne  1  St 

“y  7  3?  +  77W  (32) 


w 

w 


Y/U 

WC 


wn  =  WC 


e 


w 


wc 


vrd  -  wp 


=  J8r  -  v0(wrw0  -  jwz> 


(33) 


ve(1-V  =  -  vr(wrwe  +  iw«> 


(34) 


where 


8  8 

n  Sr  p  Sr 


n 


1  dn  _  1  S£ 

n  r  30  p  r  S0 

2 

JL_ 

UNO) 


e 

Umw 


now  with  the  aid  of  the  relation 


2  2  2  2 

W  +  Wc  +  W  -  W 

r  0  z 


(35) 
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explicit  values  of  and  vg  can  be  obtained.  By  solving  the  equations 
simultaneously,  we  find  that 


Br  (1  -  Wr)  -  BQ(WrWe  -  jWg) 
1  -  w2 


(36) 


and 


6r(WrW0  4  jWj  -  Bn[l  -  wj] 


-j 


0' 


1  ”2 
1  -  rt 


(37) 


Now 


v  3v  -  3v_  3v 
,  __r  + r  1 0  + z 

r  3r  r  30  3z 


(38) 


Hence 


3v  3vfl 
r  o 


3r *  30  “““  3z 


3v 

and  — r—  must  be  calculated. 


3V 


r 

3r 


3V 


_0 

30 


iK  ..2.  369 


7^2 1  llr  (1  -  «;>  -  17  <“rw0  - 


2JL 


3 


[rf]  pS  ».“• 


(39) 


3W 


3W  3W 

0+*V  +  8r  (~3?  W0+  Wrle  +j  if* 


3B  ,  3W 

~30  ^  "  2S0^W0  "Hf* 


] 


(40) 


3v 

z 

3z 


0 


(41) 


But 


3W 
_ r 

30 


4-W, 


(42) 


i!e 

30 


-w 


(43) 


•••  ig'+V  ’  . 


1 


0 


Thus 


3 

3  r 


2  2 
o  JGa  .  8  _ili 

n  3r3e  p  3r30 


1  _i!n  _  1  _£±  _  1_  3n  JL_  li 

3r  n  r  3r30  p  r  3r30  n  ^2  30  p  ^2  30 


!!i  _  i  A  ii!« 

30  n  r  362  p  r  362 


and  thus 


V«v 


r  /  lo  to  i 

[ (  r  n  3r  ~  r 


3* 

3r 


(1-W2) 

r 


-<  h  6n  It  -  ^  8p  I  >  (WrW3  -  JV 
r  r  r 


+(8  ^-S  4)  (1-W2) 

n  3r2  p  3r2  r 


-(6n  7  ik  -  ep  7  77i7  >  <V.  -  JV 


2  ,2 

_(g  Jl  -9„J1  _  a  1.  ■  c  —  )  my  +  -jw  ) 
^  n  r  3r30  Bp  r  3r30  ;  ^  r  0  3  z' 


+(S  \  -  B  “  )  (W  H  -  JW  ) 

n  2  3b  p  2  36  r0  z 


30 


(44) 


(45) 


(46) 


(47) 


(48) 
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+(B ~ -  B  K  )  (l-W2) 
nr2  962  P  r‘  902  6 


(a  1  3n  .lit,  ,..2  ,2  v 

-IB  —  t—  -  B  —  —  }  U'  -  W  ) 
n  r  3r  p  r  5r  0  r 


-^nT-ft-ep^!!’  <V,» 

r  r  r 


(49) 


This  can  be  reduced  to 


V. 


v  =  — M  [ (0  V2n  -  B  ?2$)  -  W2Sin2y(B  S(n)  -  0  S(*))]  (50) 

[l-W2J  n  ?  n  p 


where 


2  ^2  2 
S(x)  =  Cos20  -  2  Sin0  Cos0  -  +  Sin20  L- 

ar2  r  3r30  r2  902 


+  Sin2 9  -  +  2  Sin0  Cos0  ~ 

r  dr  z  do 

r 


(51) 


This  result  can  now  be  substituted  into  the  continuity  equation 


N(V*v)  +  join  =  0 


(52) 


and  Poisson's  equation 


V  t  =  n e/e 


(53) 


can  also  be  applied  to  yield  the  result 


N 


|jv]  ,2n ' 


77—  (ne/e  ) 
Umu)  o 


W2Sin2y(B  S(n)  -  B  S (<J>) ) 
n  p 


1 


] 


+  jwn  *  0  (54) 


rpr-  *-”**•'■ - 
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Simplifying 

V2n  =  -a2n  +  W2Sin2vlS(n)-U(uj2  /y2)(e  /e)  S(4>)  ] 

N  0 

where 

2  2  2  2 

a  =  («T  -  1)  uj/u 

2  2  ? 

4/  =  (U  -  Y")/UX 

2  2 

Ujj  =  Ne  /meo 

2  2 
X  *  ul./w 

N 

W  =  Y/U 


(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


For  either  y<<l  or  W<<1,  this  equation  can  be  solved  directly.  However, 
the  general  case  requires  the  solution  for  a  pair  of  coupled,  second  order 
partial  differential  equations  whose  parameters  are  functions  of  the  space 
coordinates.  Clearly  a  straight  forward  analytic  solution  seems  to  be  out  of 
the  realm  of  possibility,  and  it  appears  that  an  approximation  technique  will 
be  required  in  order  to  solve  the  generalized  system  equations.  In  any  case, 
the  solution  for  y  =  0  is  still  of  the  highest  interest  since  this  only  mildly 
restricts  the  applicability  of  the  theory.  If  the  restriction  y  =  0  had  been 
applied  from  the  beginning,  the  derivation  would  have  been  a  little  less  com¬ 
plicated;  however  by  developing  (55),  we  now  have  an  equation  that  can  pos¬ 
sibly  be  solved  at  a  later  time  to  provide  an  estimate  of  antenna  impedance 
for  y  f  0.  Thus  for  the  purposes  of  this  report,  it  will  be  sufficient  to 
limit  the  solution  and  only  solve  the  restricted  equation  (61). 


V2n  = 


9 

-a“n 


(61) 
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Thus  the  plasma  may  be  described  by  the  system  equations  (53)  and  (61). 
For  y  =  0  however,  symmetry  exists  in  the  6  direction  so  that  only  the  r 
coordinate  need  be  considered.  Thus  it  is  sufficient  to  solve  the  set  of 
equations 

V2$  =  ne/e  (62) 

r  o  v 

2  2 

V^n  =  -an  (63) 

where  the  r  subscript  indicates  that  only  the  r  component  is  under  considera¬ 
tion. 

Solution  of  System  Equations 

Equation  (63)  may  be  expanded  and  transformed  by  performing  the  indicated 
operations  and  allowing  x  =  ar. 


dn 


d2n  1  _ 

,2  x  dx 
dx 


+  n 


(64) 


This  is  the  familiar  Bessel’s  equation  of  zero  order.  Hence,  it  has.  as  solu¬ 
tions,  all  of  the  various  Bessel  functions  of  order  zero.  Since  x  is  in  gen¬ 
eral  a  complex  variable,  it  is  convenient  to  use  a  linear  combination  of  the 
Hankel  functions  H  (x)  and  (x)  [sometimes  H  (x)  is  written  as  (x)  ] . 

OO  0  O' 

In  general,  there  are  two  values  of  x  to  be  considered,  x  =  -v^2"  r  and 

x  =  +*a  r.  However,  it  is  a  property  of  the  Hankel  functions  [Jahnke  and 

(2)  (2) 

Emde,  1945]  that  H  (-x)  =  H  (x)  and  that  H  (-x)  =  Hv  (x) .  Thus  we  need 

0  0  0  0 

to  consider  only  the  one  solution  of  x  which  lies  in  the  first  two  quadrants 
of  the  complex  plane.  Then  for  large  values  of  x,  that  is,  at  distances 

far  removed  from  the  antenna,  we  require  that  the  solution  for  n,  the  den- 

(  2  ) 

sity  variations,  become  small.  Thus  since  Hq  (x)  approaches  infinity  as 
x  approaches  zero,  only  H(x)  can  be  an  acceptable  solution. 


Therefore, 


n  =  noHQ(ar)  (65) 

where 

*  arbitrary  constant  to  be  determined  by  boundary  conditions 

2 

a  =  the  square  root  of  a  that  lies  in  th ^  top  half  of  the 
complex  plane. 

The  solution  to  the  homogeneous  part  of  the  system  equation  (62)  may  be 
written  as 

*  -  Cx  ln(r/L)  (66) 

where  is  an  arbitrary  constant  yet  to  be  determined  and  L  is  the  length 
of  the  antenna.  This  form  of  solution  was  chosen  to  produce  the  correct 
form  of  potential  near  the  surface  of  a  long,  thin,  charged  rod  in  free 
space.  Since  (62)  is  very  similar  to  (63),  a  solution  of  the  form 

*  ”  C2Ho^ar^  (67} 

would  also  be  expected  to  satisfy  (62).  Such  is  indeed  the  case  and  the 
total  solution  may  be  written  as 


4>  *  C.ln(r/L)  -  (en  /e  o^)  H  (ar) 
i  o  o  o 

2 

where  C2  =  -n^e/e^a  has  been  chosen  so  as  to  satisfy  (62). 


(68) 


Boundary  Conditions 


General  solutions  for  n  and  $  now  exist;  however,  the  arbitrary  con¬ 


stants  C,  and  n  have  not  yet  been  determined.  These  constants  can  be 
i  o 

obtained  by  a  consideration  of  boundary  conditions  near  the  antenna  surface. 
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Our  original  pair  of  second  order  differential  equations  require  the 
specification  of  four  independent  boundary  conditions  to  determine  n  and  $. 
Two  of  these  conditions  were  fixed  in  the  derivation  above  when  both  n  and 
4>  were  required  to  approach  zero  as  the  distance  from  the  antenna  increased 
indefinitely.  Only  one  of  the  other  two  conditions  need  be  specified  if 
only  the  antenna  impedance  is  desired.  This  is  true  since  we  desire  not  a 
knowledge  of  <p  or  n  but  only  a  knowledge  of  the  ratio  of  antenna  voltage  to 
current.  Ultimately,  one  of  the  arbitrary  constants  n^  or  must  cancel 
out  of  the  impedance  expression.  Thus  if  a  boundary  condition  at  the 
antenna  surface  can  be  specified  such  that  CL  is  determined,  it  will  be 
possible  to  later  cancel  out  n^  and  completely  determine  the  antenna  im¬ 
pedance  . 

It  has  previously  been  assumed  that  the  plasma  surrounding  the  antenna 
is  homogeneous;  that  is,  that  no  sheath  exists.  Immediately  inside  the 
antenna-plasma  boundary,  the  electrons  are  no  longer  in  a  free,  gaseous  state 
and  are  thus  essentially  restrained  from  any  movement.  Hence,  it  is  assumed 
that  the  velocity  v  is  zero  at  the  antenna  surface. 

This  is  the  same  boundary  condition  used  by  Fejev  [1964]  in  a  similar 
problem.  As  he  points  out,  it  does  not  correspond  exactly  to  the  physical 
situation,  but  only  approximates  it.  More  exact  boundary  conditions  can 
become  very  difficult,  but  they  can  also  lead  to  a  more  general  solution. 
Balmain  [1965],  for  example,  has  considered  the  complex  boundary  conditions 
that  actually  exist  at  a  probe-plasma  boundary,  and  it  is  possible  that  his 
technique  could  be  applied  in  this  case.  However,  it  is  not  necessary  to 
include  all  or  these  complex  conditions  if  only  high  frequency,  small  amp¬ 
litude  effects  are  of  interest  as  is  the  case  in  this  derivation. 
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For  v  *  0  at  the  antenna  surface  (r  *  R),  equation  (25)  becomes  at  the 
surface 

(69) 


In2?11!  D 

r  r"R 


l(Ne/m)  Vr*]r=R 


2 

iL 

e 


iKL  ■•:[«].. 


(70) 


but  from  (68)  and  the  fact  that  —  [H  fx)]  =  -H.(x) 

dx  o  i 


ii 

3r 


C./r  +  (en/c  a)  H  (ar) 
1  o  o  1 


(71) 


and  similarly 


-r—  *  -n  aH  (ar) 
dr  o  1 


(72) 


Thus 


-y  (e/e  )(n  aH.(aP)) 
o  o  i 


w  (C./R  +  (en  /e  a)H  (aR)) 

NX  O  0  1 


(73) 


or 


C,  *  -(en  R/e  a)lL  (aR)  ip 
1  o  o  i 


(74) 


since 


,2  ,22.2 

4>  -  1+u  a  /w... 

N 


(75) 


Thus  the  expression  for  the  potential,  (68),  becomes 


2  2 

J>  -  -(n^e/e^a  )  (ar)  +  aRiJ;  (aR)  In  (r/L) 


- .  «  nw  11  1 1 '  ifpwi  m  ^JSMTT'  1  |  ~  » 


(76) 


Antenna  Impedance 


Impedance  is  defined  as  the  complex  ratio  of  voltage  to  current.  In  the 
case  of  the  antenna  under  consideration  here,  the  voltage  is  just  the  poten¬ 
tial  at  the  surface  of  the  antenna  and  the  current  is  just  the  displacement 
current  summed  over  the  antenna  surface.  Thus 


Z 

a 


$(R) 


or  since 


-¥ 


3D 

3t 


jwD 


(77) 


(78) 


*(R) 

2*RLju>eoEr(R) 


(79) 


but 


E 

r 


(80) 


and 


*  (n  e/e  u2)(aH  (ctr)  -  aRiJ^H  (aR)/r) 
or  o  o  1  1 

Thus  (81)  becomes 


(81) 


Z 

a 


,  K  (aR)-uR^2H. (oR)ln(L/R) 
1  o _ I _ 

jw2TTRLeQ  aHi(aR)-aRi|;2H1(aR)/R 


(82) 


or 


ln(L/R)  HQ(aR)/aRH1(aR)  ln(L/R)  -  f 


j w2iTe  L 


1  -  ip 


(83) 
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Now  the  free  space  capacitance  of  a  long  thin  rod  is  easily  shown  to  be 


C  -  2 *e  L/ln(L/R) 

O  0 

Thus 

-1  fd(aR)  /ln(L/R)-f2 
a  ujC  .  ,2 

o  L  i  - 


(84) 


(85) 


whe  re 


*(aR)  =  H  (aR) /aRH. (aR) 
o  1 


(86) 


This  is  the  impedance  of  a  monopole  antenna  in  a  magneto-ionic  plasma  with 
the  antenna  axis  aligned  with  the  magnetic  field. 

The  expression  includes  electron  collision  effects  as  a  parameter  in 
a.  If  v  *  0  and  X  M,  then  a  becomes  entirely  imaginary  and  4>  can  be  ex¬ 
pressed  in  terms  of  the  real  modified  Bessel  functions  of  the  second  kind 

K  (x)  and  K  (x) 
o  1 


Thus 


rf(*jx) 


HoOjx) 

'jxHjOjx) 


K0(x) 

xKxU) 


(87) 


We  will  find  this  result  useful  in  some  cases  where  v  =  0. 

2 

If  w  =  0,  then  ip  becomes 
H 

ip2  =  u/x 


(88) 


while  if  v  =  0  as  well,  then  a  particularly  simple  expression  results 


tp  =  oj/ojn 


(89) 
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and 


2  /  2  2w  2 

a  *  (w  -u)yJ/u 


(90) 


or 


»  r/  2  2.  ,  2, 

c*  “  l(uN  -  0)  )/u  J 


1/2 


(91) 


where 


,2  2 
a  =  -a 


(92) 


This  expression  is  useful  when  (87)  is  used  to  determine  as 


rf(a'R)  =  -r 


K  (a'R) 
o 


a'RK1(a'R) 


(93) 


Sheath  Effects 

The  first  order  effects  of  the  plasma  sheath  on  the  antenna  impedance 
can  be  taken  into  account  by  use  of  a  very  simple  physical  model.  As  long 
as  the  dc  potential  of  the  antenna  is  below  the  plasma  potential  (this  is 
usually  the  case),  a  region  of  electron  depletion  will  surround  the  antenna. 
The  mechanism  of  this  sheath  formation  is  discussed  in  Chapter  I.  The 
model  assumes  that  no  electrons  exist  within  a  region  of  S  Debye  lengths 
of  the  antenna.  Generally  S  will  vary  from  0  when  the  antenna  is  at  the 
plasma  potential,  through  about  4  when  no  net  dc  current  is  allowed  to 
flew  to  the  antenna,  and  up  to  about  10  if  large  negative  potentials  are 
applied  to  the  probe.  Thus  the  effective  radius  Rg  of  the  sheath  is 

R„  =  R.  +  SA 
S  AD 


(94) 


The  capacitance  between  the  antenna  and  the  surface  of  the  plasma  is  then 


the  capacitance  of  two  concentric  cylinders  of  radii  and  R<,. 

C_  -  2tte  L/ln(R  /R  )  (95) 

b  0  b  A 

The  effective  antenna  impedance  is  then  the  impedance  of  this  capac¬ 
itor  in  series  with  the  impedance  of  an  antenna  of  radius  Rg  since  the 
plasma  begins  at  a  radius  of  Rg.  The  general  impedance  relation  is  then 


±  ±  rf(aRs)/ln(L/Rs)-*2 

za  '  +  ! .  / 

where 

C  -  2  ire  L/ln(L/R_) 

o  o  b 


(96) 


(97) 


The  derivation  of  the  antenna  impedance  is  now  complete.  It  can  be 
seen  that  the  major  goals  have  been  met;  that  is  the  effects  of  a  hot, 
lossy,  magneto-ionic  plasma  on  antenna  impedance  have  been  calculated. 

It  should  be  noted  that  this  expression  is  valid  when  the  axis  of  the 
antenna  is  parallel  to  the  magnetic  field.  A  summary  o2  the  final  ex¬ 
pression  and  related  parameters  is  tabulated  in  Table  I. 


*** 


Summary  of  Antenna  Impedance  Theory 


Sheath  Scaling  Constant 


Chapter  IV 


COMPARISON  OF  RESULTS 


Introduction 

It  is  desirable  to  empirically  examine  the  impedance  expression 
developed  in  Chapter  III  in  order  to  compare  it  to  both  experimental 
data  and  other  theories.  Generally,  other  theoretical  expressions  are  of 
completely  different  form  and  often  either  gross  approximations  or  em¬ 
pirical  methods  must  be  applied  before  a  favorable  comparison  can  be 
discerned.  Thus,  the  theoretical  expression  developed  here  and  the  other 
theoretical  results  of  Chapter  II  have  been  programmed  for  digital  computer 
analysis  so  that  empircal  comparisons  can  be  carried  out. 

If  the  conventional  free  space  expression  for  antenna  impedance  of 
a  short  dipole  is  examined  [Balmain,  1964,  p.  55],  it  is  found  that 
the  impedance  in  the  low  frequency  limit  is  just  the  impedance  of  the  free 
space  capacitance  given  by 

Cq  =  2TT£oL/(ln(L/R)-l)  (98) 

This  expression  is  the  same  as  the  Cq  derived  in  Chapter  III  except  for  the 
-1  term.  This  term  represents  the  end  effects  that' were' ignored  in  the 
formulation  of  Chapter  III.  The  above  expression  could  be  used  to  define  C 

o 

but  it  is  even  more  preferable  to  use  an  experimentally  measured  value  of 
free  space  capacitance  so  that  all  perturbating  influences  can  be  taken  into 
account  (rocket  body,  etc.).  Thus  the  actual  mean  antenna  radius  will  be 
used  in  the  impedance  expressions,  but  the  length  of  the  antenna  in  the 
calculations  will  be  changed  from  its  true  value  toa  value  that  will 
produce  the  correct  value  of  Cq  as  determined  by  measurements. 


■  l-t 
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It  could  be  argued  that  end  effects  could  modify  more  than  just  the 
free  space  capacitance.  These  effects  should  be  small,  however,  since  the 
plasma  medium  tends  to  shield  the  antenna  from  any  external  disturbing 
influences  at  distances  greater  than  a  few  Debye  lengths  (see  Chapter  I) . 

This  is  especially  true  at  frequencies  below  the  plasma  frequency  and  at 
fremiencies  above  the  plasma  frequency,  the  end  effects  are  not  especially 
critical. 

The  empirical  comparison  of  impedance  theories  can  be  approached  from 
many  directions.  However,  it  is  desirable  to  examine  theoretical  results 
so  that  they  can  easily  be  compared  to  experimental  measurements  as  well. 

Thus  the  chosen  antenna  dimensions  are  typical  of  experimental  antennas. 

In  analogy  with  the  Plasma  Frequency  Probe,  the  plasma  parameters  can  all 
be  held  constant  while  the  exciting  frequency  is  swept  through  a  wide 
range  of  values  (.1  to  10  MHz)  and  the  impedance  determined.  The  theory  can 
also  be  explored  by  fixing  the  exciting  frequency  and  determining  the  antenna 
impedance  as  a  function  of  the  dominant  plasma  variable , plasma  frequency 
with  electron  temperature  as  a  parameter.  This  situation  will  then  corre¬ 
spond  to  the  Standing  Wave  Impedance  Probe.  It  is  also  possible  to  plot 
antenna  resistance  versus  antenna  reactance  with  both  plasma  frequency  and 
temperature  as  parameters.  Perhaps  only  the  dominate  characteristics  of  the 
impedance  are  desired,  such  as  the  series  and  parallel  resonance  frequencies 
of  the  antenna.  This  situation  actually  more  closely  approximates  the 
plasma  frequency  probe  which  determines  only  the  frequencies  of  zero  phase 
or  resonance.  Thus,  these  four  methods  will  all  be  used  to  illustrate 
antenna  impedance  throughout  the  remainder  of  this  paper. 

A  computer  program,  (SWP4MD)  was  written  to  plot  the  antenna  impedance 
as  a  function  of  either  plasma  frequency  or  exciting  frequency.  Various 
subroutines  were  written  to  calculate  the  antenna  impedance  according  to  the 
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different  theories  of  interest.  Thus,  empirical  impedance  plots  of  any  of 
these  theories  could  be  obtained  by  use  of  the  proper  subroutine  in  the 
main  computer  program.  Fortran  listings  of  these  programs  are  included  in 
Appendix  B  for  reference. 

Cold  Plasma  Comparison 

It  was  mentioned  earlier  that  most  of  the  proposed  theories  of  antenna 
impedance  are  for  the  cold  plasma  case.  Thus  the  first  test  is  to  compare 
the  derived  impedance  relations  with  other  theories  for  zero  electron 
temperature.  For  zero  temperature,  however,  the  theory  of  Chapter  III  re¬ 
verts  to  the  cold  plasma  magneto-ionic  theory  of  Appendix  A.  This  theory 
agrees  very  well  with  other  results  and  is  discussed  and  empirically  compared 
to  Balmain's  Theory  in  Appendix  A. 

For  non-zero  electron  temperature,  the  author  is  only  aware  of  two 
other  theories  for  cylindrical  antenna  impedance;  that  of  Balmain  [1964] 
and  that  of  Whale  [1965] .  Both  considered  the  case  of  finite  electron 
temperature,  but  only  for  no  magnetic  field. 

Balmain  Finite  Temperature  Theory 

Balmain's  theoretical  result  contained  several  difficult  transendental 
functions , and  he  found  it  expedient  to  assume  that  the  parameter  a'R  (see 
below)  is  large  compared  to  one.  This  approximation  greatly  simplified  his 
more  exact  result.  The  resulting  expression  is 

ln(L/R) -1  +  (K  -l)/2a'R  ,0Q* 

Z  =  - 0 -  (99) 

3  j  oj2tt e  K  L 

0  o 

where 


K  =  1-X/U 

o 

=  (X-l)  1/2o)/p 


a' 


Ub 


Whale  Theory 

Whale  (1963j,  considered  antenna  power  losses  due  to  electroacoustic 
waves.  Therefore,  his  results  apply  to  an  antenna  in  a  warm  plasma.  His 
expression  can  be  used  to  estimate  antenna  impedance  if  simple  ionic  theory 
is  assumed  for  the  antenna  reactance  while  the  power  loss  is  interpreted 
as  the  result  of  antenna  resistance.  Whale's  power  loss  equation  is 


Xj82V2(C  /L) 

n  =  - ^—2 - 

^1  4ttc  R  ^1-X 
o  s 


(100) 


where 

-  Radiated  power  per  unit  length 

8  =  A  special  function  calculated  by  Whale  and 

presented  in  the  form  of  a  graph 

=  Antenna  voltage 


If  this  expression  is  used  to  estimate  antenna  resistance  and  magneto¬ 
ionic  theory  is  used  to  find  the  antenna  reactance,  the  result  is 

2  2  1/2 

za  -  v2  •  iy2)  -  v  -  (10i) 

o 

where 


R,  =  2-ne  LR  vl-X  /XB 
L  os 


2 


(102) 


Comparison  of  Theories 

Both  of  the  above  expressions,  as  well  as  the  expressions  for  antenna 
impedance  from  the  Despain  and  magneto  - ionic  theories  were  programmed  as  com¬ 
puter  subroutines  with  the  name  ANTZ.  Thus  all  the  theories  could  be 
simulated  and  later  compared  by  using  the  desired  subroutine  ^n  the  main 


computer  program. 


Balmain's  expression  is  plotted  as  a  function  of  exciting  frequency  in 
Figure  9.  Figure  10  is  a  similar  plot  of  the  Despain  theory  with  the  same 
parameters.  For  1*0,  the  results  from  each  theory  are  identical.  However, 
as  the  temperature  is  increased,  some  discrepancies  appear.  These  are 
especially  evident  in  the  antenna  resistance  when  the  exciting  frequency  is 
above  the  plasma  frequency.  Near  T  =  1000 °K,  there  is  reasonably  good 
agreement,  but  as  the  temperature  is  increased  further,  it  can  be  seen  that 
the  shapes  of  corresponding  curves  begin  to  show  radical  differences.  This 
is  very  evident  just  below  the  plasma  frequency.  This  discrepancy  should 
not  be  unexpected  since  for  high  temperatures  and  for  exciting  frequencies 
near  the  plasma  frequency,  a'R  becomes  small  and  the  large  a'R  approximation 
made  by  Balmain  is  not  valid.  Balmain's  complete,  though  very  difficult 
expression,  must  be  used  if  his  theory  is  to  apply  under  these  conditions. 

It  should  be  noted  that  in  Balmain's  theory  (in  contrast  to  the 
Despain  theory}  no  appreciable  reactance  variation  occurs  as  a  function  of 
temperature  for  exciting  frequencies  above  the  plasma  frequency.  This  lack 
of  variation  is  viewed  with  suspicion  since  even  in  regions  where  a'R  is 
large,  it  seems  only  reasonable  that  antenna  reactance  is  always  changed 
to  some  degree  by  a  change  in  electron  temperature  and  antenna  resistance. 
Appendix  A  also  includes  a  discussion  of  a  problem  that  occurs  in  Balmain's 
theory  under  certain  special  conditions. 

It  can  thus  be  seen  that  the  general  characteristics  of  the  two  theories 
agree  in  this  particular  example  for  temperatures  below  1000°K. 

Figure  11  contains  a  comparison  of  the  Balmain  and  Despain  theories  for 
fixed  exciting  frequency  with  plasma  frequency  as  a  variable.  Again,  godd 
agreement  is  noted  for  low  temperatures  and  for  operating  frequencies  well 
separated  from  the  plasma  frequency.  Again,  somt  disagreement  is  evident 
wherever  a'R  becomes  small.  In  fact,  in  the  Balmain  theory,  the  resistance 
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Fig.  9.  Balmain  Theory  Antenna  Impedance. 
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is  negative  for  high  electron  temperatures  in  certain  regions.  This 
is  clearly  an  impossible  result  for  the  physical  system;  and  as  discussed 
above,  comparisons  in  these  regions  are  not  possible.  Again,  it  appears 
that  the  Despain  and  Balmain  theories  generally  agree  over  those  regions 
where  valid  comparisons  can  be  made. 

Although  only  a  few  sets  of  parameters  have  been  examined  here*  more 
extensive  comparisons  of  the  two  theories  have  been  made  for  several  other 
sets  of  parameters  typical  of  ionospheric  experiments.  In  each  case,  the 
conclusions  are  identical  to  those  discussed  above. 

In  addition  to  the  results  from  the  theories  of  Balmain  and  Despain, 
the  impedance  expression  as  derived  from  Whale's  theory  is  also  presented 
in  Figure  11.  It  can  be  noted  Lhat  all  the  theories  predict  finite  antenna 
resistance  for  finite  electron  temperature  over  at  least  part  of  the  fre¬ 
quency  range.  Also  the  resistance  increases  with  increasing  plasma  fre¬ 
quency  in  each  case.  The  similarity  ends  at  this  point,  however,  since  the 
Whale  theory  resistance  curve  is  clearly  in  variance  with  both  the  Balmain 
and  Despain  theory.  It  is  felt  that  because  both  the  Despain  and  Balmain 
theories  were  more  rigorously  derived  than  Whale's  work,  that  they  are 
the  more  accurate  formulations  of  antenna  impedance.  Therefore,  Whale's 
work  will  not  be  considered  further. 

The  Despain  theory  has  been  compared  to  other  theoretical  results. 

It  was  found  to  agree  with  Balmain's  work  where  comparisons  were  possible. 
It  did  not  agree  with  the  Whale  theory.  Thus  the  Despain  theory  seems  to 
be  at  least  as  good  theoretically  as  the  Balmain  theory  and  yet  it  was 
derived  from  a  somewhat  different  viewpoint  and  is  somewhat  more  general 


in  its  formulation. 
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Important  Implications  of  the  Theory 

There  are  several  very  important  implications  that  can  be  noted  in 
the  examples  of  the  Despain  theory  giv.-n  in  Figures  10  and  11.  First, 
for  no  electron  collisions  and  zero  electron  temperature,  no  antenna 
resistance  appears.  However,  if  finite  electron  temperatures  exist,  then 
antenna  resistance  appears  for  exciting  frequencies  above  the  plasma 
frequency.  Thus  one  effect  of  electron  temperature  is  to  add  a  loss 
mechanism  not  considered  in  the  usual  cold  plasma  formulations.  Second, 
it  can  be  observed  that  for  exciting  frequencies  below  the  plasma  frequency, 
a  series  (low  impedance)  resonance  exists,  whose  frequency  is  a  sensitive 
function  of  the  electron  temperature.  Above  the  plasma  frequency,  the  elec¬ 
tron  temperature  has  a  much  smaller  effect  on  the  antenna  impedance.  Third, 
the  electron  density  (plasma  frequency),  temperature,  and  collision  fre¬ 
quency  all  greatly  affect  antenna  impedance.  Hence,  it  should  be  possible 
to  directly  determine  those  parameters  from  antenna  impedance  measurements. 
This  problem  is  the  subject  of  the  next  chapter. 


Further  Examination  of  the  Theory 

The  effect  of  magnetic  field  has  not  yet  been  investigated.  There¬ 
fore,  in  Figure  12,  the  gyrof requency  has  been  taken  as  1.4  MHz.  By  com¬ 
paring  Figures  10  and  12,  it  can  be  seen  that  the  primary  effect  of  the 
magnetic  field  is  to  transform  the  frequency  scale  so  that  the  zero  reactance 
point  of  the  zero  temperature  curve  is  transformed  up  to  the  gyrofrequency 
and  the  parallel  resonance  (high  impedance)  point  is  shifted  from  the  plasma 
frequency  to  the  upper  hybrid  frequency  given  by 
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Fig.  12.  Antenna  impedance  in  the  presence 
a  longitudinal  magnetic  field. 


If  n  close  examination  is  made,  it  is  found  that  the  impedance  at  f 
with  no  magnetic  field  is  about  the  same  at  f'  with  longitudinal  mag¬ 
netic  field  where 
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The  additional  effect  of  collisions  is  shown  in  Figure  13.  It  can 
be  noted  that  the  resonances  are  not  changed  in  frequency  although  some 
change  is  evident  in  the  shape  of  the  reactance  curves.  Furthermore, 
large  values  of  resistance  are  evident  throughout  the  complete  frequency 
range;.  Thus  the  primary  result  of  electron  collisions  is  an  increase  in 
antenna  resistance. 


Comparison  of  the  Theory  to  Experimental  Results 

There  are  two  basic  methods  for  comparing  the  theory  to  experimental 
results.  One  method  is  to  use  the  theory  to  simulate  antenna  impedance 
from  the  best  independent  estimate  of  the  experimental  plasma  parameters 
and  to  then  compare  these  results  with  the  experimental  antenna  impedance 
values.  The  other  method  is  to  use  the  theory  to  determine  the  plasma 
parameters  from  the  experimentally  determined  antenna  impedance  values. 

The  resulting  plasma  parameters  are  then  compared  to  the  independent 
estimate  of  the  actual  plasma  parameters.  The  first  technique  will  be  em¬ 
ployed  in  this  section  while  the  second  will  be  used  in  the  next  chapter. 

There  are  three  crucial  impedance  characteristics  predicted  by  the 
theory.  The  first  is  a  parallel  resonance  condition  for  an  exciting  fre¬ 
quency  near  the  upper  hybrid  frequency  described  by 
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Fig.  13.  Antenna  impedance  with  longitudinal 
magnetic  field  and  electron  collisions. 
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The  second  is  the  existance  of  an  electron  tern  »rature  controlled  series 
resonance  condition  for  an  exciting  frequency  between  the  gy rofrequency  and 
the  upper  hybrid  frequency.  The  third  is  the  appearance  of  appreciable 
antenna  losses  for  exciting  frequencies  above  the  upper  hybrid  frequency 
even  when  the  electron  collision  losses  are  neglegible.  Thus  each  or  these 
characteristics  should  be  observed  in  experimental  results  whenever  the 
proper  conditions  are  met. 

Experimentally,  series  and  parallel  resonances  have  been  observed 
many  times  in  results  from  the  Plasma  Frequency  Probe.  In  fact,  whenever 
it  could  be  established  that  according  to  the  theory,  that  these  resonances 
should  occur,  they  were  indeed  observed.  In  most  cases,  independent  measure¬ 
ments  of  plasma  frequency  were  also  made  so  that  the  parallel  resonance 
condition  could  be  checked.  It  was  indeed  at  the  predicted  hybrid  frequency 
in  every  case.  Several  such  comparisons  are  discussed  in  a  paper  by 
Baker 3  Despain  and  Vlwick  [1966]  Independent  measurements  of  electron 
temperature  were  not  made  in  most  cases.  Thus  even  though  the  series  resonances 
were  consistently  observed  in  all  cases,  it  is  not  possible  to  determine  at 
exactly  what  frequency  they  should  have  occurred.  A  typical  example  from 
one  of  these  experiments  is  shown  in  Figure  14.  Each  experimental  measure¬ 
ment  is  shown  as  a  dot.  A  curve  has  been  drawn  through  the  mean  of  these 
points  as  a  representation  of  the  parallel  frequency.  There  was  evidence 
for  this  flight,  that  vehicle  potential  was  rapidly  changing  over  a  several 
volt  range.  Thus  it  has  not  been  possible  to  completely  determine  the 
antenna  potential  except  to  note  that  it  varied  repeatedly  from  about  zero 
down  to  about  minus  five  volts  during  a  period  of  about  1.6  seconds.  Sheath 
effects  that  correspond  to  a  value  of  S  from  about  zero  to  four  must  therefore 
be  effective  here.  This  accounts  for  some  of  the  spread  in  the  series  resonance 
data.  So  that  this  data  could  be  compared  to  theory,  a  pair  of  curves  were 


Fig.  14.  Experimental  series  and  parallel  resonance  frequencies 


57 


drawn  near  the  series  resonance  points.  Ihese  curves  represent  the  pre¬ 
dictions  of  the  theory  using  the  parallel  resonance  curve  and  the  indicated 
sheath  and  electron  temperature  parameters.  These  parameters  are  re¬ 
spectively,  the  highest  and  lowest  limits  that  can  be  reasonably  estimated 
for  this  flight  tor  the  sheath  variable  S  and  electron  temperature  T.  It 
can  be  seen  that  in  this  particular  case,  and  within  the  experimental  un¬ 
certainties,  there  is  good  agreement  between  theory  and  experiment.  This 
same  conclusion  was  reached  after  an  examination  or  the  other  flights  where 
good  measurements  were  obtained. 

The  observation  of  high  antenna  losses  in  regions  of  low  electron 
collision  frequency  have  been  observed  by  Whale[ 1964]  as  mentioned  earlier. 

A  recent  experiment  using  the  Standing  Wave  Impedance  Probe  also  exhibited 
a  similar  phenomena  [Baker  et  al  ,  1965).  Since  impedance  measurement^ 
of  electron  temperature  are  not  available,  it  is  again  impossible  to  pre¬ 
dict  exactly  what  losses  should  have  been  observed. 

;  t  is  possible  to  assert  that  at  least  these  experiments  demonstrate 
results  consistent  with  the  theory  and  no  contradictions  ate  yet  evident 
Such  is  not  the  case  for  the  various  other  theories  discussed  earlier. 

The  cold  plasma  theories  and  Whale's  results  do  not  predict  a  series 
resonance, and  Balmain's  warm  plasma  theory  ignores  magnetic  fields  that 
are  always  present  in  the  ionosphere.  Thus  Balmain's  theory  does  not  pre¬ 
dict  correct  values  of  series  resonance  or  antenna  losses.  It  is  not  pos¬ 
sible  to  assert  that  the  Despain  theory  is  absolutely  correct,  but  it  has 
been  proved  that  it  is  more  accurate  for  the  range  of  parameters  of  interest 
here  than  have  any  of  the  other  theories  of  which  the  author  is  aware. 


Chapter  V 


DETERMINATION  OF  PLASMA  PARAMETERS 


Introduction 

Historically,  RF  antenna  measurements  have  been  used  to  determine 
electron  density  (plasma  frequency)  and  electron  collision  frequency,  using 
one  of  the  magneto-ionic  theories  in  the  interpretation  of  the  data.  From 
the  theory  developed  in  this  paper,  it  can  be  seen  that  electron  temperature 
must  also  be  considered;  indeed,  it  is  possible  to  determine  electron  tem¬ 
perature  from  the  antenna  impedance  measurements  using  this  theory.  The 
purpose  of  this  chapter  is,  therefore,  to  demonstrate  how  electron  density, 
temperature  and  collision  frequency  can  all  be  determined  by  use  of  this 
theory. 

Antenna  impedance  is  described  by  two  parameters;  resistance  and  re¬ 
actance.  It  is  desired,  however,  to  determine  three  parameters;  that  is, 
electron  density,  temperature  and  collision  frequency.  Thus  either  impedance 
measurements  must  be  made  at  (a)  more  than  one  exciting  frequency  or  (b) 
simultaneously  on  at  least  tw>  antennas  of  significantly  different  con¬ 
figuration  or  (c)  there  must  be  some  known  relation  between  at  least  two 
of  the  plasma  parameters. 

Over  the  usual  range  of  natural  ionosphere  conditions,  there  is  a  simple 
relationship  between  electron  collision  frequency,  electron  temperature  and 
total  particle  density  ( Pfister ,  1965].  The  relation  is 

v  =  CNTNt  (106) 

3 

where  N  =  Total  particle  density  (particles/meter  ) 

C„  =  A  constant  derived  from  any  desired  model  atmosphere, 

W  -17 

taken  tc  be  2.6  x  10  in  this  paper. 
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The  total  density  can  be  found  from  a  model  atmosphere  once  the  desired 
altitude  is  specified.  The  altitude,  in  turn,  is  known  from  the  rocket 
trajectory.  Thus  the  parameters  altitude  and  antenna  impedance  can 
specify  electron  density,  temperature,  and  collision  frequency. 

Plasma  Parameters  from  the  Standing  Wave  Impedance  Probe 

The  Standing  Wave  Impedance  Probe  determines  antenna  impedance  at  a 
single  frequency  and  at  a  given  time.  Thus  for  a  given  altitude,  a  plot 
of  antenna  impedance  as  a  function  of  plasma  frequency  and  electron  tempera¬ 
ture  can  be  used  to  determine  these  parameters .  Such  a  plot  is  generated 
by  a  computer  program  known  as  IMPPLT  in  conjunction  with  the  subroutine 
ANTZ.  These  programs  are  listed  in  Appendix  B  for  reference. 

A  plot  of  antenna  resistance  versus  reactance  is  shown  in  Figure  15 
for  v  =  0,  and  for  no  magnetic  field.  It  describes  a  typical  antenna  used 
in  ionospheric  experiments  of  3  meter  length,  1/2  cm  radius  and  driven  at 
3  MHz.  The  impedance  values  shown  are  the  differential  impedance  changes 
from  the  free  space  conditions.  Curves  are  shown  for  various  constant  electron 
temperatures  and  also  for  various  constant  plasma  frequencies.  Figure  16  is 
a  similar  plot  for  v  =  1  MHz  in  which  the  effect  of  electron  collisions  can 
be  discerned.  The  primary  effect  of  the  collisions  is  to  increase  the  antenna 
resistance.  In  Figure  17,  the  collisions  are  assumed  to  follow  the  relation 
given  in  equation  (106)  at  75  km.  In  addition,  magnetic  field  effects  have 
also  been  taken  into  account  by  assuming  a  gyrof requency  of  1.4  MHz.  It 
is  evident  that  given  antenna  impedance  and  altitude,  the  plasma  frequency 
and  electron  temperature  can  be  determined  directly,  and  that  equation  (106) 
can  specify  electron  collision  frequency.  Thus  these  parameters  can  all 
be  determined  from  the  results  of  the  Standing  Wave  Impedance  Probe. 
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ANTENNA  REACTANCE  (OHMS) 


Fig.  17.  Antenna  Resistance  versus 
reactance  for  75  km. 
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Plasma  Parameters  from  the  Plasma  Frequency  Probe 

The  Plasma  Frequency  Probe,  as  explained  in  Chapter  I,  directly 
determines  the  series  and  parallel  resonance  frequencies  of  an  antenna 
that  occur  in  the  range  of  .1  MHz  to  10  MHz.  It  was  shown  in  the  last 
chapter  that  the  resonance  frequencies  are  not  modified  by  modest  electron 
collision  frequencies  so  that  v  is  assumed  to  be  zero  in  the  work  that 
follows.  The  parallel  frequency  is  not  a  function  of  electron  temperature, 
and  hence  it  defines  the  plasma  frequency.  The  relationship  for  longi¬ 
tudinal  magnetic  field  is 


For  other  magnetic  aspect  angles,  it  can  be  easily  shown  from  magneto¬ 
ionic  theory  (Appendix  A)  that  this  equation  produces  at  most  an  error 
of  about  10%  in  plasma  frequency.  For  f  »f^,  the  error  is  very  small. 
The  series  frequency,  on  the  other  hand,  is  a  function  of  plasma  freq¬ 
uency,  gyrof requency  and  electron  temperature.  Thus  it  is  desirable  to 
calculate  the  aeries  resonance  frequency  as  a  function  of  the  antenna  and 
plasma  parameters  so  that  this  relationship  can  be  examined. 


Calculation  of  Series  Resonance  Frequency 

The  condition  for  series  resonance  is  that  the  imaginary  part  of 
the  antenna  impedance  be  equal  to  zero.  Thus  from  Table  I 


or 


Imag 


1  (  1 

/  4>(aR  )/ln(L/R  )-iA 

juiC^  juiCq 

(  i-,’-  11 
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If  ip  is  found,  the  series  frequency  will  be  defined.  Thus  since  v  “  0, 

0  and  are  both  real  for 

fH  <f  <  fp  (110) 

then 

<//2  -  [ln(R  /R. )  +  <KaR  )  ] /ln(L/R  )  (111) 

S  A  S  A 

where 


C  /C  *  ln(R  /R  ) /ln(L/R  ) 
os  s  A  s 


(112) 


Equation  (111)  is  the  series  respnance  condition.  It  is  a  transcendental 

equation  and  is  therefore  not  subject  to  analytical  solution.  It  can  be 

solved  by  a  numerical  technique  such  as  Newton's  iterative  procedure. 

Subroutine  SICALC  of  Appendix  B  solves  the  equation  for  by  Newton's 

method.  The  series  frequency  f  is  then  given  by  the  definition  of  i|)  as 

s 


f 

s 


[4<2f2  -  f2] 
H  HJ 
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A  computer  program  named  TMPC1  (see  Appendix  B)  calculates  and  plots  this 
series  frequency  by  calling  SICALC  and  several  other  routines.  The  result 
is  a  plot  of  series  versus  parallel  frequency  with  electron  temperature 
as  a  parameter  for  a  given  set  of  plasma  and  antenna  parameters.  Figure 
18  is  an  illustration  of  such  a  graph  for  a  typical  antenna.  It  can  be 
seen  from  the  figure  that  given  the  series  and  parallel  resonance  frequencies, 
the  electron  temperature  can  be  determined. 

The  parameters  of  antenna  radius,  magnetic  field  and  sheath  radius 
all  have  some  influence  on  the  series  resonance  frequency.  To  illustrate 
the  effect  of  these  parameters,  the  curves  of  Figure  19  were  developed. 

Note  the  decrease  of  series  frequency  as  the  antenna  radius  is  increased 
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and  the  increase  in  series  trequency  as  magnetic  tieid  and  shea’h  widrh 
are  increased  ihus  a  knowledge  of  the  shearh  and  magnetic  tieid  are 
required  m  addition  to  the  antenna  parameters  and  resonance  frequencies 
in  order  to  define  electron  temperature. 

Electron  Temperature  from  Resonance  Frequencies 

Electron  temperature  can  be  found  from  diagrams  similar  to  Figure 
19,  or  they  can  also  be  calculated  directly  It  electron  temperature 
is  desired,  the  resonance  condition  (1C9)  is  used  to  define  <p  as  tollows 

D(aR  )  =  fii/2  +  (l-o2)C  1C  1  in{L/R  )  (lli) 

s  o  s  s 

where  all  the  parameters  on  the  right  hand  side  are  known.  i3ow  <p  is 

a  transcendental  function  of  otR  and  while  uR  cannot  be  expressed  ana- 

s  s 

iytieally,  it  .an  be  found  numerically  trom  the  value  of  4>.  Subroutine 
ETAA  determines  uRg  from  the  value  of  4)  by  a  series  of  simple  iterations. 
The  parameter  j  is  a  function  of  the  electron  temperature  T  and  rhus 

T  =  (4^2R2m/iK)  (f2  -  f2)  ( l—if 2  )  /  ( uR  )2  (115' 

s  s  H  s 

This  temperature  is  calculated  by  a  subroutine  named  TEMPSP-  The  tou' mes 
ETAA  and  TEMPSP  are  used  in  the  mam  line  program  TEP60I.  to  find  electron 
temperature  from  the  series  and  parallel  resonance  c requenc ies •  Again  the 
programs  are  listed  m  Appendix  B. 

Example  ot  Electron  Temperatures  from  Measured  Antenna  Characteristics 
To  illustrate  fhe  measurement  of  electron  temperature,  the  results 
trom  a  Plasma  Frequency  Probe  flown  <.n  an  Aerobee  ro'  ket  AC  1.60  3  were 
analyzed  This  rocket  was  tired  into  an  act i»r  aurora;  and  hence,  some 
large  variations  ot  the  plasma  parameters  occurred,  including  high  eie<tron 
temperatures.  The  complete  experimenr  has  been  described  in  detail 
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elsewhere  [Ulwick,  et  al.,  1965)  and  an  extensive  analysis  oi  the  RF  probe 
experiments  aboard  has  also  been  reported  [Baker  et  al.,  1965).  The 
important  feature  of  this  particular  flight  was  the  independent  measurement 
of  electron  temperature  by  two  other  instruments:  a  planer  retarding 
potential  analyzer  and  an  RF  electron  temperature  probe  (RF  perturbated 
Langmuir  probe,  part  of  the  Resonance  Rectification  Probe). 

Sheath  variations  occurred  with  this  experiment  However,  the  vehicle 
and  antenna  potential  could  be  defined,  and  it  was  found  that  the  antenna 
potential  was  zero  at  certain  times  during  the  flight  [Richards,  1965). 

The  experimentally  determined  potentials  are  shown  in  Figure  20  as  a  runction 
of  the  voltage  step  applied  between  the  vehicle  and  the  antenna.  It  can  be 
seen  that  at  step  5,  the  antenna  potential  is  near  zero  and  that  the  sheath 
must  be  collapsed  for  this  condition.  Thus  the  theory  tor  no  sheath  can  be 
applied  at  this  time  to  determine  electron  temperature. 

Program  TEP601  was  used  to  calculate  electron  cemperature  ior  the  step 
5  case  over  a  limited  region  of  the  flight  wheie  good  results  were  obtained. 
The  results  are  shown  m  Figure  21,  along  with  measurements  irom  the  crher 
two  instruments.  There  is  considerable  scatter  in  the  measurement  of  the 
resonance  frequencies  due  to  telemetry  noise  and  other  expeumen'al  limi¬ 
tations.  Hence  the  temperature  values  show  large  variations  irom  what  must 
be  the  true  values.  In  addition,  the  true  temperature  -alnes  must  alsc  vary 
greatly  as  well,  since  the  antenna  is  in  an  active  auroral  region  character¬ 
ized  by  disturbed  temperature  and  density  variations-  Since  each  independent 
instrument  measures  at  an  independent  time,  some  differences  between  instru¬ 
ments  is  expected  due  solely  to  the  variations  in  the  ;onospheri.  medium 
The  general  agreement  between  all  the  measurements  is  seen  to  be  quite 
good  when  the  above  arguments  are  kept  in  mind. 

In  an  attempt  to  further  observe  the  correlation  between  the  vai^^s 
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techniques,  the  data  under  all  sheath  conditions  was  analyzed  with  the 
collapsed  sheath  theory  and  the  results  were  averaged.  This  average 
curve  is  compared  to  the  other  measurements  in  Figure  22.  Again  reasonable 
correlations  were  obtained.  It  should  be  noted  that  the  experiment  was  not 
designed  to  determine  electron  temperature  and  that  further  experimental 
refinements  with  this  goal  in  mind  would  greatly  improve  the  quality  of  the 
data.  In  general,  it  is  felt  that  the  measurements  of  electron  temperatures 
by  the  resonance  frequencies  could  be  a  very  accurate  technique,  especially 
for  experimental  designs  optimized  for  accurate  determination  of  the  resonance 
frequencies. 


Fig.  22.  Averaged  elec  cron  temperatures  -  AC  3.603 


Chapter  VI 


SUMMARY  AND  CONCLUSIONS 

This  paper  first  discussed  the  plasma  state  and  the  important  parameters 
needed  to  specify  the  character  of  a  plasma.  Then  the  range  of  parameters 
encountered  in  the  ionospheric  plasma  was  examined .  Various  RF  probe  tech¬ 
niques  were  presented  and  a  need  was  shown  for  a  theory  to  describe  the 
plasma  -  antenna  interaction  so  that  measurements  from  these  probes  could 
be  interpreted  correctly. 

Various  theories  have  been  previously  presented  but  while  some  of  these 
were  suitable  for  a  restricted  range  of  plasma  parameters,  none  of  these 
theories  could  account  for  all  the  experimental  results  as  observed  in  data 
from  the  probe  experiments.  Thus  clearly  further  theoretical  work  was  re¬ 
quired  . 

In  the  theoretical  work,  the  antenna  was  assumed  to  be  cylindrical, 
electrically  short,  yet  with  sufficient  length  as  compared  to  the  radius 
so  that  end  effects  could  be  ignored.  The  surrounding  plasma  was  character¬ 
ized  as  a  homogeneous,  compressible  electron  fluid  surrounding  relatively 
fixed  positive  ions  so  that  macro  charge  neutrality  is  maintained.  The 
plasma  properties  are  described  by  the  plasma  frequency,  gyrof requency  and 
the  electron  temperature  and  collision  frequency. 

The  theoretical  approach  considered  the  forces  on  an  element  of  the 
electron  fluid.  In  addition  to  the  usual  acceleration  term,  there  were 
forces  due  to  electric  and  magnetic  fields,  electron  collisions  and  gradients 
of  pressure  in  the  electron  fluid.  The  pressure  gradients  were  related  to 
density  gradients  and  electron  temperature  through  the  ideal  gas  law  and 
equation  of  state.  The  velocity  of  the  electron  fluid  was  found  from  this 
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force  equation.  This  expression  was  then  combined  with  the  continuity  and 
Poisson's  equation  to  produce  a  second  order  partial  differential  equation 
in  terms  of  the  electron  density  variations.  Solutions  to  this  equation 
were  found  to  be  expressible  as  the  well  known  Hankel  functions.  Boundary 
conditions  then  determined  a  unique  solution.  From  a  knowledge  of  che  charge 
density,  Poisson's  equation  was  solved  and  the  potential  in  the  region  was 
then  defined.  The  fields  were  then  easily  determined  and  antenna  impedance 
was  found  as  the  ratio  of  antenna  potential  to  antenna  current  as  derived 
from  the  field  quantities.  The  usual  discontinuous  sheath  model  was  employed 
to  account  for  sheath  effects.  The  result  was  an  impedance  expression  rhat 
included  all  the  important  plasma  parameters  discussed  earlier,  especially 
electron  temperature. 

The  results  of  the  theoretical  calcuation  were  numerically  simulated 
with  a  computer,  and  the  emperical  results  were  compared  to  several  other 
theories.  The  theory  exhibited  good  agreement  with  magneto-ionic  theories 
and  also  with  Balmain's  work  for  a  range  of  parameters  in  which  it  was  valid 
to  make  comparisons.  In  addition,  a  comparison  was  made  to  experimental  re¬ 
sults.  The  agreement  was  good,  at  least  to  within  the  limits  01  the  accuracy 
of  the  experimental  parameters.  In  contrast  to  all  the  other  theories,  n: 
contradictions  were  noted  between  the  theory  and  experimental  results  Thus 
the  theory  appears,  at  least,  to  be  the  best  available  for  interpretation  of 
ionospheric  RF  probe  experiments. 

The  problem  of  determining  plasma  parameters  from  impedance  measurements 
was  considered  next.  First  it  was  demonstrated  that  plasma  frequency,  electron 
temperature,  and  collision  frequency  could  be  obtained  from  the  altitude  and 
a  single  frequency  impedance  measurement.  Next  it  was  shown  how  plasma 
frequency  and  electron  temprature  are  defined  by  the  series  and  parallel 
resonances  of  the  antenna.  Finally  an  example  of  determining  electron 
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temperature  from  experimental  Plasma  Frequency  Probe  results  was  considered. 
The  electron  temperatures  determined  with  the  aid  of  the  theory  were  compared 
to  two  other  independent  measurements  and  good  correlations  were  obtained. 

It  is  believed  that  the  experimental  electron  temperatures  derived  by 
use  of  the  theory  are  the  first  reported,  direct  electron  temperature  mea¬ 
surements  of  the  ionosphere  that  do  not  use  some  variation  of  the  Langmuir 
Probe  Technique  (see  Chapter  I).  The  theory  therefore  makes  possible  a 
valuable  new  method  of  measuring  electron  temperature. 

It  should  be  noted  that  the  theory  was  developed  to  interpret  RF  probe 
results.  Hence,  the  antennas  were  always  very  short  compared  to  the  RF 
v el eng tbs;  and  therefore,  the  electromagnetic  effects  were  dominated  by 
the  plasma  effects.  In  addition,  the  free  space  impedance  of  the  antenna 
systems  could  be  experimentally  determined;  and  thus,  the  theory  was  only 
required  to  account  for  changes  in  the  impedance  due  to  the  plasma. 

In  conclusion,  it  can  be  noted  that  the  theory  successfully  accounts 
for  all  the  experimental  phenomena  observed.  It  agrees  well  with  other 
theories  where  comparisons  are  valid.  Finally,  its  use  makes  possible, 
for  the  first  time,  the  determination  of  electron  temperature  by  an  RF  im¬ 
pedance  probe  technique. 
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APPENDIXES 


APPENDIX  A 


ANTENNA  IMPEDANCE  FROM  IONIC  THEORY 

Isotropic  Media-Capacitor  Model 

If  magnetic  fields  and  thermal  motions  other  than  collisions  are  ignored, 
then  antenna  impedance  can  be  simply  determined  from  ionic  theory. 

The  frequencies  of  interest  are  in  the  range  .1  to  10  MHz.  Thus  the 
shortest  wavelength  involved  is  30  meters.  In  almost  all  cases,  the  an¬ 
tenna  probes  under  consideration  will  be  much  shorter  than  this,  typically 
about  3  meters.  Therefore,  the  antennas  are  always  short  compared  with  the 
operating  wavelengths,  and  hence  may  be  considered  as  capacitors  whose  die¬ 
lectric  reflects  the  properties  of  the  plasma  medium  that  surrounds  the 
probes.  Capacitance  of  an  object  is  defined  as  the  ratio  of  the  total  charge 
stored  on  the  surface  to  the  electric  potential  of  the  object. 


C  =  Q/0 


(Al) 


In  terms  of  the  electric  flux  parameter  D,  the  capacitance  is  from  Gauss's 
Law 
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In  terms  of  the  free  space  capacitance  Cq  and  from  the  relation  valid 

*+■  -► 

in  free  space  D  =  e0E,  it  is  found  that 
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For  the  well  behaved  geometries  encountered  in  practice  (planes,  cylinders, 
etc.),  it  is  permissible  to  consider  both  D  and  E  invariant  in  magnitude 
and  perpendicular  with  respect  to  the  surface.  Thus  these  quantities  can 
be  taken  outside  the  integrals  and  the  following  relation  results 


C  _  £ 

C  =  e 
0  0 

where 

D  «  eE 

The  antenna  impedance  Zflcan  now  be  obtained  from 

Za=  1/jtoC 

Ratcliffe  [1962 )  derives  an  expression  for  e  expressed  by 

£'e° 

where 

2  2 

X  »  <4/1/ 

Z  =  v/w 

=  radian  plasma  frequency 
(o  =  radian  operating  frequency 

v  ■  electron  collision  frequency 

Thus  the  antenna  impedance  Z^  is 

-1 
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or 
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In  this  expression,  it  is  assumed  that  the  antenna  probe  does  not  radiate. 
In  actual  practice  however,  even  a  very  electrically  short  antenna  will 
radiate  some  energy.  This  is  often  described  by  a  radiation  resistance 
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where 


*  total  radiated  power 

I  »  antenna  current  at  terminals, 
a 

Kraus  [1930]  in  (5-55)  gives  the  following  expression  for  Rf 
dipole 
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of  a  short 
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or 

where 

Rf  »  resistance  of  monopole 

L  »  length  of  monopole 

B  =  oj/c  =  2t?A 

c  =  velocity  of  light 

*  =  wavelength 


(A10) 
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In  a  plasma  medium  is  invariant, but  as  shown  above,  e  varies  accord¬ 
ing  to  (A5).  Thus,  the  antenna  impedance  Z ^  arising  from  the  radiation 
can  be  determined  by  substituting  (A5)  into  (A10) .  (K.D.  Baker,  private  com¬ 

munication,  1964)  has  examined  this  problem  and  has  verified  that  both  the  real 
and  imaginary  parts  of  the  radiation  resistance  should  be  included  in  the 
total  expression  for  antenna  impedance.  A  similar  approach  for  the  more 
general  case  of  antenna  impedance  has  been  developed  by  Harrison  and  Benton 
[1959],  where  they  have  let  Imag  [e]  be  represented  by  -o/w  where  o  is 
the  conductivity  of  the  medium. 


The  antenna  resistance  is  thus 


Z 

ar 


(All) 


and  the  total  impedance  Z  is  then 


(A12) 


The  radiation  term  should  be  included  whenever  the  antenna  is  not  extremely 
short  or  when  very  small  changes  of  resistance  are  of  interest. 


Magneto-Ionic  Cold  Plasma 

If  a  steady  magnetic  field  permeates  the  plasma,  the  medium  is  no  longer 
isotropic.  Thus  while  (A3)  still  applies,  the  relation  between  D  and  E  in¬ 
volves  a  tensor  [ e ] .  Since  E  is  always  normal  to  the  antenna  surface,  (A3) 
can  be  expressed  as 
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where  Je  |  *  magnitude  of  E  at  the  antenna  surface, 
s 

But  by  definition  of  [e] 
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.here  t..  =  all  cross  products. 
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Aius  (A13)  becomes 
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At  the  surface  of  the  cylindrical  antenna 


E  *  |E  |  Cos  e 
X  s 


E  =  l  Sin  6 
y  's' 


E  =  0 
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Thus  since  E  *  0  and  the  t  ,  integrals  all  go  to  zero  as  they  contain 
z  lj 

terms  that  integrate  to  zero,  the  equation  becomes 
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It  is  assumed  that  E  is  not  a  function  of  0.  Thus  (A20)  becomes 
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According  to  Appleton-Hartree  Theory  [Ratcliffe,  1962] 
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or  antenna  impedance  is 
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This  expression  is  the  expression  for  antenna  impedance  from  the  magneto¬ 
ionic  theory  of  a  cold  plasma. 


Comparison  to  Balmain  Theory 


The  antenna  impedance  of  (A27)  can  be  directly  compared  to  the  results 
of  Balmain  (see  Equation  10).  As  mentioned  in  the  main  text  of  this  report, 
Balmain's  work  is  probably  the  best  published  to  date  on  cold  plasma  antenna 
impedance.  To  compare  these  two  theories,  they  were  written  up  as  Fortran 
subroutines  known  as  ANTZ  -  BALMAIN  and  ANTZ-MAGNETO- IONIC  THEORY.  These 
routines  were  successively  used  in  a  program  called  SKP4MD  that  calculated 
and  plotted  the  antenna  impedance  as  a  function  of  plasma  frequency  and 
for  fixed  exciting  frequency  (see  Appendix  B) .  The  results  have  shown  very 
good  empirical  agreement  except  when  the  exciting  frequency  is  near  the 
plasma  frequency.  A  typical  sample  of  both  theories  is  shown  in  Figure  A-l 
and  illustrates  the  areas  of  disagreement.  Stone 3  Weber  and  Alexander  [1966] 
have  also  noted  a  theoretical  discrepancy  between  Balmain's  theory  and  the 
theory  of  Herman  and  others  when  the  plasma  frequency  is  in  the  vicinity  of 
the  exciting  frequency.  In  addition,  it  appears  that  their  experimental 
evidence  confirms  that  Balmain's  results  are  not  valid  in  this  region.  Thus 
it  may  well  be  that  the  magneto-ionic  theory  presented  here  is  more  correct 
and  much  easier  to  derive  than  the  theory  of  Balmain.  It  should  also  be 
noted  here  that  this  derivation  is  similar  in  many  respects  to  the  derivation 
of  Crouse  [1964]  made  on  a  rocket  nose  cone  and  thus  should  agree  very  well 
with  his  results  if  the  cone  angle  is  set  to  zero. 

Hot  Plasma 

The  basic  assumptions  and  view  of  the  plasma  in  this  section  are  iden¬ 
tical  to  those  developed  in  Chapter  III.  The  strategy  of  deriving  the  an¬ 
tenna  impedance  differs  in  that  it  is  assumed  that  the  plasma  can  be  rep¬ 
resented  by  an  effective  permittivity  tensor  [e].  Thus  the  same  basic 
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Fig.  A-l.  Comparison  of  Balmain  and  magneto-ionic  antenna  theories. 


a-1: 

equations  are  usea  to  find  [e],  but  instead  of  developing  partial  differential 
equations  as  in  Chapter  III,  this  derivation  derives  the  antenna  impedance 
directly  from  e  by  an  integration  process  similar  to  that  used  above  for  the 
cold  plasma  case. 

Permittivity  [e] 

The  starting  point  in  deriving  [e]  is  the  force  ;quation  (25) • 
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the  force  equation  becomes 
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Now  the  usual  assumption  is  mace  [ Ratcliff e ,  1962,  p.  13]  that  the  velocity 
can  be  represented  by  a  polarization  P  as  follows 
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Thus 
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and  the  force  equation  becomes 
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where 
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Now  tensor  notation  is  introduced. 


Let 


[A] 


*x  3x3y  3x3z 


3x3z  3y3z  3z 


(A37) 


and 


-1 


[a] 


-u  -jY  jY 

J  z  J  y 


JYz  -U  -JYx 


-jY  -jY  -U 
_  y  J  x 


(A38) 


where 


Y  =  Y4  =  Yx  +  Yy+Yz 
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Then  the  force  equation  can  be  written  as 


u  £ 

eX  -  ~  [A] 

O  L 

0) 


+  -1+ 
E  =  [a]  P 


(A40) 


However 


D  =  [e]E  -  e  E  +  P 

o 


(A41) 


-  •* 


.  „t  *  ,.  .*5 


'■*- nr* 


-jr-yw-y  iv^v1 


A-ll 


P  -  [e]E  -  £qU3E 


(A42) 


Thus 


r  ,2e 

Eox--^ 
L  (1) 


[A]  |e  »  [a]  jje]  -  eo[lj]  E 


(A43) 


M  e 

e  X[o] - r2-  to]  [A]  -  [e]  -  e  [I] 

0  2  0 
0) 


(A44) 


[e]  «  eQ  X[o]  +  [I] 


-  [a]  [A]l 

w  J 


(A45) 


hM 

|U(UZ-YZ)J 


(U2-Y2)  (-jUY  -Y  Y  )  ( jUY  -Y  Y  )] 
X  J  z  x  y  J  y  x  z 


(jUY  -Y  Y  )  (U  -Y  )  (-jUY  -Y  Y  ) 
2  y  x  y  J  x  y  z 

(-jUY  -Y  Y  )  (jUY  -Y  Y  )  (U2-Y2) 
y  z  x  J  x  z  y  z 


(A46) 


[See  Ratal  if fe ,  1962;  p.  182-184] 
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where 
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then  the  elements  of  [A]  are 
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Now  let  [e]  =  e 
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Then  the  elements  of  [T 1  are  given  by 
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Therefore  [e]  is 
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where  the  elements  of  [Tl  are  given  above. 

The  permittivity  zensor  represents  the  hot,  magneto-ionic  medium.  It 
should  be  noted  that  it  contains  partial  derivatives  with  respect  to  the 
coordinates.  Thus  the  spatial  form  of  the  £  field  must  be  specified  before 
[e]  is  completely  defined.  This  implies  that  [e]  is  a  function  of  the 


coordinates. 
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Hot  Plasma  Antenna  Impedance 

The  E  field  around  the  cylindrical  antenna  is  given  by 


t  *  E(r) [Cos  (ft)  x  +  Sin  (e)y ] 


(A71) 


and  the  magnetic  field  vector  is  assumed  to  lie  in  the  y-z  plane,  at  an 
angle  y  with  respect  to  the  z  axis  (see  Figure  8  in  the  main  text) . 

Thus 


Y  -  0,  Y  *  Y  Sin  y,  Y  *  Y  Cos  y 
x  ’  y  '  ’  z  ' 


(A72) 


E  =  E  Cos  e 
x  r 


(A73) 


E  Sin  6 
r 


(A74) 


3E  3E  .  .  3E 

_ x  _ x  dr  + x  d6_ 

3x  3r  dx  36  dx 


(A75) 


3E  3E  E 

‘  Sin2  e 

3x  3r  r 


(A76) 


32E  2E 

x  o  „  3 

- j  =  — r  Cos  e 

3x  r 


(A77) 


where 


E 

E  =  — 
r  r 


Similarly 


3  E  2E  2  2  , 

_ x  =  _ o  Cos  e  (,1+Sin  8  ) 

3x3y  3 


Sin  6 


(A78) 
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31  2E 


- Z._ 

3y2 

— |  SinG 
r 

(A79) 

32E  2E  2  .  ,, 2 

y  o  Sin  0  (1+Cos  8) 

(A80) 

3x3y 

3 

r  Cos  6 

All  partial  derivatives  with  respect  to  z  are  zero  as  are  all  partial 

derivatives  of  E  . 

z 

Thus 

32F 

o 

1  X 

Ex  3x2 

Cos  e 

2 

r 

(A81) 

32E 

1  X 

2  Cos  6(1+Sin2  0) 

(A82) 

Ex  3x3y 

r2Sin6 

.  32E 

1 _ z 

Ey  3x3y 

2  Sin  0(  1+Cos2 0) 

2  a 
r  Cos  G 

(A82) 

»2F 

1 _ z 

Ey  3y2 

2  Sin  G 

2 

r 

(A84) 

and 

An  a 

U2  2Coa2e  _  (juYCosv)  2C0SB(1+Sln2e) 
r  r  Sin0 

(A85) 

A12 

Ursine  (1+Cos  0)  _  (jUYCoSy)  2Si|_l 
r  Cose  r 

(A86) 

A21 

jUYCosy  2C°^  6  +  (U2-Y2Sin2y)  .?£08  v 

r  r  Sin  0 

(A87) 

A22  “ 

jUYCosy  -2-S1.nf2^1'fCos  e^+  (U2-Y2Sin2y)  2S1^  6 
r  Cos  0  r2 

(A88) 
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It  follows  that 


11 


U2-(U/X)(U2-Y2)-VU2  2C~ -  +  jVUYCosy 


r2Sin  6 


(A89) 


12 


..,,2  2Sin0(l+Cos2  6)  .  2Sin2  6 

-jUYCosy-VU  - ^ - L  +  jVUYCosy  - 5 — 


2 

r  Cos  0 


(A90) 


T21  *  jUYCosy-V(U2-Y2Sin2y)  2CoS  e^1+Sin  -  jVUYCosy  6  (A91) 


r  Sin0 


Since 


T22  "  u2-Y2sin2Y-(U/X) (U2-Y2)-V(U2-Y2Sin2y)2S^q  6 

r 


-jVUYCosy 


2Sin0(l+Cos  0) 

2  a 
r  Cos0 


(A92) 


-e  X 

[e]  »  - [T] 

U(U  -Y“) 


(A93) 


[e]  is  now  defined.  The  antenna  capacitance  is  given  by  (A13)  and  with  (A14) 
it  is  possible  to  write 


C 

C 


-e  X 
o 

2  2 
U(U  -Y  ) 


l 

i.j 


E.T . .E. 


dS 


E  I  dS 
s 


(A94) 


now  since  E  =0,  and  r  =  R 
z  s 


-X 


2  2 

o  U(U  -Y  ) 


JL 

2tt 


2tt 

E  T  E  E  T  _E  E  T01E  E  T00E 

*  -3-2kr  +  de  W95> 

|E  I2  |E  |2  | E  |2  | E  |2 

o  s  's'  's'  's1 
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by  use  of  the  following  relations 


tdTT 

Sin2  ede  «  n 

0 

•Sir 

2 

Cos  ede  «  tt 

0 


3tt 

Sin4  ede 

0 


3tt 

4 


■dir 

o 


ede 


3tt 

4 


■dir 

Sin26Cos2  6d6  «  I1 
o 


(A96) 


(A97) 


(A98) 


(A99) 


(A100) 


and  by  noting  that  all  the  other  trigonometric  expression  integrate  to  zero, 
it  is  found  that 


'x  11  x  1  2  2  2  2 

- -  dS  *=  -7  [U  (1-3V/2R  )  -  (U/X)(U  -Y  )] 

|E  I  2 

1  s  1 


(A101) 


(•2-rr 


1_ 

2ir 


ExT12Ey  = 

| E  |2 

'  s 1 


-  j  [5U2V/2R2] 

2  S 


(A102) 


r2ir 


1_ 

2  IT 


I E  |2 

'  s 1 


7  [(U2-Y2Sin4y)(5V/2R2)] 

2  S 


(A103) 
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•2tt 

E  T  E 

■  1  —  ,y  •  i  (l/(l-3V/2R2)-(U/X)<U2-Y2)] 

|E  |2  2  s 

0  1  s 


(A104) 


The  capacitance  is  then 


7  -  l-(X/U)(l-(^2/2)Sin2Y)(l-4V/R2)/(l-W2)  (A105) 

L  s 

o 

2 

This  expression  is  identical  to  (A26)  except  for  the  4V/R  term  that 

s 

represents  the  effaces  of  electron  temperature.  Thus  (A105)  is  the  antenna 
impedance  from  magneto-ionic  theory  including  the  effects  of  electron  tem¬ 
perature. 

Simulated  Results 

The  above  expression  (A105)  was  programmed  and  the  results  were  plotted. 

An  example  is  shown  in  Figure  A-2.  If  it  is  compared  to  a  similar  curve  from 
Despain  Theory  as  developed  in  the  main  text,  the  comparison  is  very  poor. 

In  fact ,  the  results  shown  are  physically  so  unlike  other  results  that  they 
cannot  be  accepted  even  as  an  approximation.  Thus  this  theoretical  technique 
cannot  be  used  when  finite  electron  temperatures  are  under  consideration. 

The  exact  reason  for  the  discrepancy  is  unknown  to  the  author;  however, 
it  is  his  opinion  that  the  use  of  a  polarization  vector  P  to  represent  the 
hot  electrons  is  in  error.  This  formulation  does  not  allow  for  electroacoustic 
(longitudinal)  waves  to  exist  in  the  plasma  as  does  the  theory  developed  in  the 
text.  It  may  well  be  that  for 


4V/R2  «  1 
s 
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Fig.  A-2.  Magneto-ionic  finite 
temperature  theory  impedance  plot 


that  this  formulation  would  produce  correct  results.  This  would  be  the  case 
for  either  very  high  plasma  frequencies  or  low  electron  temperatures.  It  is 
not  suitable  for  the  usual  range  of  ionospheric  parameters. 


APPENDIX  B 


COMPUTER  PROGRAMS 


Introduction 

All  of  the  following  programs  were  coded  in  Fortran  IV  and  processed 
on  an  IBM  7044  computer  at  the  University  of  Utah  Computer  Center.  The 
computer  system  includes  a  CalComp  digital  plotter.  Extensive  use  was 
made  of  this  device  to  display  simulated  and  experimental  results.  There 
are  many  programs  that  are  part  of  the  IBM  7044  and  CalComp  system.  These 
system  subroutines  will  not  be  included  here,  but  can  be  obtained  from  the 
University  of  Utah  Computer  Center.  Only  the  special  external  routines 
developed  for  the  work  reported  in  this  paper  will  be  given. 

There  are  four  mainline  programs  and  a  number  of  subroutines,  some 
of  which  are  used  in  more  than  one  of  the  main  line  programs.  Hence  if  a 
given  routine  cannot  be  found  after  the  main  program  in  which  it  is  called, 
then  it  will  either  be  cataloged  under  one  of  the  other  main  line  programs 
or  is  a  subroutine  included  in  the  7044/CalComp  system.  Each  subroutine 
called  by  a  program  can  also  call  other  subroutines  and  read  input  data 
cards.  Thus  before  an  attempt  is  made  to  use  the  programs,  all  subroutines 
should  be  checked  for  additional  call  statements.  The  Fortran  IV  program 
listings  follow  the  descriptions  of  the  routines. 

Main  Program  SWP4MD 

SWP4MD  is  the  main  program  for  plotting  theoretical  antenna  impedance 
as  a  function  of  frequency.  The  theory  is  contained  in  subroutine  ANTZ 
of  which  there  are  four  versions: 


(a)  Magneto-ionic  theory 

(b)  Balmain  theory 


(c)  Whale  theory 

(d)  Despain  theory 
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Any  one  of  these  can  be  used  with  the  main  program.  The  program  produces 
a  tabular  listing  of  the  frequency,  impedance,  and  plot  coordinates  with 
a  heading  block  on  each  page,  it  can  also  write  a  tape  suitable  for  data 


storage  and  external  plotting.  It  produces  a  CalComp  plot  of  impedance 

g 

versus  frequency.  Each  plot  is  8"  x  15",  the  ordinate  ranges  from  -10  to 

g 

+10  ohms  on  a  logarithmic  scale.  The  abscissa  begins  at  zero  and  ends  at 
a  specified  maximum  frequency  (MHz).  There  are  many  options  and  input 


variables.  These  are  explained  below. 


Variable  Name 
NPL0TS 
NRESLV 

NTAPE 

NSKIP 

NPRINT 

NPLT 


MTHY 


ND0 

KGRID 

FMAX 


Function 


Number  of  separate  figures. 

Number  of  separate  data  points  per  curve 
(resolution) . 

Identification  of  tape  drive  unit  that  out¬ 
put  is  written  on.  If  negative  or  greater 
than  4,  no  tape  is  written. 

Number  of  files  skipped  before  data  is 
written  on  tape  unit  NTAPE. 

Ratio  of  total  data  points  plotted  and 
written  on  the  tape  to  the  data  points 
printed  on  the  listing. 

A  branching  constant.  If  set  to  zero,  no 
plot  is  produced.  If  set  to  1,  a  plot  is 
produced  including  labels.  If  set  equal 
to  2,  the  plot  figure  is  labeled,  but  the 
curves  are  not. 

A  branching  constant  such  that  if  set  to  2, 
the  impedance  is  found  as  a  function  of 
plasma  frequency,  otherwise  impedance  as  a 
function  of  operating  frequency  is  determined. 

Number  of  curves  per  figure. 

A  dummy  variable.  May  be  set  to  any  value. 

Maximum  desired  frequency  (MHz)  on  the  ab¬ 
scissa  of  the  figure. 
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Variab le  Name 


Function 


FX  The  plasma  frequency  (MHz)  is  MTHY  is  2, 

otherwise  it  is  the  operating  frequency  (MHz) . 

FH  Electron  gyrof requency  (Hz). 


FNU 


Electron  collision  frequency  (Hz) . 


T  Electron  temperature  (°K). 

R  Antenna  radius  (m) . 


AL  Effective  antenna  length  (m) . 

ANGLE  Angle  between  the  antenna  axis  and  the  magnetic 

field  vector  (degrees). 


RINIT 

XINIT 


M0DE 


Resistance  to  be  added  to  the  antenna  impedance 
and  its  shunt  capacitance. 

Reactance  divided  by  the  applied  frequency  that 
is  to  be  added  to  the  antenna  impedance  if  MODE 
is  1.  Otherwise,  it  is  the  fixed  inductance  in 
series  with  the  antenna  and  its  shunt  capacitance. 

A  branching  constant.  If  M0DE  is  set  to  1, 
jXINIT*F  and  RINIT  are  added  to  the  antenna 
impedance.  If  M0DE  and  MTHY  are  set  to  2,  the 
difference  between  free  space  and  actual  antenna 
impedance  is  produced.  If  MTHY  is  not  2  and  M0DE 
is  set  to  2,  the  actual  antenna  impedance  is  pro¬ 
duced.  If  M0DE  is  set  to  3,  first  the  shunt  cap¬ 
acitance  and  then  the  series  resistance  and  in-’ 
ductance  are  added  to  the  antenna  impedance.  If 
set  to  4,  the  difference  impedance  from  free  im¬ 
pedance  is  plotted,  with  first  the  shunt  capaci¬ 
tance  and  then  the  series  resistance  and  induc¬ 
tance  taken  into  account. 


The  following  cards  are  called  by  the  program: 


Card  No. 

1 

2 

3 

4 


Format 

7110 

13A6 

13A6 


Variables 

NPL0TS ,  NRESLV,  NTAPE,  NSKIP,  NPRINT,  MPL0T,  MTHY 
Major  Title 
Minor  Title 


2I10,2F10. 3, 

3E  12.6,14  ND0,  KGRID,  FMAX,  FX,  RINIT,  XINIT,  MODE 


6E13.6 


FH,  FNU,  T,  R,  AL,  ANGLE 


•ML  .-*  « 


X 
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There  is  one  card  similar  to  card  5  for  every  separate  curve  (the  number 
of  curves  is  given  by  ND0) .  Thus,  if  ND0  should  equal  3,  cards  6  and  7 
would  be  similar  to  card  5.  Each  separate  figure  requires  the  groups  of 
cards  that  follow  card  1.  Thus,  if  NPL0TS  should  be  2,  all  the  cards  after 
card  1  would  be  repeated  for  the  second  figure. 

The  following  subroutines  are  called  by  SWP4MD: 


Subroutine  Name 

Location 

Function 

C0NSTS 

External 

Calculate  physical  constants 

SKPFLS 

System 

Skip  over  tape  files 

ANTZ 

External 

Calculate  antenna  impedance 

IDPL0T 

System 

Generate  plot  identification 

PL0T 

System 

Generate  plot  pen  movements 

SYMBL4 

System 

Generate  plot  symbols 

NUMBER 

System 

Generate  plot  numbers 

FINI 

System 

Generate  plot  ending 

Subroutine  ANTZ 

There  are  several  versions  of  ANTZ,  each  of  which  calculates  antenna 
impedance  according  to  several  different  theories.  The  primary  variables 
are  transferred  through  a  labeled  common  statement  as  follows: 

C0MM0N / IMPC0M/ Z ,  F,  FP,  FNU,  FH,  ANGLE,  T,  R,  AL 


where 


Variable 

Z 

F 

FP 


Function 

Complex  antenna  impedance 
Exciting  frequency 
Plasma  frequency 
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Variable 

FNU 

FH 

ANGLE 

T 

R 

AL 


Function 

Electron  collision  frequency 
Electron  gyrof requency 
Magnetic  aspect  angle  (degrees) 
Electron  temperature 
Antenna  radius  (m) 

Antenna  length  (m) 


Physical  constants  are  transferred  through  another  labeled  common  statement 
/C0NSTA/.  These  constants  can  be  obtained  by  a  CALL  C0NSTS  statement  in 
the  main  line  program.  Further  definitions  of  the  constants  are  included 
under  the  description  of  subroutine  C0NSTS. 

The  Despain  theory  version  of  ANTZ  also  includes  an  additional  labeled 
common  statement  /FPHIC/AR,  PHI  where 


Variable  Function 

AR  Complex  variable  oR  in  the  Despain  theory 

PHI  Complex  variable  0  in  the  Despain  theory 

The  Despain  theory  version  of  ANTZ  calls  subroutine  APHICF  that  calculates 
PHI  from  AR  and  transfers  the  variables  through  the  above  described  com¬ 
mon  statement.  Subroutine  APHICF  must  be  included  along  with  the  main 
line  and  other  programs  when  the  Despain  theory  versions  of  ANTZ  is  used. 

The  Whale  theory  of  ANTZ  requires  an  interpolation  subroutine  INTERP, 

and  this  subroutine  must  be  included  when  the  Whale  theory  version  of  ANTZ 

is  used. 


Subroutine  APHICF 

This  subroutine  calculates  the  function  PHI,  given  by 


-  Www  <•***'  * 
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PHI  -  H  (AR) / (AR*H. (AR) ) 
o  1 

where 

H  »  Hankel  function  of  the  first  kind,  zero  order 
o  * 

■  Hankel  function  of  the  first  kind,  first  order 
and  PHI  and  AR  are  complex  variables  transferred  through  a  labeled  common 
statement  as  follows: 

C0MM0N/APHIC/AR,  PHI 


Subroutine  INTERP 

INTERP  is  a  interpolation  subroutine  that  produces  a  interpolated  value 
of  Y0UT  for  a  given  XIN  from  the  arrays  of  data  XA  and  YA  where  NPTS  is  the 
size  of  the  array. 

Subroutine  CflNSTS 

This  subroutine  calculates  several  physical  constants.  These  constants 
arc  transferred  through  a  labeled  common  statement  as  follows: 

C0MM0N/C0NSTA/PI,  PI2,  EMASS,  Q,  CK,  EO,  G,  UO,  ALPHAK 


The  constants  are  as  follows: 

Constants  Value 


PI 

PIZ 

EMASS 

AMASS 

Q 

CK 


3.14159265 
2.0  *  PI 
9.1083E-31 
EMASS*1836. 12*26. 5 
1. 60206E-19 
1. 38044E-23 


EO 


8.85434E-12 
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Constants 


Value 


G 

UO 

ALPHAK 


3.0 

PI*4.0E-07 

4 . 0*PI*PI*EMASS/ (G*CK) 


Main  Program  TMPC I 

TMPCI  Is  the  main  program  for  calculating  and  plotting  the  series 
resonance  frequency  of  an  antenna  versus  the  parallel  resonance  frequency 
with  electron  temperature  as  a  parameter.  The  program  calls  the  following 


subroutines  in  the  given 

crder: 

Subroutine  Name 

Location 

Function 

INPUT 

External 

Read  in  input  data 

GRID 

External 

Generate  a  plot  grid 

SICALC 

External 

Calculate  <l>  function 

CURVEL 

External 

Plot  array  of  data  as  a  line 

ARR0W 

System 

Generate  a  plot  arrow 

PL0T 

System 

Generate  plot  pen  movements 

NUMBER 

System 

Write  numbers  on  the  plot 

Output 

External 

Generate  output  listing 

FINI 

System 

End  plotting 

EXIT 

System 

End  computing 

Subroutine  INPUT  (TMPCI) 

INPUT  reads  in  the  input  data  for  main  program  TMPCI.  The  following 
cards  are  called  by  this  routine: 

Card  No.  Format  Variables 


1 

2 

3 


2I5.6E11.5  NT,  NFP ,  R,  AL,  FPMAX,  FPMIN,  SHEATH,  FK 

5E16.8  T(l) ,  T(2) ,  T(3),  T(4),  T(5) 

5E16.8  T(6) ,  T(7) ,  T(8),  T(9),  T(10) 
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where 


Variable  Name 
NT 
NFP 
R 
AL 

FPMAX 

FPMIN 

SHEATH 

FH 

T 


Function 

Number  of  curves  per  figure 

Number  of  calculated  points  per  curve 

Antenna  radius  (m) 

Antenna  length  (m) 

Maximum  desired  value  of  parallel  frequency 
(MHz) 

Minimum  desired  value  of  parallel  frequency 
(MHz) 

Desired  size  of  sheath  in  Debye  lengths 
Value  of  electron  gyrofrequency 
Electron  temperature 


Subroutine  SICALC 

SICALC  computes  the  if/  function  of  the  Despain  antenna  impedance  theory. 

It  is  used  in  conjunction  with  main  program  TMPC1  and  the  input  variables 
are  described  under  the  description  of  TMPC1.  Newton’s  method  is  used  in 
an  iteration  process  to  find  tp  from  the  imput  parameters.  The  following 
subroutines  are  called  by  SICALC: 

Subroutine  Name  Location  Funct ion 

C0NSTS  External  Calculation  of  constants 

K  External  Calculation  of  Bessel  Function 


Subroutine  K 

Subroutine  K  is  a  function  subprogram.  The  call  statement  is  K(N,X). 
It  calculates  the  modified  Bessel  function  of  the  second  kind,  of  order  N 
and  argument  X  for  N  ■  0,1.  Subroutine  I  is  called  by  subroutine  K.  X  is 


a  real  variable. 
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Subroutine  I 

Subroutine  I  Is  a  function  subprogram.  The  call  statement  is  I(N,X). 

It  calculates  the  modified  Bessel  function  of  the  first  kind  of  order  N 
and  argument  X  for  N  *  0,1.  X  is  a  real  variable. 

Subroutine  GRID 

Subroutine  GRID  is  a  plot  routine  that  draws  a  system  of  labeled 
coordinates  9"  high  and  of  variable  length.  There  are  four  options  con¬ 
trolled  by  TITLE,  the  first  variable  in  the  call  list.  Option  1  takes 
all  the  control  constants  from  the  call  list.  Option  2  reads  the  control 
constants  from  input  cards  and  puts  the  values  in  the  call  list  and  in  the 
labeled  common  statement  /PLTCMN/.  Option  3  reads  the  constants  and  puts 
them  only  in  /PLTCMN/.  Option  4  reads  no  cards  but  takes  the  variables 
from  /PLTCMN/.  The  call  list  specifies  the  desired  option,  and  the  various 
possibilities  are  given  below. 

Option  Call  List 

1  (TITLE,  SCX,  SCY,  XXB,  XXE,  YYB,  YYE,  NXX,  NYY,  XX15) 

2  (6HVREADR,  SCX,  SCY,  XXB,  XXE,  YYB,  YYE,  NXX,  NYY,  XX15) 

3  (6HbREADb) 

4  (6HC0MM0N) 


The  input  variables  are  explained  below: 

Variable  Name  Function 

TITLE  and  HEAD  These  are  dimentioned  arrays  of  A  format 

data.  These  data  are  the  labels  for  the 
plot  coordinates. 


SCX  and  SCALEX 


The  number  of  major  divisions  along  the 
abscissa. 


SCY  and  SCALEY 


The  number  of  major  divisions  along  the 
ordinate. 


1 
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I 


Variable  Name 

■ 

Funct ion 

1  XXB  and  XB 

A 

The  smallest  value  of  the  variable  to  be 
plotted  along  the  abscissa. 

[  XXE  and  XE 

The  largest  value  of  the  variable  to  be 
plotted  along  the  abscissa. 

|  YYB  and  YB 

The  smallest  value  of  the  variable  to  be 
plotted  along  the  ordinate. 

|  YYE  and  YE 

The  largest  value  of  the  variable  to  be 
plotted  along  the  ordinate. 

^  NXX  and  NX 

The  number  of  significant  figures  in  the 
abscissa  number  labels. 

NYY  and  NY 

F 

The  number  of  significant  figures  in  the 
ordinate  number  labels. 

XXI 5  and  X15 

* 

The  desired  length  of  the  abscissa  of  the 
graph  in  inches. 

Cards  are  read  where 

options 

2  and  3  are  used. 

Card  No. 

Format 

Variables 

r 

13A6 

HEAD  (Label  for  abscissa) 

2 

«r~ 

13A6 

HEAD  (Label  for  ordinate) 

3 

•* 

13A6 

HEAD  (Major  title  of  figure) 

4 

13A6 

HEAD  (Minor  title  of  figure) 

5  6F10.5,2I5,F10.5  SCALEX,  SCALEY,  XB,  XE,  YB 

YE,  NX,  NY,  X15 

No  external  subroutines  are  called.  The  result  is  a  labeled  graph 
X15"  by  9",  divided  into  SCALEX  x  SCALEX  divisions  with  10  minor  divisions 
indicated  along  the  edge  of  the  figure. 


Subroutine  CURVEL 

Curvel  is  a  plot  routine  used  in  conjunction  with  subroutine  GRID. 
Curvel  generates  a  curve  of  line  segments  that  connect  the  array  of  points 


■r 
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given  by  AX,  AY,  and  N.  The  abscissa  values  are  AX,  the  ordinate  values 
are  AY  and  N  is  the  number  of  points  or  array  values.  AX  and  AY  may  be 
dimensioned  at  any  desired  value  in  the  main  program.  The  remaining  argu¬ 
ments  in  the  call  list  are  scaling  constants  provided  by  subroutine  GRID 
and  have  the  same  meaning  as  explained  above. 

Subroutine  OUTPUT  (TMPC1) 

OUTPUT  prints  a  listing  of  the  results  of  subroutine  SICALC  in  a  tabular 
form  with  a  heading  block  on  each  page.  0UTPUT  is  called  by  TMPTC1.  No  sub¬ 
routines  are  called  and  no  cards  are  read  by  0UTPUT.  All  data  is  transferred 
through  the  labeled  common  statement  /INPTS/. 


Main  Program  TEP601 

TEP601  is  the  main  program  for  the  calculation  of  electron  temperature 
from  the  series  and  parallel  resonance  frequencies  of  the  antenna  impedance. 
TEP601  reads  the  input  data,  calls  subroutine  TEMPSP,  and  punches  and  prints 
both  the  input  data  and  calculated  temperature.  The  program  variables  are: 


Variable  Name 


Function 


R 

AL 

PERCNT 


LIMIT 


CC 


S 


Antenna  radius  (m) 

Antenna  length  (m) 

Allowable  error  in  calculation 
of  temperature  (10%  error  would 
be  put  in  as  .10). 

Maximum  number  of  int era t ions 
allowed  in  the  recursive  cal¬ 
culation  of  electron  temperature. 

Value  of  any  coupling  capacitor 
(farads).  If  none  exists,  CC 
should  equal  1.0. 

Sheath  width  in  number  of  Debye 
lengths. 


I 
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Variable  Name 


TSTART 


FMT 

F 


Function 


Iteration  starting  value  of 
electron  temperature  (typ¬ 
ically  1000.0). 

Format  specification  for  main 
input  data  card  deck. 

Series  resonance  frequency 
(MHz). 


r 


FP 


Parallel  resonance  fre¬ 
quency  (MHz) . 


1 


SEC 


Time  in  seconds  or  any  other 
data  identification  variable. 


r 

l 


t T 


An  input  variable  carried 
over  from  the  input  data  to 
the  new  punched  cards  and 
printed  listing. 


The  following  cards  are  called  by  the  program: 

Card  No.  Format  Variables 


1 

3F10.2,I5,E15. 

7.2F10.2 

R,  AL, 

PERCNT ,  LIMIT,  CC,  S,  TSTART 

2 

13A6,  A2 

FMT 

3 

FMT 

F,  FP, 

SF.C,  0T 

4 

FMT 

F,  FP, 

SEC,  0T 

etc. 


The  following  routines  are  called: 


Subroutine  Name 

Location 

Function 

CRLIDr. 

System 

Identify  punched  cards 

TEMPSP 

External 

Calculate  electron  temperature 

CRDIDF 

System 

End  card  punching 
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Subroutine  TEMPSP 

Subroutine  TEMPSP  calculates  the  electron  temperature  from  the  series 
and  parallel  resonance  frequencies  of  an  antenna.  Sheath  effects  can  be 
taken  into  account  by  specifying  the  width  of  the  sheath. 

The  calculation  begins  by  assuming  an  initial  electron  temperature. 
This  is  provided  externally  along  with  the  resonance  frequencies  and  an¬ 
tenna  parameters.  The  Debye  length  is  then  calculated  from  these  values. 
The  Debye  length  is  then  used  to  find  the  capacitance  between  the  antenna 
and  the  outside  edge  of  the  sheath.  The  value  of  4  in  the  Despain  theory 
can  now  be  calculated.  External  subroutine  ETTA  calculates  the  argument 
of  4  from  the  value  of  4-  The  electron  temperature  can  then  be  determined 
from  this  argument  (a'R  in  the  Despain  theory).  This  value  of  temperature 
is  compared  to  the  initial  value,  and  if  it  is  within  the  set  error  limits, 
the  value  of  temperature  is  printed  and  punched  and  a  new  set  of  data  is 
considered.  If  the  temperatures  differ  more  than  the  allowed  error,  the 
initial  test  temperature  is  set  equal  to  the  newly  calculated  temperature, 
and  the  calculation  is  repeated. 

The  variables  for  this  routine  are  transferred  through  the  labeled 
common  block  /KEN/  and  they  are  defined  above  in  connection  with  program 
TEP601. 

Subroutine  ETAA 

ETAA  calculates  the  inverse  function  n  of  the  4  (a'R)  of  the  Despain 
theory  for  real  numbers  only.  Given  the  value  of  (A,  the  subroutine 
determines  a'R  by  a  iteration  process.  A  value  of  a'R  is  assumed.  Thi 
corresponding  value  of  4  is  then  determined.  This  value  is  compared  to 
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the  original  value  of  0,  and  according  to  the  difference  between  the  values, 
the  program  either  returns  a  value  for  a'R  or  makes  a  better  estimate  of 
a'P  and  repeats  the  process  until  the  difference  is  less  than  .00001. 

ETAA  calls  subroutine  K(N,X)  discussed  above  in  connection  with  main 
program  TMPTC1. 


Main  Program  IMPPLT 

IMPPLT  is  the  main  program  for  plotting  antenna  resistance  versus 
reactance  with  electron  temperature  and  plasma  frsquency  as  parameters. 
The  actual  antenna  impedance  is  calculated  in  subroutine  ANTZ.  As  dis¬ 
cussed  above,  several  versions  of  ANTZ  have  been  developed  and  any  of 
these  can  be  used  in  this  main  program.  The  program  variables  are: 


Variable  Name 


Function 


NPL0TS  The  number  of  separate  figures  to  be 

produced - 

NRESLV  The  number  of  calculated  points  per 

curve . 

NTAPE,  NSKIP  Dummy  variables,  may  be  any  integer. 

NPRINT  The  desired  ratio  of  calculated  and 

plotted  points  to  the  number  of  points 
printed  on  the  tabular  listing. 

MPLT  The  desired  ratio  of  calculated  and 

plotted  points  to  the  number  of  points 
indicated  by  a  short  line  segment 
drawn  normal  to  the  curve  and  through 
the  indicated  point. 

MTHY  A  branching  constant  such  that  if  it 

is  set  to  1,  the  electron  collision 
frequency  is  specified  by  the  altitude 
and  electron  temperature,  otherwise 
collision  frequency  is  taken  from  the 
input  data. 
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Variable  Name 

Function 

ND0,  KGRID 

Dummy  variables,  may  be  any  integer. 

FMAX 

The  maximum  value  of  plasma  frequency 
(MHz)  to  be  plotted. 

FX 

The  value  of  the  antenna  exciting 
frequency  (MHz). 

I-:; 

Electron  gyrofrequency  (Hz). 

Fm'U 

Electron  collision  frequency  (Hz). 

R 

Antenna  radius. 

AL 

Antenna  length. 

ANGLE 

Magnetic  aspect  angle  (degrees). 

FNUC0N 

A  constant  that  relates  collision 
frequency  to  altitude  and  electron 
temperature. 

DEN 

Total  particle  density  of  atmosphere 
at  the  given  altitude. 

ALT 

Altitude 

following  cards 

are  called  by  the  program: 

Card  No. 

Format 

Variables 

1 

7110 

KPL0TS ,  NRESLV,  NTAPE,  NSKIP, 
NPRINT,  NPLT ,  MTHY 

2 

2110,  2F10.3 

DN0,  KGRID,  FMAX,  FX 

3 

6E13.6 

FH,  FNU,  T,  R,  AL,  ANGLE 

4 

2E20.8,  F10.3 

FNUC0N ,  DEN,  ALT 

5,6, 7,8,9 

are  called  by 

subroutine  GRID  discussed  above. 

Cards  2  through  9  are  called  once  for  each  individual  figure.  The  number 
of  figures  is  specified  by  NPL0TS. 
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The  following  routines  are  called  by  IMPPLT: 


Subroutine  Name 

Location 

Function 

C0NSTS 

External 

Calculate  physical  constants 

GRID 

External 

Generate  a  plot  grid 

ANTZ 

External 

Calculate  antenna  imped¬ 
ance 

PL0T 

System 

Generate  plot  pen  move¬ 
ments 

FINI 

System 

End  plotting 
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MAIN  PROGRAM  FOR  PLOTTING  ANTENNA 
IMPEDANCE  AS  A  FUNCTION  OF  FREQUENCY 

JIBFTC  SWP4M0 

C  PRuGRAM  for  plciting  antenna  impedance 

01 NCNSIUN  XXI1410),  HEADI50),  XYZISO),  XRI15I0I 
CUPP  LEX  2«ZFIRST,2SERIE  <2 SHUNT ,ZPAR,Z1NTSH 
COHMUN/IMPCOM/Z  ,F ,F P ,FhU,FH, ANGLE • T ,R ,AL 
LOGICAL  NTFRST  ,  UTAPE  ,  FIRST 
UTAPfc  *  .FALSE. 

CALL  CONSTS 

RC  AO  1.  NPLOTS, NRESLV f  NT  APE ,  NSKIP,  NPRINT,  HPLT .  HTHV 
PRINT  1,  NPLOT  S*  NRESLV,  NT  APE  ,  NSKIP,  NPRINT,  MPLT,  MTHV 

1  FORHAT (7110) 

NR  1  *  NRESLV  *1 
NR  X  *  1 500/NRE  SL  V 

IF  I  INTAPE.GE.OI.ANO. (NlAPE.LE.A) I UTAPE*. TRUE. 

IF  I .NOT. UTAPE)  GO  TO  6 

6  REM  I  NO  NTAPE 

7  CALL  SKPFLSINTAPE, NSKIP, 0) 

B  00  1'13  JNP  «1,  NPLOT S 

NTFRST  *  .FALSE. 

READ  9,  ( HE  AO ( J ) , J* 1  ,  I  3 ) 

REAO  9.  I  HE  AO  I  J )  ,  J*  1 A  ,  26 ) 

9  FORMAT  1 1 3A6 ) 

REAO  10, NOO ,KGRI 0  «FHAX,FX«RIN1T,XINIT , C SHUNT, NODE 
PRINT  3,NOO,KGR1D,FMAX,FX,R|NIT,XINIT,CSHUNT,MODE 

10  FORHAT  (2I10,2F10.3,3E12.5,I4) 

3  FORMAT I1H1, 211 0.2F10. 3, 3E 12.6,19) 

FXM  =  F  MAX*  1.0E+06 
FXX  *  Fx»  1.0E*06 
SCALE  *  15.0/FHAX 
SCXCON  »  .  1E*05/ SCALE 

11  FP  *  FXX 

12  F  *  FXX 

13  DO  111  H*l,  NOO 

REAO  14, FH.FNU.T.R.AL, ANGLE 
PRINT  2.  FH,FNU,T,R,AL, ANGLE 

2  FORMAT!  IHI  ,6E 13.6) 

14  FORMAT  1 6E I  3 .6 ) 

FHX*FH»l.0E-06 

FNUX*FNU*1.0E-06 

I F I UTAPE )  WRITtlNTAPE.15)  F , FP.FXM.FH, FNU 
IF (UTAPE )  WRITE (NTAPE, IS)  T , R, AL .ANGLE , SCALE 

15  FORMAT! SE 16.B) 

KPAGC  *  52 

JK  *  NPRINT 

IF (MTHY.EQ.2)  GO  TO  26 

16  00  25  K  *  1 ,  NR1 
AK  *  (K-l ) • (NRX) 

IFPNT  *  AK/ 10. 0 
F  *  AK»SCXCON 

ZSHUNT*CHPLX(0.0, l-1.0/CSMUNT»F»PI2>) 

2SFR IE “CMPLX I RINIT,XINIT»F»PI2) 

FPNT  *  F».1E-0S 
CALL  ANT2 

GO  TO  (551, 552, 553, 552), MOOE 
551  RAA*REAL(2(*RINI T 


■  S"  '  *v 
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XAA» A I MAGI /)*XINIT»F 

552  CO  TO  67 

553  ZPAR«Z*Z$HUNT/I Z+ZSHUNT) 

RAA«RE AL ( ZPAR )  *REAL  I ZSCRI E ) 

XAA>AIMAG( ZPAR I *A INAGl ZSERI E ) 

954  GO  TO  67 

67  RA  ■  RAA 
XA  ■  XAA 
AH  af.O 

IF  I  XA.EQ.  O.O)  XA«1.  OE-37 
IF ( XA.GE. O.OIGO  TO  17 
AH  a  -1.0 
XA  ■  -XA 

17  XT  ■  AL0G10 I XA ) 

IF  IHOOE.EO.il  GO  TO  66 
IFIXT.LT. 1-2. 01)  XT«0.0 
XXIX )*XT»AH 
GO  TO  69 

68  IFIXT.LT. O.OIXT  *0.0 
XX IK )  a  XT* AM/2. 0 

69  IX  ■  lOOO.OaXXIKI 
ARa  *1.0 

IFIRA.EO.O.OIRA-l. OE-37 
IFIRA.GE. O.OIGO  TO  16 
AR  a  -i.o 
RA  ■  -RA 

18  RT  ■  AL0G10IRA) 

IF  ( HOOE.EO. 1 )  GO  TO  70 
IFIRT.LT.I-2. 0))  RT-0.0 
XRIK)aRT»AR 
GO  TO  71 

70  IFIRT.LT. 0.0)  RT «0. 0 
XRIK)  «  RTa AR/2.0 

71  IR  «  IOOO.OaXRIKI 

IFIUTAPE)  WRITE  (NT APE t  5)  IFPNT,  IR  .IX.RAA.XAA,  F 
5  FORMAT  ( 31 I0|2X| 3E16 .8 ) 

JK  »  JK*1 

JJ  »  JK/NPRINT 

IFIJJ.LT. I IGU  TO  25 

XP AGE  -  XP AGE  ♦.  1 

IF (KPAGE.LE.50IG0  TO  2J 

KP AGE  =  I 

PRINT  19 

19  FORMAT  1 1  HI  1 
PRINT  20 

20  FORMAT!  IH  ,’.32H - 


1  - 

2  -  j 

PRINT  4,  FX,  FHX,  FNUX,  ANGLE.  T.R.AL 


4  FOKMATI  2X.4HFP  a,F5.2,14H  HHZ  ,  WH  -,FA.2,  14H  Hl.i  ,  NU  •, 
IF7.4,  I7H  HHZ  ,  ANGLE  «.F5.2,  13H  DEG  ,  T  «,  F7.1,  1AH  DEG  K 
2  ,  R  »,  F7.5,  UH  H  ,  L  •,  F5.2,  2H  H) 

PRINT  37ftH£AO{J)fJal,13) 

37  FORMAT  I  III  1 2X  1 1  3A6.  50H  -  IMPEDANCE  AS  A  FUNCTION  OF  DRIVING  FREQ 
IUENCY  I 
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PRINT  20 
PRINT  20 
PRINT  22 

22  FORMAT  I  1 10M  FREQUENC  V  (MHZ )  RESISTANCE  IOHMSI 

1  (OHMSI  PLOT*  (INCHESI  PLOTR  (INCHES) 

PRINT  20 

23  PRINT  24 1  FPNTfRAAfXAAtXR(K) ,XX(K) 

24  FORMAT ( I H  fF 10 . 3 ,2E25. 8 ,2F20. 3 I 
JK*0 

25  CON ■ INUE 
GO  TO  36 

26  ZSERI£-CMPLX(RINIT,XIN(T»FXX»PI21 
ZSHUNT-CMPLX(0.0,I-1.0/(CSHUNT»FXX«PI2) ) ) 

DO  35  K  «  It  NRI 

AX  •  (X-l !• (NRX! 

IFPNT  -  AK/iO.O 
FP  «  AK-SCXCCN 
FPNT  «  FP-.1E-05 
CALL  ANTZ 

IF(K.EO.i)  ZFIRST-Z-ZSERIE 
GO  TO  (441, 442, 443i444) , MODE 

441  RAA»REAL(Z)+RINI T 
XAm*A(MAG(Z)»X1NIT«FXX 
GO  TO  87 

442  RAA«REAL(Z)-REAL(ZF(RST) 

XAA- At MAG(Z)-AIMAGIZFIHST) 

GO  TO  87 

443  ZPAR«Z»ZSHUNT/(Z*ZSHUNT) 

RAA  »RE AL ( ZPAR) 4 REAL ( 2SERI t  > 
XAA«AIMAG(ZPAR)*AIMAG(ZSERIE ) 

GO  TO  87 

444  ZPAR>Z«ZSHUNT/(Z«ZSHUNT) 

IF(K.EO.l)  ZINTSH-ZPAR 

RAA  -REAL (ZPAR)-( REAL (ZINTSH)-REALIZSERl Ell 
XAA-AIMAG (ZPAR )- ( AI MAG ( ZI NTSM J- A I HAG( Z SERI E ) ) 

87  RA  *  RAA 
XA  *  XAA 
AM  »1.0 

IF ( XA.EU.O.O)  XA-l.OE-37 
IF(XA.GE.O.O)GO  TO  27 
AM  *  -1.0 
XA  *  -XA 

27  XT  *  ALOGIO(XA) 

IF (MODE. EO. I)  GO  TO  78 
IFIXT.LT. (-2. 011  XT-0.0 
XXIKI-XT-AM 
GO  TO  79 

78  IF  (XT.LT.O.O)Xl  -0.0 
XX ( K )  -  XT-AM/2. 0 

79  IX  »  1 000. 0-XX (K ) 

AR-  ♦1.0 

IF (RA.EO.O.OIRA-I.OE-37 
IF (RA.GE.O.OIGU  TO  28 
AR  «  -1.0 
RA  »  -RA 


REACTANCE 
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28  AT  -  AL0&10IRA) 

IFINOOE.EQ.il  CO  TO  80 
IFIRT.LT. (-2. 01)  RT-0.0 
XAIK)*Rf»AR 
CO  TO  81 

ao  IFIRT.LT. 0.01  RT-0.0 
XRIKI  ■  &T-AA/2.0 
at  IR  -  1000.0-XAIKI 

IFIUTAPE)  WRITE  INTAPE. SI  IFPNT, IR, IX,RAA,XAA,FP 

JK  -  JK+1 

JJ  »  JK/NPRINT 

IFIJJ.LT. I  ICO  TO  3S 

APACE  »  KPAGE  ♦  1 

if  I APACE  .Lb. SOI  GO  TO  33 

KPAGE  *  ! 

PRINT  19 
PRINT  20 

PRINT  2I.IHEA0IJI ,J*i ,131 

21  FORNATI 1H  ,  2X.13A6.4BH  -  IMPEDANCE  AS  A  FUNCTION  OF  PLASMA  FREQUE 
INC  V  I 

PRINT  666,  FX,  FHX,  FNUX,  ANCLE,  T.R.AL 
666  FORMAT!  2X.4HF0  -.F5.2.14H  MHZ  ,  MH  -,F4.2,  14H  MHZ  ,  NU  *, 
1F7.4,  17H  MHZ  ,  ANGLE  -,F5.2,  13H  DEG  ,  T  •,  FT.l,  14H  OEC  K 
2  ,  R  «,  F7.S,  IIH  M  ,  L  »,  F5.2,  2H  Ml 

PRINT  20 
PRINT  20 
PRINT  22 
PRINT  20 

33  PRINT  24,  FPNT,RAA,XAA, XRIKI, XXIK1 
JK-0 

35  CONTINUE 

36  IF IUTAPE IENO  FILE  NT APE 
IFIMPLT.EQ.OIGO  TO  Ill 
IFINTFRSTIGQ  TO  106 
SCFP  -  SCALE-FX 
M1B-ISCALE-FX-1. 51/2.0 
WIE-W18+1.5 

M2B  -  (13.5*SCALE*FX)/2.0 

U3B  -  U2B  ♦  .Of 

AIB-W1B+1.55 

A2B-M2B+1.6 

TVS  ■  0.0 

7YM  -  0.0 

KI  *  K 

CALL  IDPLOT 

CALL  PLOT  12. 0,5. 3, -3) 

FPP«FX»SCALE 
L  ■  150. /SCALE 
LL  -  L/10 
00100J-1.L 
AJ  -  J 

AJ  -  AJ-SCALE/10.0 
CALL  PLOT  I AJ, 0.0,2 1 
CALL  PLOT  IAJ,*0.1,3I 
CALL  PLOT  I AJ,-0. 1,2 1 
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100  CALL  PLOT  IAJ.0.0,  31 
IFIMPLT.NE.lt  60  TO  101 
CALL  PLOT  If PP,*0.1 f  31 
CALL  PLOT  IFPP,-0.1,2I 

CALL  SYM8L4IFPP-.05  » »0. 25 »0.15»1HF, 0.0,11 
CALL  SVNBL4IFPP4.02,40.2,0.1,lHP,0.0,il 

101  DO  1 02 J  »1,LL 
AA  *  LLM-J 
BJ  »  AA* SCALE 

CALL  PL0T|BJ,*.15,3I 
CALL  PL0TI8J,-. 15,21 

102  CALL  NUMB£R(8J-0.17,-0.30,0.15,AA,0.0,lt 
IFIMPLT.NE.lt  GO  TO  103 

CALL  SYMBL4I 6.35.-.65.0.2 ,12HFREQUENCV-MC,0.0, 121 

103  CALL  PLOT  I0.0,-4.0,31 
IF  I  MODE  .NE.ll  GO  TO  222 
00  1 05  J  >1.17 

JC  «  J-9 
CJ  ■  JC 
CY  «  C J/2.0 
Cl  «  CJ 

IFICJ.LT. 0.01  Cl  *  -CJ 
CALL  plot i  o.o,  cy  ,21 
CALL  PLOT  I *0. 1 «CY  ,21 
CX  •  *10.0 

IFICJ.LT. 0.01  GO  TO  114 

CALL  NUMBERI-O.l.CY-O.l  ,.15,CX  ,90.0,01 

104  CALL  NUHB£RI-.25,CY40.1,0.1,CZ.90.0,01 

105  CALL  PLOT  I0.0.CY.3I 
GO  TO  406 

222  DO  405  JJJ-1,13 
JC-JJJ-7 
CJ*JC 

CY-CJ*. 6666666b 
C2-CJ 

IFICJ.LT. O.Ot  CZ— CJ 
CZ-CZ-2.0 

CALL  PLOTIO.O,  CY  ,21 
CALL  PLOTIaO.I.CY  ,21 
CX  •  tlO.O 

IFICJ.LT. 0.01  GO  .0  114 

CALL  NUMBERI-O.l.CY-O.l  ,.15,CX  ,90.0,01 

223  CALL  NUMBER I -. 25,CYtO. 1,0.1,CZ,90.0,0) 

405  CALL  PLOT  I0.0,CY,3I 

406  CALL  SVMBL4 1-.5.-1.8 ,.2 .24HANTENNA  IMPEDANCE  I  OHMS  I ,90.0, 24 1 
CALL  SYMBL4I0.0  ,-4.5  ,  0.2 .HEAD 1 1 ) ,0.0, 78 1 

CALL  SYMBL4I 3. 35,-4. 7  ,  0. 1 ,HEAOI 14 1 ,0.0,78  I 

CALL  PLOT  1 40.0 ,-4.0, 3 1 

CALL  PL0TI15. 0,-4. 0,21 

DO  118  J  -  1  •  17 

JC  -  J  -  9 

CJ  ■  JC 

CY  -  C J/2.0 

CALL  PLOTI 15.0,CY«2I 

CALL  PLOTI 14. 9,CY,2I 
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IIS  CALL  SLOT! 15.0 *CYr3l 
CALL  SLOT!  0.0, ♦4.0, 2) 

106  NTFRST  -..TRUE. 

F1KST  -.TRUE. 

DO  109  K  -  1,  NR1 
AK  >  IK-l  )»NfiX 
AK  -  AK/100.0 
XT  -  XX(KI 

IF ( XY.GT.4.0)  XV  -  4.0 
IFIXY.LT. (-4. 0)1  XY  *  -4.0 
IF  (  FIRST)  CALL  PL0TtAK,XY,3) 

FIRST  -  .FALSE. 

CALL  PL0T(AX,XV,2) 

IFtMPLT.NE.il  GO  TO  109 
IFIAK.LT.W1E)  GO  TO  107 
lFtAK.GT.AlS)  GO  TO  lof 
IFtXY.GT.YVSIYVS  •  XV 
IF  t XV.LT. VVHI YVN  -  XY 

107  IF (AK.LT.W2B)  GO  TO  108 
lFtAK.GT.W3S)  GO  TO  108 
IFtXY.LT.XYSIXYS  >  XY 
IF ( XY.GT.XYNI XYM  -  XY 

108  IFIKI.EQ.O)  GO  TO  109 

IFt ( XY.GE .0.0 ) . ANO. ( XYY.LE.O.O) I  GO  TO  115 
XYV  -  XY 

109  CONTINUE 
FIRST  -  .TRUE. 

00  110  K-1.NR1 
AK  -  (K-1MNRX 
AK  -  AK/100.0 
RR  -  XR(K) 

IF (RR.GT.4.0)  RR  -  4.0 
IFtRR.LT. (-4. 01)  RR  •  -4.0 
IFt  FIRST  I  CALL  PLOT  IAK,RR,3> 

FIRST  -  .FALSE. 

CALL  PLOT  (AK.RR.2) 

IF IMPLT.NE. 1 1 GO  TO  110 
IFtAK.LT.W28)  GO  TO  110 
IFtAK.GT.W38)  GO  TO  110 
IF IRR.GT.RYSIRYS  -  RR 
IF IRR.LT.RYNIRYN  -  RR 

110  CONTINUE 

111  CONTINUE 

IF  (MTHY.NE.  1  )G0  TO  112 
IFtMPLT.EQ.OI  GO  TO  113 
IFIMPLT.NE.1IG0  TO  112 
IFtSCFP.LT.l.OIGO  TO  112 
IFIYYS.GT.3.0IYSW  -  YYM  -  .5 
IFt YYS.LE. 3.0) YSW  -  YYS  ♦  .5 
IFt YYS.LT. 0.0)YSW  •  YYM  -  .5 
IFtVYS.LT.t-3.0) ) YSW  -  YYS  ♦  .5 

CALL  SVMBL4( WIB, YSW»0. 1, 17HANTENNA  REACTANCE, 0.0, 17 ) 
IF  1 1  SCFP-15.0) .GT. (-1.6)1  GO  TO  112 
IFt XYS.LT. t-3.0) ) XSW  -  XYM  ♦  .5 
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XVK  ♦  .5 
XVS  -  .5 
RVN  -  .5 
XVS  ♦  .5 
RYN  -  .5 
IF(RYS.LE.(-3.0))RSW  -  XVS  ♦  .5 

CALL  SYN0L4(  M2B,XSK,0.1,17KANTf:NNA  MACTANCE.0.0,17) 
CALL  SVML4(M2B,RSM,0.1,1BHANTENNA  RESISTANCE, 0.0«1>) 


IFUYS.CE.C  O.ODXSM 
IF (XYS.GE •  (  3.CHXSM 
IFUYS.CT.I  3.0IIRSH 
IFCRYS.LE. (  3.0) )R$M 
IFIRVS.LE.I  0.0) )RSM 


112  CONTINUE 

CALL  FLOT  (18.0, -5.3, >3) 

113  CONTINUE 

I F ( UTAFE ) RE MI NONT AFE 
IF(NFLT.EQ.O)  STOF 
CALI.  FINI 
STOk* 

114  CX  «  -10.0 

CALL  NUNBERI-O. 1,CV-. 22, .15 »CX ,40.0,0) 
IMHOOE.NE.  1 )  CO  TO  223 
GO  TO  104 

115  CALL  SVNBL4( AX-.03,. 17,.07,  3 XT  -,40.0,3) 


TT  -  T 

IF(T.EQ.O.O)  TT  -  1.0E-36 
KAT  -  ALOCIO(TT) 

TAX  -  KAT  *7 

TX  -.17  ♦  TAX*. 05636363 

CALL  SVNBL4IAK-. 10, TX-.05*. 15, 1H. ,90.0,1) 

CALL  SYN8L4  (  AX  -  .03  ,  TX  ,  .07  ,  INK,  90.0,1) 
CALL  NUNBERUK-.03,. 41,. 07, T, 90.0,0) 

XI  -  0 

CALL  PLOT  (AK,XY,3) 

GO  TO  109 
ENO 
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SUBROUTINE  TO  CALCULATE  ANTENNA 
IMPEDANCE  ACCORDING  TO  MAGNETO¬ 
IONIC  THEORY 


O  tlBFTC  ANTI  MAGNETO  -  IONIC  THEORY 

C  FUNCTION  SUBPROGRAM  FOR  C ALCLAT I  ON  OF  ANTENNA  IMPEDANCE 
SUBROUTINE  ANTZ 
COMPLEX  U#U2,  YUtOtZ 

CONHON/CONSTA/P1, PI  2, EHASS, AMASS, 0,C*,EO,G,UO, ALPHAS 

COMMON/ INPCON/2,F,FP,FNUtFH, ANGLE, T.R.AL 

CL-ALOGIAL/R) 

lF(F.LE.O.O)  GO  TO  1 

MF»PI2«F 

CO«PI2»EO» AL/CL 

IF  (FP.LE.O.O)  GO  TO  2 

F2*F««2 

FP2-FP»»2 

X-FP2/F2 

Y-FH/F 

Z-FNU/WF 

U«(1.0,0.0)-(0.0,1.0I»Z 

Y2«V«Y 

U2-U»U 

VU-V2/U2 

VR*(  (G*CMTI/(2.0*EMA$S*R*R)*FP2) 

0»(0.0,1.01»WF»CO»( U*( YU-l.O)*! l.-YR)»! X-( X»YU»SIN( ANGLE )**2) /2» )  I 
IF (REAL (0 1 .EQ. O.  O.ANO. A1MAG’ OI.EQ.O.OI  GO  TO  U 
Z-U»( YU-1. 01/0 
RETURN 

1  Z»(0.0,-. 1 E ♦ 37 1 

8FTII0M 

2  Z  —  (0.0,1. OI/(WF»C01 
RETURN 

11  PRINT  12,YU,UtX 

12  FORMAT! 1 AH  ERROR  IN  ANTZ  ,  5E16.8) 

RETURN 

END 
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SUBROUTINE  ANTZ  ACCORDING  TO 
BALMAIN  THEORY 


« IBFTC  AUTZ  L 1ST 

C  FUNGI  ION  SUBPROGRAM  FOR  CALCULATION  OF  ANIENNA  IMPEDANCE 
SUBMOUfINE  ANTZ 

COMPLFX  U,AZ,EXO,EKP,AA2,AA,FM,FR  , AR , J , TT , CSQKT , CL OG, TNHAC, A A  1 , Z 
COhMON/CONSTA/PI ,PI2,EMASS,AMASS,O.CK,EO,G,UO, ALPHA* 

COMMON/ I MPCOM /Z  *  f  *FP  i  F  NUi  FH  .ANGLE  , I , R, AL 
WF  *  P 1 2»F 
F2  *  F»F 
f P2  *  FP»FP 
X  *  FP2/F2 
V  *  FM/F 
1  *  FNU/MF 

U  «  II. 0,0.0)  -I0.0,1.0)»Z 

EKG  «  I 1.0, C.O)  -  X/U 

FKP  *  I 1.0,C.0)-X»U/IU»U-V»Y) 

AA2  «  EKP/EKU 

tt  *10.0,0.0) 

IFIT.LT. 1.0) 00  TO  2 
A2  *  I  F2» ALPMAK/T  )•  I  X-U) 

AR  *  R.CSJRT. A2) 

TT  *  IfcKO  -  I  1.0,0.  0)1/112.0,0.  0)*AR) 

2  AA  «  CS0RTIAA2) 

RAO  *  P  I • ANGL  I:/ 1 80. 

SN2  *  S INI  RAO  )»«2 
CN2  *  COSI  RAO  )»»2 

FM  .  CMPLX I SN? , 0. 0)  ♦  AA2»CMPLX|CN2,0.0) 

FR  *CSORTIFM) 

AA  1  *  IAA  +  FR  )/(( 2.0,0. 0)»FM) 

TNMAC  =  TT 

IFIIREALIAAM.tO.O.U).  AND.  I A  I  MAO  I  AA  1  )  .EO.O.ODGO  TO  3 
TNMAG  *  -ClOGIAAl)  +TT 

1  Z»AA.|  |  ALOGI  AL/R  )-l  .  C )  ♦  TNMAG )  / 1 1  10.  0 ,1 . 0 )  »E KP«FR)  •  I  P 1 2 •WF«EO» AL )  ) 
RETURN 
fcNO 


aTssasar’1 
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SUBROUTINE  ANTZ  ACCORDING  TO 
WHALE  THEORY 


SUOROUTINt  ANTZ 
COMPLEX  Z,ARC,A2C 

COMMON/CONSTA/PI  ,P12  .EMASS, AMASS, Q,CK,EO,G, HO,  ALi'HAK 
COMMON /IMP COM/  Z.F, FP,FNU,FH,ANGLfc,T,R, AL 
DIMENSION  XROUOOl,  BKRO(IOO) 

LOGICAL  FIRST 
REAL  LAND 

IFIFIRST!  GO  TO  11 
FIRST*. TROF. 


DA 

TA  8KR0/.50, .52, .535 

,.55 

,.565 

9 

.58, 

.60,  . 

61,. 62, .635,. 645, 

.655 

,.66 

15, 

.67, 

.675, .685, .695,. 

70,. 

705,. 

7 

10,. 

715,. 

72, .725, .728, .73, 

.733 

,.73 

25, 

.738 

,.740, .745, .748, 

.750 

,.755 

.758 

,.759 

,.760, .762, .763,. 

764, 

.765 

3,. 

77,. 

772,. 77,. 765,. 755,. 75,. 74 

.735 

,.72, 

•  7159*709  *695 9  *665 9*6 

5, .6 

43. 

.615 

9 • 600 » • 5A0  *  *565  » 

.555 

,.545 

.51, 

.52,. 

51, .505, .495, .485 

,.47 

5, .4 

565 

,.455, .450,. 440,. 435 

,.430,.  42 

.415 

,.410 

,.405, .402, .395,. 

390, 

.385 

6,. 

380, 

.375, .370, .368,. 

365, 

.360, 

355, 

.350, 

.340, .325, .315,. 305,. 

295, 

7.290,. 

?80, .275, .270, .265/ 

DATA  XR0/.10,.I1,.12,. 

13,. 

14, .15 

,.16 

, .  I  7, 

.18, .19, .20, .21,. 

22,. 

23,. 

124 

,.25 

,.26, .27, .28, .29 

,.30 

9*319 

• 

32,. 

33, .34, .35, .36, .37, .38 

,.39 

,.40 

2.. 

41,. 

42, .43, .44, .45,. 

46 1  • 

47, .48 

,.50 

,.55, 

*609*659*709  *759  • 

80,. 

C5,. 

4. 2. 5, 2. 6. 2.  7,2.8 ,2. 9, 3. 0,3. 1,3.2,;- 


3, 3. 4 ,3. 5, 3. 6, 3. 7,3. 8, 3. 9,4.  0,4 


).lv4.2i4.3t4.4i4.$t4.6i4.7|4»8v4.,|S.OiS.$t6.0i6.$i7.0t7.$i8.0i8.5 
6,9.0,9.5.10.0  / 

11  CL  *  ALOGUL/R  I 

IF (F .IE. 0.0)  GO  TO  1 
WF  -PI2*F 
CO  *  PI2»E0»AL/CL 
XO  »  -1 .0/ I MF«C0 I 
IFIFP.LE.O.OI  GO  TO  2 
F2  *  F**2 


FP2*  FP*»2 
X  *  FP2/F2 
I F ( X.EO. 1.0)  GO  TO  1 
ZX  *  XO/II.O  -  XI 
IFIT.LE.O.OI  GO  TO  3 
V  *  SORT(G»CK»T /EM ASS  I 
LAMD  *625. Q» SORT  1 T) /FP 
RO  *R  ♦  4.0»LAMD 
C2  *(C0/ALI»»2 
IF ( X.GE . I .0 1  GO  TO  3 
A2  *11.0  -X I »F P2  •  ALPHAK/T 
AR  *  SORT ( A2 I »R0 


2.14  *4 . 0»  (  ZX» • 2 1 

CALL  INTFRP  IXRO,BXRO,AR,WHETA,lOO) 

PVI*PI2«EO*RO»SORT( l.O-XI /(X»V»C2»AL«MHETA»*2I 

PVI 2*  PV  I  »«2 

IFIZX4.GT.PVI2IGO  TO  3 

KI*  .5* IPVI  ♦  SORT! PVI 2-2X4)  I 

R2*  .5* IPVI  -  SORT ( PVI 2-ZX4 I  I 

ZR  *R1 

IF(RI.GT.R2)ZR*R2 
/  *  CHPLXIZR.ZXI 
RE  TORN 


1  Z  *  CMPLXIO.O,-IC*371 
RE  1 ORN 

2  Z  *  CMPLXI O.O.XOI 
RE  TORN 

1  Z  *  CHPLXIO.O.ZX) 

RE  TORN 
tNO 
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8UM0UTME  TO  CALCULATE  ANTENNA 
IMNCCANCE  ACCONOtNC  TO  OESPAM  THEONY 


tllFTC  ANTZ  DESPAIN  THE  OUT  12/0«/65 

C  SUtROUTINE  TO  CALCULATE  ANTENNA  1MPEOANCE  AC  CORO  INC  TO  DETAIN  THEORY 
SUtROUTINE  ANTZ 

COMPLEX  U.A2.U2.S2. A2U2DU.AR.CSQAT.PH1.Z.XO 

COX AON/ CONST A/P J .PI  2.EHASS. AMASS, QtCK.EO,  C.UO.ALPHAK 

CON  HON/ INf’':ON/Z,F,FP.FNU,FH, ANCLE,  T.R.AL 

COMMON  /FPHIC/AR.PH1 

CL-ALOCIAL/RJ 

IF(F.LE.O.O)  CO  TO  1 

NF-PI2»F 

CO-P12*EO«AL/CL 

XQa-IO.O.  l.OI/IWFKOI 

IF  (FP.LE.O.OJ  CO  TO  2 

F2-F»*2 

FP2-FP««2 

X-FP2/F2 

Y-FH/F 

Z-FNU/MF 

U-l l.O.O.OI-IO.O.l. OI«Z 

Y2»Y«Y 

U2-U*U 

S2-IU2-V2I/IU*XI 
1FIT.LE.O.OI  CO  TO  3 
A2* I S2- I l.O.O.OI l*FF2»ALFMAR/T 
AR-R«CSQRTIA2I 
3  IF(REALIARI)  6.7.E 
«  1FIA1MACIARIIU. 17.10 
7  IF  I  AIN AC(  AP  I 19.17.10 
•  IFI AIHACIARI 19.9*11 
3  AR— AR 

17  IFIAIHAGIARI.LE.O.OI  AR-AR*(0.0,.1E-16I 

10  CALL  APHICF 

Z-X0»(PHI/CL-S2)/((  1.0.0. 0I-S2I 
RETURN 

1  Z-I0.0.-.1EO7I 
RETURN 

z-xo 

return 

3  Z~X0*S2/I(1.0,0.0)-S2I 
RETURN 

11  PRINT  12.AR.FH.FINS.ALPUAK.T.R 

12  FORNATt  1AH  ERROR  IN  ANTZ  .  7E16.A) 

RETURN 

ENO 
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SUBROUTINE  FOR  CALCULATION  OF 
PHYSICAL  CONSTANTS 


MIFTC  CONSIS 

SUM  OUT  1  ME  COMSTS 

COMMON  /CONSTA/MI,MI2,ENASS,ANASS,Q,CK.EO,G,  UO.ALMMAR 

Ml  -  3. 14194269 

MI2  -2.001 

ENASS  ■  9. 1083E-31 

AMASS  ■  EMMS- 183*. 12*26.3 

Q  ■  1 . 60206E- 1 9 

CR  -  1.38044E-23 

EO  -  S.SS434E-12 

C  -  3.0 

UO  -  Ml  •A.OE-OT 

ALMHAK  -  4.0OIOI-ENASS/IG-CK) 

RETURN 

ENO 


INTERPOLATION  SUBROUTINE 


SUBROUTINE  I NTFRMIXA, YA, XIN. VOUT.NPTS ) 

DIMENSION  XAIIOM.YACIOO) 

I»I 

1  IFIXA(I«I).GE.XIN.AND.XAII).IT.XIN)  GO  TO  2 
IF!  I.GE.INPTS-  2) )  GO  TO  2 

;•!♦! 

GC  TO  l 

2  YCUT-IIXII  ~X  Al  1*1  >)*(XIN-XAt  1*21) -YAm/C  <XAI!)-XAI  !♦!))•(  XAf  I  )>XA 
II  I«2)m«IIXIN-XA(I  ))*(XIN-XAIIM2n»YAII«l)/l(XA(I«l)>XA(n;*(XA(I 
2+H-XAI  M21)  >)■»((  XI F9-XAT  I) ) •( X1N-XA! I* 1) ) -YAI I«2> /(  IXA(I«2)-XAI  1M 
3)>«<XAI!*2)-XA<!)))) 

RETURN 

ENO 
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SUBROUTINE  FOR  CALCULATION  OF 
PHI  FUNCTION 

SIBFTC  APHICF  DECK 

SUBROUTINE  APHICF 

LOGICAL  NTFftST.VALl.VAL2.VAL3.VAL4 

DOUBLE  PRECISION  XR.XI.PHIR.PHI I .PI.PI I. PI I 2.CP I2R.CPI2 I .P14I «P 13 
14.CP4R.CP4I.CP3ft.CP31. Jft.Jl iCNULTft.CMULT l.COVOR.OCVOI. AK.Al.Bft.Bl 
2.VR. VI • V2R.V21 .BOft.BOl .Blft. B1 I. BJOft ,BJO 1 .BJlft.BJl l.BNOft.BNOI. BNlft* 
3BN1I .TJOft.TJOI.TJlft.TJlI .TNOft.TNO 1 .TNlft .TNI  I . F.St 60RS* BO IS.Blft S.Bl 
4IS.F1.SH.BN1RS.BN1IS.BJ lft S.BJ US. BFft.BF 1.0LOC.OATAN2.H0R.HO1 .Hlft.H 
S1I .XBft, XBI . XBGR.XBQI .P9R.P0I .OOft.OOI.Plft.Pl I.Qlft.Ql I.TPOft.TPOl ,TQD 
6ft.TOOI.TPlft.TPlI.T01ft.T01 1 .0.02 .0H2. 04. 040. 04. 00. 020.  F2.XBlft.XB 

71 1 . XBOI ft .XBOI I .  ftPO.ROO.ftPl.ftOl.POftS.POIS.QOftS.QOlS.PlftS.PllS.Q 
SIRS. 01  IS. 00ft SI .OOlSl.OlftSl.QllSl 
CGHPLEX  Aft. PHI 
COHHON  /FPHIC/Aft.PHI 
Xft  -REAL  I  Aft  | 

XI 'A  I  MAGI  Aft  I 
VAL1-. FALSE. 

VAL  2  s.  FALSE. 

VAL3-. FALSE. 

VAL4-. FALSE. 

IF  (NTFRST )  GO  TO  1 
NTFRST-.TRUE. 

PI  -  3. 1415926535897932 

PII-1. 000/PI 

PII2-16. 000-PI  1 

CPI2R-PII2 

CPI2I-0.000 

P 1 41  -PI/4.000 

PI  34-PI  41 *3.000 

CP4R-P14I 

CP4I -0.000 

CP3R-PI34 

CP  3 1 -0.000 

Jft -0.000 

JI -1.000 

CHUL  T«  (AK.Al.Bft.Bl)-  I  AMBR-AI-BI  > 

CNULTI (AK.Al.Bft.Bl)-! BR-AI -AK-B 1 1 
COVOft ( AK.AI ■ 8ft .81 )-( AK-BR-AI-BI I / IBR-BR-BI-BI I 
COVOI (AK.AI .Bft.BD- (AI-BK-AK-BI 1/ (Bft-Bft+Bl-Bl 1 
1  Vft-CHULTft I > SOO  .O.OOO.XR.Xl) 

VI-CHULTlI.SOO.O.OOO,Xft.XI) 

K-l 

IF ( (Vft.EO. 0.000) .AND. IYI.EQ.0.0) )  K-2 
GO  TO  ( 3.4) ,K 
4  YR-.  10-17 
VI -0.0 

3  V2R-CNULTRI Y«, YI i Yft , VI I 
V2I-CNULTI I VR.VI .VR. VI ) 

IFIIXR**2+X1**2) «GE .64.000)  GOTO  13 

BOft-l.OOO 

BOI-O.OOO 

BIR-BOft 

81 1 -B  01 

BJOft -B  Oft 

BJOI-BOI 


1 
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BJ1R-80R 

BJ1I-BOI 

BNOR-O.OOO 

BNOI -0.000 

BN1R-BOR 

BNII-BOI 

T  JOB -BOR 

TJOI-BOl 

TJ1R-BOR 

TJ1I-BOI 

TNOR-BNOR 

TN0I-6N0I 

TN1R-BOR 

TNI I-BOI 

F-1.000 

S -0.000 

DO  9  K-t.100 

60RS-B0R 

B01S-B0I 

B1RS-B1R 

B1IS-B1I 

Fl-F-1.000 

BOR  —  (BORS-Y2R-BOIS-Y2II/(F-FI 
BOI  —  (B01S-Y2R-B0RS-Y2II/IF-F) 

B1R  — (81RS-Y2R-B1  IS- Y2I 1/ IF-F1 I 

BlI  — IBM  S-V2R-B1RS-Y2 II/IF-F1I 

S-S-I.OOO/F 

SH-2.O0O-S  -I.  >00 /F 1 

BJOR-BJOR-BOR 

BJOI -BJOI-BOI 

BJIR-BJIR-B1R 

BJI I-BJI I-BI I 

BNOR-BNOR-BOR-S 

BNOI-BNOI-BOI-S 

BN1R-BNIR-BIR-SH 

BNII-BNH-Bll-SN 

IF (8 JOR.NE. TJOR.OR.B JOI .NE.  T  JOI I  CC  TO  B 
IFIBJ1R.NE.TJ1R.0R.BJII.NE.TJII)  GO  TO  7 
IF  (BNOR.NE  •  TNOR.  OR.  BNOI  .NE.TNOI )  GO  TO  6 
IFIBN1R.NE.TN1R.0R.BN1I.NE.TN1I)  GO  TO  5 
GO  TO  10 

5  TN1R-BN1R 
TN1I-BN1I 

6  TNOR-BNOR 
TN01-BN0I 

7  TJ1R-8J1R 
TJ1I-BJ1I 

B  TJOR-BJOR 
T  JOI  -BJOI 

9  F-Fl 

10  BNOR-2. ODO-BNOR 
BNOI -2 . OOO-BNOI 
BN1RS-BN1R 
BN1IS-9N1I 
BJ1RS-BJ1R 
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BJl I S«H J1 I 

BNIR*CHUITRIBN1RS,BNI IS, YR.YI )*CDVDR(  1 . 000,0. ODO, YR, YI I 
BN1 i*CHULTI (BN1RS.BN1IS, YR.YI) *COVDI (  1 .000,0.000, YR.YI ) 
BJ1R*CHULTRI8J1RS,BJ1 IS, YR.YI) 

BJ! I »C MULT I (BJIRS,BJIIS,YR,YII 

6FR-2. 000*1. 577215*669015328  ».5D0*DL0GI VR*VR*VI*YI ) I 
8F  I  «2. 000*0 AT AN2  ( YI  ,  YR ) 

BNOR-ICHULTRIBFR.BFI ,B JOR ,B JO  1 ) -BNOR ) *PI I 
BN0I*(CHUITI(BFR,8FI ,B JOR, B JO I )-8N0 I ) *PI I 
BN1R*(CNULTR|BFR,BFI ,B J1R.BJ1 I )-BNlR)*PI I 
BN1I*(CHUITI IBFR.BFI ,8J1R,8J1I)-8N1 I l*PI I 
11  HOR*BJOR*CNULTRI JR, JI  , BNOR.BHOI  1 
HOI *B JO  I *CHUL1 1 ( JR.JI ,  BNOR.BHOI  ) 

H1R«BJ1R*CHUL TRI JR, J1 .BN1R.BN1 1 1 
HI  I «8 J1 1 *CMULTt I  JR, J I ,BN1R,BN1 1 1 

PHI R«CDVDR IHOR.HOI , I CHUITR ICHUITRI YR , YI ,H1R,H1I ),CNULTI ( YR, YI ,H1R, 
1H1I ) ,2.000,0.0001), ICHULTI ICNULTRI YR.YI ,H1R,H1 I ) .CNUIT It YR, YI ,H1R, 
2H1I 1,2.000,0.000111 

PHI  I -CD VO  I IHOR.HOI, ICRULTRICMULTRI YR.YI , H1R.H1 1 ) ,CMULt 1 1 YR,YI ,H1R, 
1 H 1 1 ) ,2.000,0.000 II, ICHULTI ICNULTRI YR.YI ,H1R,H1I I, CHIU  T IIYR.YI.H1R, 
2H1I I, 2. 000,0.000)1) 

PHI*  CHPL XI PHI R , PHI  I ) 

RETURN 

13  XBR-CMUl TRI1*. 000,0. 000, YR,Yl) 

XBI-CMULTI 116. 000,0. 000, YR.YI ) 

X8QR*CMULTR  I  X8R,  XBI  ,  XBR.X8I  ) 

X80I -CHULTI IX8R.X8I .X8R.XBI ) 

N*2 

IF! (X80R.EQ.0.0). AND. I X60 I.E0.0.0))H*1 
GO  TO  116. 15). H 
16  X8QR-. 10-17 
XBQI *0.000 
15  P0R*1.000 
P0I-0.000 
QOR-POR 
QOI-POI 
P1R-P0R 
P1I-P0I 
01R-3.000 
011*0.000 
TPOR-POR 
TPOI-POI 
TOOR*POR 
TQOI-POl 
TP1R*  POR 
TP1I-P0I 
TQ1R-Q1R 
1011*01 I 
0*1.000 
F*1.000 

X8IR-C0V0RI1. 000, 0.000, X8R, XBI) 

X8I I *CDVD 1 1 1.000 ,0.000, X8R.X8 1 ) 

X8QIR*CDVDR( 1. 000,0. 000, X8QR , XBQI ) 

X80II-C0V0I II. 000,0. 000, XBQR, XBQI) 

00  20  K*1 ,500 
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DN2-0-2.000 
02-0+2.000 
04-0+4.000 
040-04+04 
06-0+6.000 
00-0+0 
020-02+02 
Fl-F+1.000 
F2-F+2.000 
RPO— D2Q+00/IF+F11 
*00— 020+040/ 1 F1+F21 
RPI— 0N2+0+02+D4/1F+F1) 

*01— 0+D2+D4+06/IF1+F2I 

F0*S-( IX8Q1R+RP0I+TP0R-IXBQI I +RPO) +TPO I 1 

POI S- 1  (  X8QI R+RPO) +TPOI +1X8QI I +RPO )  +TPOR ) 

QORS-1 1X8QIR+RQ0I +T00R-1X8Q1 l+ROOI+TQOI I 
001  S-l 1 X8Q1R+RQ0I+TQ0I+IX8Q1 1 +*00 1  +  100*1 
P1RS-1 1X8QIR+RP1 )+TPlR~(X8QI 1+RPll+TPlI I 
P 1 1 S«( ( X801R+RP1 1 +TP1 1 +1X801 1 +RP1 1 +TP1R I 
0 IRS- 1  I X80IR+RQ1 I +T01R- (X801 1 +RQ1 1+T011 I 
Oil S-M  X80 1 * +R01 1+TQ1 1+1X801 1+RQ1I+T01RI 
1FIVAL1)  GO  TO  32 

IF  I IQ1RS+01RS+Q1 1S+Q1 IS ) . LE. (TQ1R+TQ1R+T011 +TQ1 1 ) I  GO  TO  19 
VAL 1-. TRUE. 

32  1FIVAL2I  GO  TO  33 

IF!  (P1RS+P1RS+P1IS+P1  ISl.LE.  ITP1R+TP1R+TP.11+TP11)  I  GO  TO  18 
VAL2-. TRUE. 

33  1FIVAL3I  GO  TO  34 

IF  11 Q0RS+Q0RS+Q01S+Q0 1  SI . LE. ITOOR+TQOR+fOO I+TQOI ) )  GO  TO  17 
VAL3-.TRUE. 

34  IFIVAL41  GO  TO  35 

IF((PORS+PORS+POIS+POIS;.lE.ITPOR+TPOR+TPOI+TPOIII  GO  TO  16 
VAL4-.TAUE. 

35  IFIVAL1.AN0.VAL2.AN0.VAL3.AN0.VAL4)  GO  TO  21 
GO  TO  693 

16  POR-PORS+POR 
POt-POIS+PO! 

IF1TP0R.E0.P0R.AN0. TPOI.EO. POI)  V4L4-.TRUE. 

TP  OR  -P  OR  S 
TPOI-POIS 
GO  TO  35 

17  OOR-QORS+OOR 
001-001 S+QOI 

IFITOOR.EO.OOR.ANO.TOOI.EQ.OOII  VAL3-.TRUE. 

TOOR-OORS 
TOOI  -001  S 
GO  TO  34 

18  PIR-P1RS+P1R 
P1I-P1IS+P1I 

1FITPIR.EO.P1R.ANO.TP1I.EO.P1I)  VAL2-.TRUE. 

TPIR-PIRS 
TPII-Pl IS 
GO  TO  33 

19  Q1R-Q1RS+Q1R 
01I-Q1I S+01I 
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IF  I  TQ1R.EQ. QIR. ANO. TQII.EO.OII)  VALl-.tRUE. 

TQ1R-Q1RS 
TQ1I-Q1IS 
GO  TO  32 
693  0-04 

20  F-F2 

21  QORSi-QOR 
QOI SI-QOI 
OIRSI-QIR 
QIISI-QII 

QOR— CMULTR I QORSt ,Q0 1 S  i ,  X8I R ,  X8 1 1 1 
QOI— CMULTI I00RS1 •Q01S1,X8IR,X8II) 

Q1R«CNULTR(Q1RS1,Q1IS1,XBIR,XBII I 
Oil -CMULTI (Q1RS1,Q1IS1,XBIR,X8III 

PHIR  -COVORICMULTRUR,  JI,  I POR-CNULTR I  JR,  J I ,  OOR,  001 1 1 , 1  POI -CMULTI  I  J 
IR,JI «QOR,QOI III, ICNULTI I  JR, Jl , IPOR+CMULTR I  JR, J I .OOR.OOI II, I POI *CHU 
2LTI I JR.JI .OOR.OOI II 1 1 .CMULTRI2. 000, 0.000,  CMULTR! VR , VI , IP 1R-CMULTR 
31  JR, J I .QIR.QI III, IPII -CMULT II JR.JI, QIR, 01 III 1 ,CMULT 1 1 VR, VI , IP1R-CN 
4ULTRI JR, JI.OIR.QII) I, <PII»CMULT1 1  JR, Jt,0lR,01l III  I, CMULTI 12.000,0 
3. 000, CMULTR I VR ,V| c(PlR -CMULTR! JR, Jl , QIR, 01 1 1 1 , ( PI  I -CMULTI I  JR, J 1,01 
6R, 01 1 1 1 1 (CMULTI I VR, VI , IPIR-CMULTRt JR, Jl ,QiR, 01 1 1 1, t  PI I -CMULTI I  JR, J 
TI .01R.011 1 1 1 1 1 

PHI  I  -CO VO 1 1 CMULTR. JR.JI, (POR-CMULTRI JR, Jl , OOR.OOI I  I , I POI -CMULTI I J 
IR, JI. OOR.OOI 1 1 1, ICNULTI I JR.JI, IPOR-CMULTRI JR.JI, OOR.OOI II, IPCIKMU 
2LTIIJR.JI .OOR.OOI 1 1 1 1. CMULTR 12. 000,0.000,  CMULTR ( VR.VI, I PIR -CMULTR 
31 JR.JI, QIR, 01  III, IPI1-CMULT I (JR.JI, QIR, Qll I  1 1 tCNULT 1 1 VR, VI , I P1R+CM 
4UL TRI JR, Jl, 01R, Oil  1 1, IPII-CNULTI I JR, JI.OIR.QII I II  I, CMULTI 1 2. 000,0 
3.000.CNULTRIVR,VI.IP|R«CNULTRIJR,JI,QIR,OIIII,IPII-CI«ILTIIJR,JI,01 
6R,QIIIII.CNULTI<YR,VI,IPIR»CNULTR|JR,JI,QIR,01III,IP1I«CNULT||JR,J 
71 , QIR, 01 1 1 1 1 1 1 
PHI-  CMPLXl PHI R, PHI  1 1 
RETURN 
END 
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MAIN  PROGRAM  FOR  CALCULATING  AND 
PLOTTING  SERIES  AND  PARALLEL  RES¬ 
ONANCE  FREQUENCIES  WITH  ELECTRON 
TEMPERATURE  AS  A  PARAMETER. 


$1 BFTC  TMPTC1 

COMMON  /INPTS/  T( 10) ,NTfFP( 1 500 1 ,  NF P, R,  AL  ,  S I ( 1500), SSI tO.1500), 
1  SHEATH, FH.FSI 1500) 

DO  6  KKK-1,6 
CALL  INPUT 

CALL  GRIDI6HVREADR,SCALEX,SCALEY,XB,XE,YB,YE,NX,NY,X’5) 

PRINT  3, I  TIKI ,K-1 >10) 

PRINT  *,NT,NFP,NX,NY,(FP(K> ,K*I,5) 

PRINT  3, SCALEX, SCALE  Y«  XB,XE  »>B  «YE X 15, R»AL«  SHEATH, FH 
DO  2  J  -l.NT 

CALL  SICALC I T I Jl ,R, AL.FP.S I ,NFP, SHEATH, FH) 

00  1  K  « 1 ,  NF  P 

FS(K)«S0RT!FH»»2aSI IK)«»2*IFPIK)*»2-FH**2II 

1  SSIJ.KI  -SI  IK) 

CALL  CURVEL!FP,FS,NFP,XB,XE,YB,YE,X15) 

FFF-IFS1NFP)-YB)»9.0/I YE-YB) 

FF-I FPINFPl-XB ) »X15/ I XE-XB) 

AJ-J-1 

CALL  ARROM  IFF.FFF.0.0  I 
CALL  PLOT  I IFF41.0I iFFF.2) 

CALL  PLOTIIFF«2.0),AJ,2  I 
CALL  PLOT  I  I  FF  +  3.0) , A J,2) 

2  CALL  NUMBERUFF+l. 75), AJ..12S, TUI, 0.0,0) 

3  FORMAT! IH  , 1  IE  12 .6) 

A  FORMAT! 1H  , Al 10,5616. fl) 

CALL  OUTPUT 
CALL  FINI 
6  CONTINUE 
CALL  EXIT 
ENO 


SUBROUTINE  FOR  READING  INPUT  DATA 
FOR  MAIN  PROGRAM  TMPTCI 


» I  BFTC  INPUT 

SU6RCUTINE  INPUT 

CCPMCN  /INPTS/  T 1 10)  ,NT,FPI  1500  )  ,  NFP,  R,  AL  ,  S 1 1  1500),  SSI  10,1500, 
1  SHEATM,FH,FS(1500) 

RE  AO  1 ,  NT  ,  NF  P,R  ,  AL ,  F  PMAX  ,  FPM  IN,  S  HEAT  H,  FH 

1  FCRMATI2I5.6E11.5) 

REAO  2,  ITU)  ,J  -1,10) 

2  FCRMATI  5E16.8) 

AN  -  NFP  -1 

OELF  - 1 FPMAX-FPMI N )/ AN 
OC  3  J  -  1 ,NFP 
AJ-J-1 

3  FP I  J)  -FPKl  N»  AJ»OELF 
RETURN 

ENO 
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SUBROUTINE  FOR  PRINTING  DATA 
OUT  OF  MAIN  PROGRAM  TMPTC1 


SUBROUTINE  OUTPUT 

CCMMCN  / INPTS/  T(I0J.NT,FP(I500),NFP,R,ALiSI(i500),SS(l0,l500>, 

1  <FEATH,FH,FS(1500) 

DIPENSION  1(10) 

LCGICAL  FIRST 
FIRST  -.TRUE. 

R  «  R«I.0E*02 
13  KPAGE  •  52 

CC  3  J  *  I.NFP 
KPAGE  •  KPAGE  ♦  1 
IFIKPAGE.GT.50)  GO  TO  A 

1  PPINT  2  .  FPIJIiISSIK.J) «K  -I, NT) 

2  FCRMATI1H  ,FII.3, I0F12.3) 

IFI. NOT. FIRST)  GO  TO  3 

CC  3  K* 1 , NT 

SS(K(J)<>SCRT(  FH»»2»SS(Kt  J)»»2*IFPI  J )  •»2~FH»»2 ) ) 

3  CCNTINUE 

IFI  .NOT  .F  1RST  (RETURN 
FIRST  -.FALSE. 

GC  TO  13 
A  KPAGt-1 
PRINT  5 

5  FCRMATI1HI) 

PRINT  t 

ZIF(FIPST)  PRINT  7 
IFI. NOT .F IRS  T  JPRINT 12 
7  FCAMATI1H  ,  90H 

1C  VL INDR  1C AL  ANTENNA  PS1  CALCULATIONS  I 

12  FCRHATI 1H.105H  CYLINDRICAL  AN 

1TENNA  CALCULATIONS  OF  SERIES  RESONANCE  FREQUENCY  I 

PPINT  6  ,R.  AL.FH 

6  FCPPATI1H  .2AX.I0H  RADIUS  -  .F6.A.17H  CM  ,  LENGTH  «  ,F6.A, 

12BF  METERS  .  GYRO  FREQUENCY  -,F  5.3  ilIH  MEGACYCLES  I 

PRINT  t 
PRINT  9 

9  FORMAT!  92H  PARALLEL  FREQUENCY 
1  ELECTRON  TEMPERATURE  I 

CC  1C  L-1,10 

10  1 (L I  -  TIL) 

PRINT  lI.IKLItL-liIO) 

11  FCRMAT I 16H  (MEGACYCLES)  .  I6,AH*K  ,  91  16,AH»K  )) 

PRINT  t 

GC  TC  1 
ENC 
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SUBROUTINE  TO  CALCULATE  SERIES 
RESONANCE  FREQUENCY 


ilBFTC  S 1  CALC 

SUBROUTINE  SICALC(T,R,AL,FP, SI, NFP, SHEATH, FM  I 
DIMENSION  FP(ISOO)  ,  SI (1500 > 

COMMON  /CONSTA/  PI , PI  2 ,EMASS , AMASS , O.CK, EO, G,UO, ALP MAK 
HEAL  K 

CL  •  ALOCIAL/R ) 

CALL  CCNSTS 

OBYCON  »SQRT(CK/(4.0»P12»EMASSI )»1.0E-06 
S12  ». I 

DC  7  J  -l.NFP 
L  -  NFP+l-J 

IFIFP(L).LE.FH)  CO  TO  10 
ALPFT  »ALPHAK»IFP(L) •»2-FH»»2)/T 
DBYLNT  -  08VC0N* SORT ( T l/FPI L) 

RS  •  M  08YLNT»SHEATH 
SHTCON  »  ALOGIRS/R) 

N  -  0 
GO  TO  3 

1  SI2«AN».04 

2  SI2*.5»SI2 

3  S2  -1.0  -SI2 
N  -  N*l 

AN  *  N 

IF(N.GT.IOO)  GO  TO  8 
A2  -  S2»ALPFT 

IF ( A2. LE . l.OE-18 )  A2  »  1.0E-18 
AR  «  R S» SORT ( A2 ) 

RKAR  «  MO.ARI/Kll.AR) 

KKM  -  1 • 0-RKAR»»2 
S72  -  2.0»S2 

SI2P  •  ( S12»RKM*S22» I SHTCON tRKAR/AR 1 1 / 1 RKM*S22*CL ) 
SLOPE  -  1.0  ♦  RKM/(S22*CL) 

DEL  -IS12P  -SI2) 

IF  (SLOPE  1 2, 2, 4 

4  IFISI2P  12,2,5 

5  IF ! SI2P-  1.016,1,1 

6  SI2  -  SI2P 

IFIA6SI0EL1 .GT..0011G0  TO  3 

7  SIILI  *SQRT( S 1 2 1 
RETURN 

B  PRINT  9,UEL,SI2,AR,T,FPtJ) , SLOPE 

9  FORMAT  I 3IH  SICALC  HAS  FAILEO  TO  COVERGE.  ,6E16.8I 
GO  TO  7 
10  SI2-0.0 
GO  TO  7 
FND 
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SUB  FUNCTION  FNOQftAM  FOB  CALCULATION 
OF  THE  BCSSCL  FUNCTIONi  K0  AND  K, 

REAL  FUNCTION  K(N.X) 

REAL  I 

K  -  0.0 

IF(N.LT.O)  N  — K 
IFCN.GE.2 I  RETURN 
IFlX.eO.O.OUO  TO  1 
IF(N.EQ.l)  60  TO  3 
IFU.GT.2.0)  SO  70  2 

V  -  X/2.0 
T2  -  V*Y 
Y4  ■  Y2*V2 
Y*  •  Y4*Y2 
VI .•  Y4*V4 
Y10  -  V6*V4 
¥12  ■  V6*Y6 

K  -  -ALOCI  Y J*I(N* XI-.97721966  >.42278420*92  t.23069796*Y4 
14  .03488990*96  ♦  .00262698*V8  >.00010790*Y10  4  .^0000740*Y12 
RETURN 

1  K»1.0E«37 
RETURN 

2  Z  -  2.0/X 
12  •  Z*Z 
Z  3  -  Z2*Z 
Z4  -  Z3*Z 
Z  3  -  Z4*Z 
Z6  -  Z9*Z 

K  -  1.29331414  -.07032398*2  ♦  .02189968*22  -  .01062446*23 
14  .008878  72*  Z4  -.00291940*29  *.0009320I*Z6 
K  -  K*EXF(-X1*(X**(-.91 1 
RETURN 

3  IFIX.ST .2.0)60  TO  4 

Y  •  X/2.0 
Y2«  Y*Y 
Y4«  Y2*Y2 
Y6-  Y4*Y2 
Y9-  Y6*Y2 
910-98*92 
Y12*V10*92 

R  -  ALOC(YI*I(N.XI  4(1.0  4. 19443144*92  -. 67278979*94- .18194097*96 
1-  . 01919402* Y8  -  .00110404*910  -  . 0000468 • VI 2 >/X 
RETURN 

4  Z  -  2.0/X 
Z2*  Z*Z 
Z3*  Z2*Z 
Z4*  Z3*Z 
Z9-  Z4*Z 
Z6-  Z9*Z 

K  *  1.29331414  4.23498619*Z  -  .03699620*Z2  4  .01904268*23 
1  -  .00780393*Z4  4  .00329614*Z9  -  .00068249*26 
X  -  R*EXF (-X1*(X*»(-.9I  I 
RETURN 
END 


1-38 


SUBFUNCTION  PROGRAM  FOR  THE  CALCULATION 
OF  THE  BESSEL  FUNCTIONS  Iq  AND  I, 

REAL  FUNCTION  IIN.X1 

OCUSLE  PRECISION  Y, FACN, FACN.FJ, SUM, YEAN.FK.FKN 
V-X/2.0 
I  -C.0 

IF  I  N.LT .0 )  N  •  -N 
FACN  -.1.0 
FACN  -1.0 
IFIN.EQ.O I  CO  TO  4 
DC  U  *  ltN 
FJ  -J 

1  FACN  -  FACN-FJ 
SLR  -{YaaNJ/FACN 
TERN  *  SUN 

2  00  )  K  >1  i  SOO 
FK  -  K 

FKN  •  K  ♦  N 

TERN  -  TERM-Y-Y/I FKN-FK) 

SUN  -  SUN  ♦  TERN 
SS-SUM 

IFI I.EQ.SS  )  RETURN 

3  I  -  SS 
RETURN 

4  SUN  a  1.0 
TERN  -  i.o 
CO  TO  2 
END 


SUBROUTINE  FOR  DRAWING  A  CURVE 
FROM  AN  ARRAY  OF  DATA 


AIBFTC  CURVEL 

SUBROUTINE  CUR VEl (A* ,AY ,N,X8 ,XE, VB, YE ,X1SI 
OINENS ION  AXIIOI.AYIIOI 
LOGICAL  FIRST 
FIRST  -  .TRUE. 

OX  -XE-XB 
OY  »VE-YB 
SCX-XI5/DX 
SCV-9.0/0V 
DO  l  J  -l.N 
X  -IAXIJI-  XBI*SCX 
Y  alAYIJI-  YB ) »SCY 
IFIX.GT.XIS)  CO  TO  3 

IFIX. LT.O.OI  GO  TO  3 
IFIV.GT.9.0I  CO  TO  3 

IFI Y . L T .0.0)  CO  TO  3 

IFI  F  IRST  (CALL  PL0T(X,Y,3I 
CALL  PLOTI  X,  V,  21 
FIRST  -  .FALSE. 

GO  TO  1 

3  FIRST  ■  .TRUE, 
t  CONTINUE 
RETURN 
ENO 
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SUBROUTINE  TO  PROOUCE  A  LABELED 
QRIO  FOR  PLOTTED  DATA 

SUBROUTINE  CR ID( T ITLE ,SCX,SCY,XXB,XXE, VYB ,VYE,NXX,NVY , XX15J 
C  ALTERNATE  ARCUNENTS  CRI0I6H  RE AO  )  •  CR ID! 6HVRE ADR)  ,  CRIDI6KCONNON) 

DIPENSICN  T ITLEI 52 ) 

COPHON/PLTCHN/HEADI 65), SCALEX, SCALEY, XB,XE,VB, YE, NX, NY, X15 
CALL  HOLt  CHN«6HC0NM0N ) 

I F ( TITLE! I)  .EO.CNN)  CO  TO  106 
CALL  NOLI READIN,6H  READ  ) 

CALL  H0LIREADER.6HVREA0R) 

IFIITITLEIU.NE. READER). AND. ITI TLEI1) .NE. READIN ) I  CO  TO  104 
READ  I01«IHEADI J),J  -1,52) 

101  FCRHATI 13A6) 

READ  102,  SCALEX,  SC  ALEV,  XB,XE  ,VB« YE,NX,NY,X15 

102  FCRHATI  6F10. 5,215, 1F10.5J 
PRINT  1000 

1CC0  FCRHATI  1H1) 

PRINT  1C2,SCALEX  ,  SCALE Y,XB,XE,V8,VE,NX, NY.X15 
PRINT  1COO 

IFITITLEI1I.EQ.REA0IN)  CO  TO  103 

SCX  «  SCALEX 

SCY  -  SCALEV 

XXB  «  XB 

XXE  -  XE 

VYB  -  VB 

YYE  -  YE 

NXX  -  NX 

NYY-  NY 

XX15  -  X15 

103  CALL  IDPLOT 
GC  TO  106 

104  CALL  IDPLOT 

CC  105  J  -  1,52 

105  HEADIJ)  -  TITLEIJI 
SCALEX  -SCX 
SCALEY  -SCY 

XB  -  XXB 
XE  -  XXE 
YB  -  YYB 
YE  -  YYE 
NX  -  NXX 
NY  -  NYY 
X15-  XX15 

106  CCNTINUE 

IFISCALEX.LT. 0.0)  SCALEX  -  SCALEX 
I F I SCALEX .GT . 100.01  RETURN 
IFISCALEV.LT. 0.0)  SCALEY  -  -SCALEY 
IFISCAIEY.GT. 100.0)  RETURN 
IFIX15.LT .0. 0)  X15  -  -X15 
IFIX15.GT. 100.0)  RETURN 
X  -  0.0 
Y  -  0,0 
AN  -  N 

CALL  PLOT  I  4.0,1.45 ,-3) 

12  CALL  PLOT  10.0  ,9.0,  2) 

CALL  PL0TIX15  ,4.0,  2) 

CALL  PL0TIX15  ,0.0,  2) 

CALL  PLOT  10. 0,  0.0,  2) 

CX-IXE  -XBl/SCALEX 
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i 


2 

3 


4 

3 


* 

7 


C 


199 

2C0 

201 

210 

211 

212 


BY- 1  YE  -  VBI/SCALEV 
AX  -  X13/SCALEX 
AX10  -AX/10.0 
AY  -  9.0/SCALEY 
AY10  -AY/10.0 
V-V«AY10 

IFIV.6E.  9.0  I  60  TO  2 
CALL  nOTIO.OtY.21 
CALL  nOTIO.l.V.21 
CALL  nOTtO.O.V.3} 

60  TO  1 

Y  -  9.0 

CALL  n0T(0.0t9.0t21 
X  •  XAAX10 

IFIX.6E.X19  1  60  TO  4 
CALL  nOTlX  ,9.0.21 
CALL  nOTtX  ,0.9,21 
CALL  OLOTIX  ,9.0,31 
60  TO  3 

CALL  n0T(X15  ,9.0,21 

Y  -Y  -AVIO 
IFIV.LE.O.OI  60  TO  * 
CALL  nOTlXIS,  Y,2  I 
CALL  FLOTIXll-.l.Y.B  1 
CALL  nOTIXIS  ,Y,3  1 
60  TO  9 

X  -X19 

CALL  FL0TIX19  ,0.0,21 
X  -  X  -AX  10 
IFIX.LE.0.01  60  TO  6 
CALL  FLOT  IX, 0.0,21 
CALL  FLOT  IX, 0.1, 21 
CALL  FLOT  IX, 0.0,31 
60  TO  T 
X  -  0.0 

Y  -  0.0 

XS  -  1X19-13.41/2.0 
CALL  nOTIO.O,  0.0,  21 
NNX  -  0 
NNY  >  0 

I  FI  XE.1 200  ,200,201 
IF( YE  1210,210,211 
XA  -  ABSIXB1 
GC  TO  199 
XA-XE 


60  TO  199 
YA  -  ABSIVBI 
GC  TO  212 
YA  -  YE 
CONTINUE 
IFIXA.LT.  1000.1 
IFIYA.LT. 1000. I 
IFIXA.LT.  100.1 
IFIYA.LT.  100.1 
IFIXA.LT.  10.1 
IFIYA.LT.  10.1 
IFIXA.LT.  1,1 
IFIYA.LT.  1.1 


NNX  •  NX  - 
NNY  •  NY  - 
NNX  •  NX  - 
NNY  •  NY  - 
NNX  •  NX  - 
NNY  •  NY  - 
1'NX  -  NX 
NNY  -  NY 


IFINNX.LT .0)  KMX  -  0 
1F1NNV.LT .0)  NNY  ■  0 
XSN  •<  NX 
VSN  -  NY 

SPACEX  -  .07+XSN 

SPACEY  -  .07+VSN 

1 F (  XE.LT.0.0)  SPACE X  »  SPACEX  ♦  .0T 
IF1YE  . LT.O. 0) SPACE Y  »  SPACEV  ♦  .07 

xn  -  xe 

VN  •  YB 

CALL  NUMBER! X-SPACEX.-.25.. 15.XN,  O.O.NNX) 

CALL  NUN8ER!-.10.V-SPACEV..13.VN.90.0,NNV) 

9  V  -  Y  PAY 
YN  •  VN+DV 

IF! Y .ST .9.0)  60  TO  10 

CALL  NUMBER!-. 10, Y-SPACEV,.1S.VN. 90.0. NNY) 
CALL  PLOT (0. 0.V.3) 

CALL  PLOT  1X15  .V.2I 
V  -  V  ♦  AY 
YN  •  VN  ♦  OY 
IFIV.6T.9.0I  60  TO  10 
CALL  PL0TIX1S.V.3) 

CALL  PLOT I0.0.V.2) 

CALL  NUMBER!-. 10, V-SPACEV,. 15. YN. 90.0, NNY) 

CC  TO  9 

10  X  ■  X  PAX 
XN  »  XN+DX 

IFIX.CT.X15  »  60  TO  11 

CALL  NUMBER! X-SPACEX, -. 25, .15.XN.0.0.NNX) 

CALL  PLOT IX. 0.0.31 
CALL  PLCT1X, 9.0.21 
X  •  X  ♦  AX 
XN  a  XN  +  OX 
IFIX.GT.X1S  1  60  TO  11 
CALL  PLCTIX. 9.0.31 
CALL  PLOTIX. 0.0.21 

CALL  NUMBER! X-SPACEX,-.25».15»XN,0.0.NNX> 

CC  TO  10 

11  CALL  SVMBL4I 1 . 65+ US.*. 55. .15, HE AD  1 1 1.0.0. 78) 
CALL  SV+5L9I XS.  -1.00. .20  .HEADI27). 0.0.78) 
CALL  SVPBL4I 1.43+XS,-1.30,.13,HEA0I40),0.0,78) 
CALL  SVMBL4I -.35,-. 35. . 15  .HEAOI 14) ,90.0. 78) 
RETURN 

ENC 


MAIN  PROGRAM  FOR  CALCULATION  OF 
TEMPERATURE  FROM  SERIES  ANO 

parallel  resonance  frequencies 


S-A2 


SIBFTC  TEP601 

DIMENSION  FMTIIA) 

IS  FURMAT(12X.2F8.2.F10.2,2F<).0,10X.I6> 

COMMON  /KEN/  R,  XL.  PERCNT,  LIMIT.  CC.  S,  TSTART. F.FP.T.KK.N 

1  FORMAT  I3FI0.2,  15.  E15.7.  2F10.2I 
RtAO  1.  R. XL. PERCNT.LINIT.CC. S. TSTART 
PR  INTI,  R. XL, PERCNT. LI MIT.CC .S .TSTART 
X- 1 

READ  2.  FHT 
PMINT2.  FMT 

2  FORMAT  1 13A6.  »2I 
CALL  CROIDS 
PRINT  3 

3  FORMAT  ( 1H1 ,  SX.  AKTIHE,  2X,  6HSERIES,  5X,  3MPAR,  5X,  MSEC,  2X, 
1  8HNEM  TEMP,  2X,  8H0L0  TEMP,  5X,  SHALT) 

8  READ  FMT,  F,  FP,  SEC.  OT 
CALL  TEMPSP 

PRINT  15.  F.FP.SEC.T.OT.K 

PUNCH  IS,  F.FP.SEC.T.OT.K 

X-  K*1 

GO  TO  8 

CALL  CRD1DF 

CALL  EXIT 

END 


SUBROUTINE  FOR  CALCULATION 
OF  THE  ETA  FUNCTION 


fIBFTC  ETAA 

SUBROUTINE  ETAAI  Y  .  ETA.  XL) 

REAL  K 
N-0 

AJ  -0.0 
AL  -0.0 
XL-  0 

IF  IY.GT. 0.0)60  TO  1 
ETA  -1.0EO7 
RETURN 

1  X-1.0 
DX  -  .5 

DYH  ■  V  4.00001 
DYL  -  Y  -.00001 

2  I  FI  X.LE  .0  .0 )  X  -1.0E-37 
XL-  KIM 

I  FIX  .GT.  I.OEOO)  GO  TO  5 
IF  (N.GT.SOO)  GO  TO  5 
PHI  -  X  10  ,X )  /  ( X*K( 1  ,X) ) 

I F I  PHI .GE .DYH )  CO  TO  3 
IFIPHI.LE.DYL)  GO  TO  A 
ETA  -  X 
RETURN 

3  AL  -  AL  *1.0 
AJ  >  0.0 

OX  -  OX  *AL/2.0 
X  -  X*DX 
GG  TO  2 

A  AJ  -  AJ  *1.0 
AL  -  0.0 

nx  -  dx»aj/2.o 

X  -  A  -OX 
GO  TO  2 

5  PRINT  6  .  X  ,V  ,DX  .PHI 

6  FORMAT! 1«H  ETA  HAS  FAILEO  .AE20.B) 
ETA  ■  X 

RETURN 

END 
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SUBROUTINE  FOR  CALCULATING  ELECTRON  TEMPEATURE 
FROM  SERIES  AND  PARALLEL  RESONANCE  FREQUENCIES 

SIBFTC  TENPSP 

SUBROUTINE  TENPSP 

CON  NON  /K«/  R.  At*  PERCNT,  UNIT,  CC,  S,  1ST  ART,  F,  FP,  T,  KK.  N 
DATA  P 1/3 .14 159265/,  EN/9.1083E-31/,  EK/1 .38044E-23/,  6/3.0/, 

I  A/0.0/,  TH/1.4/,  EOO/I. I5434E-12/ 

PI2-  PI*PI 

TCON*t EN/(6»EK) )»4.0£412»PI»PI 
0EVCCN-(EK/(4.0>  PI2*EH>> 

OETCCN  -  SORT  ( DBVCON)  »1.  OE-06 
CKCN»  0.0 

IFICC  .6£.  1.0)  60  TO  22 
CK0N-(2.0aPI«AL/CC)»E00 
22  FH2-  FM»FM 
N»0 
KK-0 

FP2  -  FP*FP 
F2  -  F»F 
FN2-  FP2-FH2 
FC2  -  FP2-F2 

IF  (FN2  .IE.  O.Ot  60  TO  2B 
IF(  F  .IT.  FH)  60  TO  2B 
FN-  S0RTIFN2) 

SI2*  CF2-FH2  J/FN2 
S2-  1.0-SI2 
T-  TSTART 
32  TEST  -  T 

0LEN6T-  OBVCONaSORT (TI/FN 
RS«  R*OLENCT»S 

PHI-  SI2-AL0CI AL/RS )  -  S2»C AL061RS/R)  ♦  CRON) 

IF  (PHI  .IE.  0.0)  60  TO  28 
CALI.  ETAAIPHI,  ETAS,  III 
KK-KK4KL 

T-  TC0N-RS-RS-FC2  /CETAS*E7.«) 

N-N*l 

IF  (N  .6T .  UNIT)  60  TO  36 
IF  (S  .EQ.  0.0)  RETURN 

IF  ( ( ABS( T-TESTl/T)  .6T.  PERCNT)  60  TO  32 
28  T«  0.0 
RETURN 

36  PRINT  37,  F,  FP,  T,  N,  KK 

37  FCRHAT  UX1,  21HT00  RANT  ITERATIONS  (,  2F3.2,  FIO.O,  217,  IN)) 
RETURN 

END 
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MAM  MOQRAM 
ANTENNA  REMTANGC 


KM  PLOTTM8 
VERSUS  REACTANCE 


8I8FTC  IttHT  MAIN  PROGRAM  FOR  PLOTTING  ANT  Z  VS  R  AND  X 
C  PROGRAM  FOR  PLOTTING  ANTENNA  1NFE DANCE  AS  A  FUNCTION  OF  R  AND  X 
OINENSION  HE AO! 50 », TEN  (HI 
CONFLEX  Z.ZO 
LOGICAL  FIRST 

COMMON/ IHPCON/Z.F.FP.FNU.FH. ANGLE .T.R.AL 

COMMON  /CONSTA/FI ,PI 2 .EMASS. AMASS.Q.CX.EO.C,  UO,  ALPHA* 

CALL  CONSTS 
alfh-alfhak 

READ  I  • NPLO  TS ,  NRESL  V  »  NTAFE • NSK I  Ft  NFRI NT  * NFLT *  NTHY 
FRINT1  .NPLDTS.NRESLV.NTAPE.NSKIP.NPRINT.NPLT.HTHV 

1  FORNATI 71101 

00  32  JNF  ■  l.NPLOTS 
READ  10  .NOO.KCRID.FMAX.FX 
FRINT  3  .MOD, KCRID.FHAX.FX 
3  FORNATI 1H1.21 I0.2F10. 3 1 
10  F0RNATI2110.2F10.3t 
FX*FX*. 1E*0T 
FMAX*FNAX*. IE* 07 
REAO  14, FH.FNU. T.R.AL. ANGLE 
FRINT  2, FH.FNU, T.R.AL, ANGLE 

2  FORNATI 1H1.6E1 3. 41 
14  FORNATI 6£ 13.61 

20  FORNATI2E20.8.F10.3I 

21  FORNATI 1H1.2E20. 61 

REAO  20.FNUC0N  ,  DEN, ALT 
FRINT  21  ,  FNUCON  .  OEN 

CALL  SR10I6NVREA0R, SX.SY.XB, XE.YB.YE.HX.NY.XISI 
SCX  -  X15/IXE-XBI 
SCT  -  9.0/1 VE-YBI 

C  FLOT  FOR  10  CONSTANT  TENFERATURES 

OATA  TEN/250., 350., 500., TOO., 1000., 1400., 2000., 2600. .4000.. 5600. 

1  ,  4000. 

F  •  FS 
FF  -  0.0 
T  -  0.0 
CALL  ANT2 
ZO  -  2 

RES  -NRESLV 
REF  -FNAX/RES 
FRINT  1000 
1000  FORNATI 1H1 I 

FRINT70, SCX  ,SCT  .F.FP.T.Z.XES.REF 
TO  FORNATI 9E 14. 6) 

FRINT  1000 
0030JT  -1,11 

89  FCRNATI 127H  FLOT  X  FLOT  Y  OF  FREQ  FLA  FRE 

10  TENFERATURE  RESISTANCE  REACTANCE  FUN 

2  REF  //.19H  TENFERATURE  EQUALS.F20.1.8H  DEGREES. . 12HALT.  EQUALS. 

3  F10.2,  3H  KN///I 
FRINT  86 

88  FORNAT  (1H1) 

FRINT  89.T.ALT 

IF1NTHY.EQ.1)  FNU  -  FNUCDN*DEN»T 
IFIFNU.GT.2,*FI«FXI  ALFHAK*ALPH«5. 0/9.0 
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PXO  -  0.0 
PVO  -  0.0 
FIRST  -.TRUE. 

JL  -  0 

00  28  JF-l.NRESLV 

FJ  •  JF-1 

FF  -  FJ-REF 

CALL  ANTZ 

Z  -  Z  -  ZO 

PXO  -  PX 

PYO  -  PV 

PV  -(REAL! ZI-YB) *SCY 
PX  -U(NAG(Z)-XB)*SCX 

IF  I NPRI NT . E 0.0. AND. JF/5-5.EQ.JF)  PRINT  70.PX,PY,F,FP.T,Z.FNU,REF 

IFIPX.GT.X15)  PX  -  X 15 

IFIPX.LT. 0.0)  PX  -  0.0 

IFIPV.LT.O.OIPY-O.O 

IFIPT .GT.9.0)PV*9*0 

IFIFIRSTICALL  PL0TIPX.PY.3) 

FIRST  -.FALSE. 

CALL  PL0T(PX,PY.2I 
JL  -  JL  *1 

IF  I JL. LE.  MPLT)  GO  TO  28 

JL  -  0 

OXX  -  PX-PXO 

OYY  •  PV— PYO 

100  IF (OXX I  103.101.103 

101  IFCDYV)  102. 20. 102 

102  OXX  -  .1E-0A 

103  S  ■  OYY/OXX 
S2  •  S-S 

OV  -.05- SORT ( S2  /  ( S2+1 .0) ) 

OX  -.05*SQRT( 1.0/IS2  >1.01 3 

IF(S)  104.105.105 

104  OX  -  -DX 

105  CALL  PL0T(PX-DV.PY»DX.3I 
CALL  PLOT  I  PXpDY.PY-0X.2I 
CALL  PLOT (  PX,  PY,  31 

28  CONTINUE 

29  T  -TEN(JT) 

IFlALPHAK  .NE.ALPH)  ALPHAK-ALPH 

30  CONTINUE 

32  CALL  PL0TCX15P4. 0.-1. 45,-3) 

CALL  F INI 

STOP 

END 
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